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ABSTRACT 


The objective of this research is to investigate the dynamic elastic properties of materials and 
their temperature and frequency dependence using an acoustic resonance based method. In this 
technique, the torsional, flexural, and/or longitudinal vibrational modes of a “free-free" bar are 
selectively excited and tracked as a function of temperature. The resonance frequency of the torsional 
mode is used to obtain the dynamic shear modulus. The dynamic Young’s modulus is obtained from 
either the resonance frequency of the flexural or longitudinal mode of the bar. The quality factor, 
Q, of each mode is measured to obtain the damping properties of the material. A two channel phase- 
locked loop (PLL) is used to track the resonance frequency as a function of changing temperature for 
the particular resonance mode selected. Using this technique, the storage modulus and loss tangent 
may be obtained in a continuous fashion. Materials tested in this thesis include: Polyurethane PR- 
1592, a common sonar encapsulant, Polymethy] Methacrylate (PMMA) or plexiglass, and 


Polycarbonate. This research encompasses the theory, accuracy, limitations, and applications of this 


measurement technique. 
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I. INTRODUCTION 


A. BACKGROUND AND MOTIVATION 

Elastomeric materials are important in many engineering 
applications. This research, in part, was undertaken to 
support continuing efforts to develop fiber optic acoustic 
sensors. The design and fabrication of fiber optic 
hydrophones is particularly dependent upon the mechanical 
properties of the materials used in the hydrophone. In 
order to enhance the sensitivity of fiber optic devices, 
various compliant transducers have been designed to generate 
large amounts of strain or phase shift within the arms of an 
optical fiber interferometer [Ref. 1,2]. The transducers 
have been produced in the shapes of plates, mandrels and 
ellipsoids. These techniques often involve bonding the 
fibers to an elastomeric material in such a way as to induce 
strain within the fiber itself. More common interest in the 
acoustic properties of materials originates from the 
automotive and industrial communities in which the 
elastomers play an important role in the damping of noise, 
shock and vibrations. 

It is the static elastic moduli that determines the 
deformation of a body to a constant force. For example, in 


a hydrophone the pressure at a particular operating depth 








produces a force on the transducer that results ina 
corresponding deformation. However, the dynamic moduli are 
used to determine the acoustic sensitivity or responsivity 
of a transducer to time-varying loads. 

The technique used here for measuring the dynamic 
Young's and shear modulus is a refinement of one first 
developed by Barone and Giacomini [Ref. 3] with 
modifications made by Rudnik at UCLA. Additional revisions 
to the technique and apparatus have been made by Garrett and 
Brown at the Navy Postgraduate School [Ref. 4,5]. This 
method was recently used by Wetterskog, et. al, [Ref. 6] to 
measure the temperature dependence and dynamic moduli of 
materials considered for use in their fiber-optic 
hydrophones. The technique was also used by Tan [Ref. 7] 
while investigating the viscoelastic transitions of 


elastomeric materials. 


B. METHOD OF MEASUREMENT 

The dynamic moduli of materials are often measured with 
forced vibration techniques which are generally classified 
as either resonant or non-resonant procedures. Here, a 
resonance-based method employing a freely supported rod is 
used. The resonant frequency is proportional to the square 
root of the modulus of the material. In general, the 
displacement of a sample at resonance is larger by a factor 


equal to the quality factor (Q) than it is under a static 


load. Accordingly, when the Q>1, resonant techniques have a 
higher signal-to-noise ratio than nonresonant techniques. 
The quality factor is the ratio of the energy stored to the 
energy lost within a material during one cycle and is 
inversely proportional to the loss tangent. The loss 
tangent, tan 6, is defined as the tangent of the phase angle 
between the applied stress and resulting strain. The Q, 
therefore, is a measure of the absorptivity of the material. 
The complex dynamic moduli (storage and loss properties of 
the material) can be determined using the resonant frequency 
of each mode, the physical dimensions and the mass of the 
sample rod. 

The method discussed here employs a forced vibration 
created by electrodynamic transduction. The force is 
generated on a small wire transducer coil mounted ona 
sample rod placed between a pair of permanent magnets. The 
detection of the vibration at the receiver is also by means 
of electrodynamic transduction. A receiver coil is mounted 
on the opposite end of the sample, and located between an 
identical pair of permanent magnets. An important feature 
of this technique is that the torsional, flexural, or 
longitudinal mode of the sample rod can be selectively 
excited depending upon the orientation of the transducer 
coil within the magnetic field. Thus, both the shear and 
Young's modulus can be determined independently using this 


single apparatus. The resonant technique utilized in this 





investigation is discussed and illustrated in more detail in 


References 4 and 5 as well as Chapter III of this Thesis. 


Cc. SCOPE 

The scope of this project is to further refine this 
measurement technique to permit independent measurement of Q 
as a function of temperature, and develop new methods to 
test very compliant and high loss materials. Further, the 
capabilities and limitations of this method as to 
restrictions on damping and stiffness properties of 
candidate materials were to be investigated. Materials 
tested in this thesis include; Polyurethane PR-1592, PMMA, 


and Polycarbonate. 


II. THEORY 


A. INTRODUCTION 
1. Classical Theory 
The classical theory of elasticity describes the 
behavior of materials categorically as either ideal elastic 
solids or viscous fluids. The mechanical properties of the 
elastic solid are described by Hooke's law: For small 
deformations, stress is directly proportional to strain and 


independent of the rate of strain. Algebraically that is, 
o«<6€, O0=Ge; (1) 


where o is the applied stress, ¢ is the resulting strain, 
and G represents the appropriate modulus for the system (in 
this case G is the shear modulus). 

This implies that the moduli of elastic materials 
have no frequency dependence. Ideal viscous fluids are 
governed by Newton's viscosity law: The applied stress is 
proportional to the rate of strain and independent of the 


strain itself. Algebraically that is, 


vo. av 


on’ = ‘On i (2) 


ox 





where o is the applied stress, @V/dn is the rate of strain 
and p is the coefficient of viscosity. The energy 


transmitted to a viscous liquid is dispersed as heat. 


Contrastingly, in an ideal elastic solid the energy of 
deformation is entirely recoverable. [Ref. 8,9,10] 

Under certain stress and temperature conditions many 
materials exhibit a combination of elastic and viscous 
behavior. This suggests that materials may exhibit either 
ideal Newtonian viscosity or ideal elasticity properties in 
certain temperature regions. These materials are classified 
as being viscoelastic. The mechanical properties of these 
materials are specified by a frequency and temperature 
dependence of the complex elastic moduli and loss modulus. 

2. Response of a Driven Harmonic Oscillator 
If a harmonic oscillator is driven with a sinusoidal 


force such that, 
F(t) = Fsin(wt) , (3) 
then the displacement for the system is of the form, 
x(t) =X sin(wt-o) . (4) 


g@ is the phase difference with respect to the driving force. 
The amplitude of displacement (X) in terms of the resonance 


frequency and the quality factor can be expressed as, 


eS =e ee, 


(5) 
K4| (1-£,932)? + (222 


where f,=f/f, (f, is the system resonance frequency}. The 


dimensionless normalized displacement response of the system 


is therefore expressed as, 


Xx 1 


F/K : 
/ 2)2 + ( Zze1)2 
The phase difference between the driving force and 


displacement is, 


(6) 


(7) 


Equations 6 and 7 are plotted verses the frequency ratio for 


various Q's in Figures 2.1 (a) and (b). 
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Figure 2.1 (a) Normalized displacement response vs. 


f.. for Q's of 10, 5, 2, 1 and .5. 
(b) Phase (displacement) vs. f 
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Note that the maximum displacement of the system 
does not necessarily correspond with the resonance frequency 
(f,,=1). This is especially true for low Q systems (Q< 5). 
However, it is seen from figure 2.1(b) that the maximum 
rate-of-change of the phase always occurs at resonance, 
which corresponds to ¢~=90° for the displacement response. 

The corresponding velocity of the system is obtained 
by taking the derivative of the displacement. The 
amplitude of the velocity response (v) is, 


F/K . 
iE, (8) 


(1+£,,.7)" + — L* 
The normalized amplitude of the velocity response 


can be expressed as, 


Vege 
F/K F (9) 
aed 2 re ( rel ) 2 


The phase of the velocity response then becomes, 


).= 90" (10) 


Equations 9 and 10 are plotted verses the frequency 


ratio for various Q's in Figures 2.2 (a) and (b). 
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It can be seen from these two plots that the maximum 
velocity amplitude does occur at resonance (f=f, or f,,=1), 
no matter what the Q for the system. Further, the maximum 
rate-of-change of the phase always occurs at resonance which 


corresponds to ¢=0° for the velocity response. 


B. DYNAMIC MODULI 
1. Complex Moduli 

The technique used in this thesis determines both 
the shear and Young's complex dynamic moduli by measuring 
the resonance frequency and Q of the fundamental modes of 
vibration of a rod. If the dimensions of the rod are known, 
the shear modulus can be calculated from the speed of 
torsional waves in the rod and the Young's modulus can be 
calculated using the speed of either the longitudinal or 
flexural waves in the rod. 

If a periodic shearing stress is applied to an 
elastomeric material, then a periodic strain having the same 
frequency will develop in the material. Figure 2.3 isa 
sketch of an oscillatory stress (0) and the resulting strain 
(e€), between which there exists a phase angle (6). The 
phase angle is between 0 and 90° and is dependent upon the 
amount of damping present in the material. 

The complex shear modulus (G’), which relates the 


complex stress and strain, can be represented by two vector 


components [Ref. 8,9] and is illustrated in Figure 2.2. 





Gt=G'+G" . (11) 


The storage or elastic shear modulus, G', is defined 
as the ratio between the stress and the in-phase component 


of the strain [Ref. 8]. The storage modulus is proportional 


to the energy stored during one cycle of the applied stress. 





Figure 2.3 Time profile of applied periodic shear 
stress and related strain. 

The viscous or loss shear modulus (G") is defined as 
the ratio between the applied stress component and the 
quadrature component of the strain [Ref. 9,10]. The loss 
shear modulus is a measure of the energy dissipated or lost 
as heat during one cycle. The storage and loss modulus are 


related by the tangent of the phase angle (6) as shown in 
Figure 2.4. 
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G* 
G" 


Figure 2.4 Vector diagram of the complex shear 
modulus. 


If a sinusoidal tensile stress is applied to the 
elastomeric sample instead of a shear stress, then the 
resulting longitudinal strain will be out of phase with the 
applied tensile stress. The resulting complex Young's 


modulus (E’) is related to its storage and loss components 


by, 
E*=E/+E! . (12) 


The complex modulus can be explained using a forced 


damped harmonic oscillator model as shown in Figure 2.5. 








Figure 2.5 Damped simple harmonic oscillator with 
applied force. 


The equation of motion for this system is 


P(t) = MoX() 4p OX) ae xcey (13) 


Ft ™ dt 
where X is the displacement from equilibrium, M is the mass 
of the oscillator, R, is the mechanical resistance, K is the 
spring constant, and F(t) is the applied force [Ref. 11]. 
If the force applied and resulting motion is sinusoidal, 
then the expression, after appropriate differentiation, can 


be written in complex variable form, [Ref. 5] 


F = -w?MX* + jwRX* + KX’ (14) 


’ 


or 


F = -w?Mx* + K(1 + 9%) x . (15) 
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It can be seen from Equation (15), that the 
coefficient of the second term is the effective stiffness 
for the system. This effective stiffness is analogous to 


the complex modulus and can be expressed as follows, 


K*=K(1 + 54) = K(1 + jtans) . (16) 
When the driving frequency coincides with the 
mechanical resonance of the system (the frequency at which 


the mechanical reactance vanishes) then, w=w,-VK/M. The 


storage and loss components of the effective stiffness are, 


K'= 02M ; K! = ix (17) 


2. Quality Factor and Loss Tangent 
The Q, as seen in Figures 2.1 (a) and 2.2 (a), is an 
indicator of the sharpness of the system response at 
resonance. The Loss Tangent is the ratio of the energy 
dissipated (E" or G") to the average rate of energy stored 


(E’ or G’) per cycle of vibration, or the reciprocal of Q, 
tan = & = EU. a : (18) 


The quality factor for a system can be defined ina 


variety of ways, as follows: 


1. The Q is the ratio of the reactance of the 


system to the resistance, 
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2nf£ 3m 
R 


m 





oi (19) 

2. As seen in Figure 2.1 (a), Q is the ratio of 
the displacement amplitude at resonance to that at zero 
frequency, 


_ X (f=£,) 
~ ~X(#=0) 


(20) 
Equation (20) is true for simple harmonic oscillators 
independent of the value for Q, including Q’s less than 1, 
assuming the stiffness of the system is frequency 
independent. For the case of a viscoelastic material, a 
multiplicative factor, E(w=0)/E(w=w,), is required to make 
Equation 18 valid for dynamic mechanical systems. Here, E 
is the modulus of elasticity corresponding with the mode of 
excitation. 

3. The Q is also proportional to the rate of 


change of the phase at resonance, and is determined by 


taking the derivative of ¢ (Eq. 11) with respect to 


frequency: 
£, ab 
= 2s 21 
0 2 aE ltte sae) 
The derivation [Ref. 11] of this definition is found in 
Appendix A. 
4. In terms of the relaxation time (tT) of the 


system the Q is, 
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OQ=nf.t , (22) 


where t=2m/R, or the time that it takes for the free decay 
amplitude to decrease to 1/e of its initial value, where 


ln e=1 [Ref. 10]. 





Sis Q can also be approximated by, 
Le 
Q= BFE’ (23) 


where f, is the resonant frequency and Af is the full 
bandwidth over which the displacement has dropped to V2 of 
its value at resonance (valid for Q>5). 

6. The Q can also be found from a pole and zero 
plot of the mechanical admittance of the oscillator. If the 
characteristic roots (poles) of the complex admittance of 
the oscillator are expressed as y = a+jb and 7" = a-jb, 


from, 


ne yb? +a* _ b . (24) 


-2a 2a 
The details of this definition of Q are in Appendix B. 
hes The Q shown in terms of the critical damping 


is 


oO (25) 


= 
26. 


where ¢ is the fraction of critical damping for the system. 


3. Influence of Temperature and Frequency on the 
Properties of Materials 


The mechanical properties of viscoelastic materials 
are greatly influenced by temperature and frequency. In 
general, as the temperature of a material increases the 
storage modulus decreases. 

The transition of a material between ideal elastic 
solid behavior to that of a viscous liquid is typically 
designated by the glass transition temperature (T,). At 
temperatures well above T,, the molecular chains ina 
viscoelastic material are able to slip past one another and 
begin to flow due to their increase in energy. The 
increased temperature causes the material to behave more 
like a viscous liquid provided the molecular weight of the 
material is low. For materials which have a high molecular 
weight or molecules which are slightly crosslinked, the 
material will behave in a rubber-like fashion. In general, 
the higher the molecular weight for a material the higher 
the degree of intermolecular bonding. Therefore, more 
energy will be needed to break those bonds resulting ina 
higher melting temperature and/or less ability for the 
molecular chains to flow past one another. [Ref. 8,9,12,13] 

The rubber-like behavior occurs because the 
molecular chains within the material are extremely twisted 
in their natural (unstressed) state, and become untwisted 


whenever a stress is applied to the material. The chains 
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revert back to their natural (twisted) state when the stress 
is removed. This process will continue whenever the stress 
is applied or removed. Therefore, the material is capable 
of highly elastic deformation when in this rubber-like 
state. [Ref. 8,9,12,13] Below T, the elastomer is said to 
be in a glassy state. This glassy state of the material 
originates in the atomic structure of the material. At low 
temperatures, the thermal and free volume energy of the 
material is low and the molecular chains can no longer slip 
past one another. In other words, the molecular chain 
configurations within the material are frozen into place 
causing the material to behave more like an ideal solid. 
(Ref. 8,9,12,13] 

The temperature and frequency dependence of the 
material's properties may be significant in the rubber and 
glass regions. However, this dependence is most pronounced 
when the viscoelastic material is in transition between the 
glass-like (elastic solid) and rubber-like (viscous liquid) 
states. [Ref. 9,12,13] 

A viscoelastic material typically goes through four 
different regions or transitions of dynamic mechanical 
behavior as a function of temperature and frequency. These 
regions are defined as [Ref. 9,12]: 

A glassy region 
A glass-rubber transition region 


A rubbery region 
A flow region 


1.) 
2.) 
a4 
4.) 





These four regions are illustrated in Figure 2.6 below. 
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Figure 2.6 Temperature and frequency dependence of the 
storage modulus and loss tangent for 
viscoelastic materials. 

Figure 2.6 shows that in the glassy region, the 
storage modulus has a maximum value and is relatively 
constant. In the glass-rubber transition region, the 
storage modulus decreases rapidly and the loss tangent 
reaches a maximum value. In the rubbery region the loss 
tangent decreases rapidly and then levels off. The storage 
modulus stabilizes and is relatively independent of 
temperature and frequency. The storage modulus decreases 
rapidly in the flow region; however, only those materials 
which have low molecular weight and virtually no cross 


linking will actually begin to flow [Ref. 9,12]. 
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C. MODES OF A "FREE-FREE" BAR 

A long, thin rod made of an isotropic, homogeneous solid 
with length (L) and constant cross-sectional area will 
propagate vibrational waves in three independent modes, 
provided the wavelength of vibration (A) is much greater 
than the rod diameter. The boundary conditions existing on 
the rod, coupled with the physical dimensions, mass, and the 
elastic moduli will determine the particular resonance 
[Ref. 4]. 

In this application, each end of the rod is "free", 
which corresponds to that of zero moment and zero stress at 
the boundaries [Ref. 10]. The physical dimensions of the 
rod are normally easy to measure. By measuring the 
resonance frequency, the elastic properties of the bar can 
easily be determined as will be explained next. 

1. Longitudinal Mode 


The phase speed for longitudinal waves is [Ref. 8] 


c= |= ; (26) 


where p is the mass density of the material and E is the 
Young's modulus of elasticity. 

The longitudinal resonances of a rod having free- 
free boundary conditions are easily found by noting that the 
speed (c), wavelength (A) and frequency (f) are related by, 


c=\f and that the fundamental wavelength is equal to twice 
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the length of the rod. The higher order modes for a 
nondispersive material are harmonically related and can be 
expressed as integer multiples of the fundamental. Thus, 


the longitudinal resonances of a "free-free" rod are: 


nc, 


Y= = 27 
fn ar 2 LQ Bac. es 3 (27) 





where n is the mode number which corresponds to the number 
of displacement nodes in the bar standing wave. The dynamic 
Young's modulus (E,) can be obtained from the previous two 


equations: 


oF? 
E,=4pL7(—) (28) 


2. Torsional Mode 
The torsional resonances of a rod having free-free 
boundary conditions are similar to that of the longitudinal 
mode, in that both are nondispersive and higher order modes 
will also be integral numbers of the fundamental. The 


torsional resonance frequencies of a "free-free" rod are 


Oty Be tow 6s (29) 
2L 


f7 = 
The phase speed for the torsional waves is 


eee ae al 
” p 


where G is the shear modulus for the material. 
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The dynamic shear modulus can be obtained from the 


previous two equations: 


f- 2 
G = 4pL?(—*) ‘ (31) 


From equations 26 and 29 it can be seen that, for the 
longitudinal and the torsional modes, both the dynamic 
Young's and shear moduli are independent of the cross- 
sectional shape of the rod. 
3. Flexural Mode 
The phase speed for the flexural wave (c,;) is 
dispersive and varies with the square root of the frequency. 


The phase speed for the flexural waves is [Ref. 10] 


Gp = fanixc, . (32) 


where the radius of gyration (x) is given by [Ref. 10] 
= @k 
w= (2) fzds , (33) 


where S is the cross-sectional area of the rod and z is the 
distance of an element above the neutral axis in the 
direction of flexure. For a circular rod of diameter d, x 
is d/4. For a rectangular bar of thickness t, x is t/V12 
[Ref. 10]. For a hollow rod the radius of gyration is 
V(od?-id?)/2, where od is the outer diameter and id is the 


inner diameter of the rod. 
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The application of the free-free boundary condition 
for flexural vibrations leads to a series of resonances 
given by [Ref. 8] 


mn?cx 
—— i 


—_— 
fa = or 


n = 3,.0112,4.9994,7,9,11,... , (34) 


where L is the length, c,; is the flexural wave speed, n is 
the mode number and x is the radius of gyration. This 
result is accurate at low frequencies where the effects of 
rotary inertia and shear deformations associated with the 
flexure can be neglected. [Ref. 4,14,15] 

The dynamic Young's modulus can be obtained from 


Equations 32 and 34: 


1024 pL’, fn. 
ae ge ee * si 

Both the longitudinal and flexural modes are 
proportional to the Young's modulus. Typically, the 
resonant frequencies corresponding to the longitudinal 
(compressional) modes are an order of magnitude higher than 
the resonant frequencies of the flexural mode. Because the 
redundancy in the measurement technique is easily exploited, 
using both modes provides both an internal check of the data 
and/or an extended frequency range over which the Young's 


modulus can be determined. 
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4. Correction for Added End Masses 

The coils and the epoxy add additional mass and 
stiffness to the sample which is under investigation. The 
additional stiffness is usually small and can be neglected 
for those samples which are stiff enough to support their 
own weight and shape. The added mass of the transducer 
coils is generally only a few percent of the total mass of 
the rod, but causes slight changes in the resonance 
frequency. This change can, to first order, be accounted 
for by the substitution of an effective length (L.) term 
into the equations for the dynamic moduli [Ref. 4]. 

The effective length term can be found from 
Rayleigh’s method of equating the maximum kinetic and 
potential energies of a vibrating system to obtain the 
resonance frequencies of the rod [Ref. 16]. The rod is 
considered to be a simple harmonic oscillator with potential 
energy of 1/2Kx,? and with kinetic energy of 1/2mv,?. Noting 
v,=wx, [Ref. 10] and equating the potential and kinetic 


energies yields, 
+Mw?x? = 1Kx? = w= = : (36) 


where K is the spring constant, M is the mass of the system 
and w is the angular frequency. The change in resonance 
frequency in terms of the change in the kinetic energy is 


expressed as [Ref. 4] 
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Sf_d0__1 KE (37) 
2 KE ~ 


Garrett has shown [Ref. 4] that the relative 


shift in the kinetic energy is 


SKE _ 2m | (38) 
KE  M 


Thus, the effective length for the longitudinal mode 


of vibration having mass added to each end (1/2 m) is 
Lete = L(1+2) (39) 


and the Young’s modulus is now expressed as, 


ge 2 
E = 4p (Leze)” (—) . (40) 


The same procedure is used to find the effective 
length in the torsional case. The difference is that the 
vibration in the torsional mode is generated by the 
restoring force within the material and the moment of 
inertia due to the mass of the rod. The kinetic energy 


related to the rotation of the rod is 
KE = 41(w0,)? : (41) 


It has also been shown [Ref. 16] that the relative change in 


kinetic energy for the torsional mode is 
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SKE 2 Am (42) 


As was the case with the longitudinal mode, the length of 
the rod is inversely proportional to the torsional resonant 
frequency. The relative shift in the frequency and 


consequently the change in the relative length of the rod is 
Ee os ON os SORE eM (43) 
and the effective length for the torsional mode becomes, 
Ee, 6 ttiasey , (44) 
M 


The corrected shear modulus thus becomes, 


G = 4p (Lee) (f7/n)* . (45) 
The relative shift in the kinetic energy for the 


transverse (flexural) mode of vibration is [Ref. 4] 


8KE _ 4.1m, 4m 
s+ z (46) 


The fundamental resonant frequency for the flexural mode is 
proportional to the square root of the energy ratio and 
inversely proportional to the square of the length of the 


rod. The relative change in length is 


8L 
L 


am 
2 KE M 





df _ -1, <1 3KE) 
“2a 
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This leads to an effective length correction term for the 


flexural mode which is identical to the longitudinal mode, 
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Les L(l + 2B) (48) 
and the Young's modulus for the flexural mode becomes, 
z= 1028 ee cay. : (49) 
The effective length correction accounts only for 
the added mass presented by the addition of the transducer 
coils. For compliant samples the additional stiffness added 
by the coils would need to be accounted for as well. 
However, to a first approximation, the argument could be 
made that the added stiffness would offset the correction 
for the added mass, making neither correction necessary. 
Details on how the transducer wire coils are attached to the 
rod in order to excite and detect the vibrations in the 
sample will be presented in Chapter III. 
D. MASTER CURVES: PRESENTATION OF DATA AS A FUNCTION OF 
REDUCED FREQUENCY 
1. Introduction 
When a sinusoidal stress is applied to a material, 
the frequency and temperature dependence of both the storage 
moduli and loss tangent are related. This relationship is 
known as the time-temperature superposition principle, and 
is illustrated by the fact that the same effect upon the 
dynamic material properties can be generated either by 


varying the temperature and holding the frequency constant 
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or by varying the frequency in an opposing manner while 
maintaining constant temperature. 

It is difficult to distinguish between the frequency 
and the temperature dependence of the viscoelastic 
properties of the material. As was seen earlier in Figure 
2.6, the dependence of the material properties upon 
temperature and frequency changes dramatically in the glass 
transition zone. This shows that the relationship between 
frequency and temperature is even more dramatic in this 
region. A method of reducing the variables (frequency and 
temperature) was developed to simplify the separation of the 
frequency and temperature dependence of the material 
properties to a dependence upon a single variable. This 
dependency is normally expressed in terms of a shift factor, 
a;, which combines the temperature and frequency into one 
variable. The shift factor is a monotonically decreasing 
function of temperature. [Ref. 8] 

2. Williams-Landel-Ferry Equation 

For many viscoelastic materials, the general form 

for the shift factor (a;) is expressed by the Williams- 


Landel-Ferry (WLF) equation [Ref. 8] as, 
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i) 
c, (T-T, 
10589" Eee oe 
2 ™ £6 


where c,° and c,° are empirically determined constants, T is 
the measured absolute temperature in kelvin (K), T, is a 
chosen reference temperature also given in K. 

The constants c,° and c,° have been determined 
previously for Polymethyl methacrylate (PMMA) by Ferry 
[Ref. 8] and for PR-1592 by Capps [Ref. 9]. The equation 
for the corresponding shift factor, as given by Ferry, for 


PMMA is 


-21.5 (7-211) 
Log a, = —22:5(7-211)) 51 
94" 43.1+T- 211 (51) 


and as given by Capps for PR-1592 is 


-12.9 (T-283.15) 
Log ap =o 3 52 
9 4@r~ 307 + T- 283.15 (82) 


The approximation represented by the WLF equation 
works well for materials whose glass transition temperature 
happens to fall within the frequency range of interest such 
as, PR1592. For other materials such as, PMMA and 
polycarbonate, this is not the case. Without T,, the 
reference temperature becomes somewhat arbitrary and leads 
to confusion when analyzing data. The other problem with 
the WLF equation can be seen in Equations 51 and 52 above. 
That is the reference temperatures used in the two equations 


are different. Even if the materials had been the same, the 
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data presented by each could not be readily compared. The 
data are presented here for all three materials tested in 
both raw form and in terms of the shift factor for the 
broadest possible comparison. 

There are alternative methods for approximating a 
material’s dependence upon temperature and frequency. One 
such method is that of a fractional calculus model as 
described by Bagley in [Ref 17]. However, until an 
approximation is demonstrated to be suitable for the 
majority of viscoelastic materials, it may be better to 


simply present the raw data. 
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III. EXPERIMENTAL METHOD 


A. INTRODUCTION 

In the following sections the technique for measuring 
the dynamic properties of materials will be outlined. It 
should be noted that the sample rod is assumed to be a 
lumped acoustical system, and all the equations for the 
resonance and Q described in Chapter II pertain to the 
following discussion. 

This measurement technique is both precise and 
relatively inexpensive to perform, and since it is resonant- 
based, improves the signal-to-noise ratio. The resonant 
frequency, combined with knowledge of the dimensions and 
mass of a homogeneous and isotropic sample, allow 
computation of the material's dynamic moduli. 

The resonant frequencies of the longitudinal, torsional 
and flexural modes of vibration are measured independently 
using a single sample. Thus, both the dynamic shear and 
Young's modulus can be obtained. This is quite convenient, 
since only two independent moduli are needed to completely 
define the entire set of the mechanical properties of an 
isotropic, homogeneous elastomeric material (Ref. 11]. 

The basic experimental instrumentation consists of the 


following components: 
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A test sample made from material of interest. 


A transducer coil mounted on each ends of the 
test sample. 


A suspension apparatus used to support the test 
sample in magnetic fields produced by a 
pair of permanent magnets. 
Electronic instrumentation used to excite and 
detect acoustic vibration of the sample. 
B. TEST SAMPLES 
1. Sample Preparation 


Any nonferrous material which can be cast, drawn, or 


machined into the shape of a long, thin rod with constant 


cross section can be investigated. This is true regardless 


of what that cross section is. 
2. Transducer Coils 

The transducer coils mounted on either end of the 
sample are identical. Reciprocity dictates that either end 
of the rod can be used as the driver or the receiver. The 
coils are comprised of No. 32 gauge insulated copper wire 
which is tightly wound into a 4 turn coil. Each 4 turn coil 
weighs approximately 0.3 g. While the actual coil 
parameters are not of particular importance, the wire used 
should be small to reduce mass loading effects on the end of 
the sample but have sufficient current capacity (typically 
50-200 Ma,,,) [Ref. 4]. 

The coils are shaped to fit the sample and aligned 


in a similar orientation on each end of the sample. The 





coils are then attached using a general purpose five-minute 
epoxy. The ends of the insulated coil are stripped of their 
insulation and solder is applied to ensure good 
conductivity, and a simple continuity check is performed. 
The 4 turn coils typically have a resistance of 1.1 Qf. 

It should be noted that there are compromises 
associated with the size and number of turns used in the 
transducer coils. For example, if the turns ratio is 
increased by a factor of N, the received signal will be 
increased by a factor of N’. This results from a factor of 
N more force delivered to the sample and a factor of N 
higher receive sensitivity. This will be accomplished, 
however, at the expense of added mass to the ends of the 
sample. The larger the mass of the coil, the greater the 
end mass correction. [Ref. 5] 

3. Physical Properties of Test Samples 
The physical properties for each of the finished 


test samples are listed in Table 3.1. 
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TABLE 3.1 PHYSICAL PROPERTIES OF TEST SAMPLE 


L oe Oe | a bala: A M “we () | je Cem) | 
[___rowselsons [eee [25 


= M (rod) + m (coils) 





























C. SUSPENSION APPARATUS 

The test apparatus used in this research is a 
modification of one used previously [Ref. 6,7], and is shown 
in Figure 3.1. It provides for sample excitation, position 
adjustment, and suspension. The excitation system is 
comprised of one set of permanent U-shaped magnets located 
at each end of the apparatus, between which the driving and 
receiving coils are placed. The gap between the magnet 
pairs is adjustable from 0 to approximately 3 cm. The 
magnetic field strength within the gap is approximately 2.4 
+ 0.1 KOe (0.24 + 0.002 Tesla) over a temperature range of - 
17°C to 85°C [Ref. 5]. The magnetic remanence for these 


magnets increases at lower temperatures [Ref. 18]. 








Figure 3.1 Resonance based dynamic moduli measurement 
apparatus. 


D. ENVIRONMENTAL CHAMBER 

An environmental chamber is used to house the samples, 
thereby permitting temperature control. The BHD-408 Bench 
Top Temperature and Humidity Test Chamber (Associated 
Environmental Systems) used [Appendix C] has been modified 
to accept an external CO, canister used to lower the minimum 
temperature to approximately -65°C. The upper limit of this 
system is 145°C. The temperature control stability of this 
system is specified at +0.5°F (=%0.3°C). The internal 
dimensions of the chamber are 24"x 24"x 24". The 
environmental chamber is capable of either manual or 


automated operation using an IEEE 488 interface. 
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E. SELECTIVE EXCITATION AND DETECTION OF DESIRED RESONANT 
MODES 


One of the most important features of the free-free bar 
method of measuring the dynamic modulus is the ability to 
excite the longitudinal, torsional and flexural modes of 
resonance independently using the same inexpensive 
transducer. The sample rod of the material is positioned in 
the measurement apparatus such that the transducers mounted 
on the ends of the rod lie within the gap created between 
the permanent magnet poles. This places the transducers in 
the vicinity of the maximum magnetic field strength during 
measurements. The orientation of the transducer coils 
within the magnetic field determines which resonance mode is 
excited. The induced voltage which results from the 
movement of the coils within the magnetic field can be 


expressed mathematically as follows [Ref. 4] 


= -4 f z. (53) 
ae ge) Bada ’ 


where Da is in an incremental area subtended by the 
transducer wire and B, the magnetic field. The induced emf 
for a small segment of wire, length 1, moving with velocity 
u, in B is 


eee (54) 
emf = B-lxu 
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1. Torsional Mode 

The torsional mode of vibration is excited when the 
transducer coil lies within the magnet pair as shown in 
Figure 3.2, and the coil driven sinusoidally. This creates 
a force which produces a torque on the rod, thus generating 
torsional waves which propagate within the rod. 

A similarly placed transducer at the opposite end of 
the sample responds to the torsional vibration and generates 
an electro-motive force (emf) within the magnetic field 
created by the permanent magnet. The flux within the coil 
will vary as a result of the change in the angle of the coil 


relative to the magnetic field. 





Figure 3.2 Orientation of the transducer coil within 
gap of the magnet for excitation and 
detection of the torsional mode. The arrows 
indicate the direction of motion created by 
the electromagnetic force on the coil. 
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2. Flexural mode 

The flexural mode is excited by rotating the sample 
by 90° and translating the rod up or down approximately a 
distance equal to one half of the bar diameter such that one 
section of the long transducer coil sections lies ina 
stronger magnetic field as shown in Figure 3.3. The forces 
generated by the two transducer coil sections will be in 
opposite directions; however, the positions of the coil 
sections within the magnetic field gradient are such that 
the forces will not be equal. The two opposing forces will 
cause flexural vibrations to propagate along the axis of the 
rod. The flexural motion of the vibration will cause the 
similarly-oriented receiving transducer coil at the opposite 
end of the sample rod to move up and down within the field. 
The change in flux through the transducer coil, induced by 


this movement, generates an emf. 
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Figure 3.3 Orientation of the transducer coil within 
gap of the magnet for excitation and 
detection of the flexural mode. The arrows 
indicate the direction of motion created by 
the electromagnetic force on the coil. 

3. Longitudinal Mode 

The longitudinal mode of vibration is excited by 
increasing the distance between the magnets so that the 
magnetic field is concentrated on the small, vertical, 
section of coil that spans the diameter of the rod at the 
end of the sample as shown in Figure 3.4. The currents 
within this small section of the transducer coil generate a 


longitudinal force on the end which will propagate as 


longitudinal vibration along the length of the rod. 
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Figure 3.4 Orientation of the transducer coil within 
gap of the magnets for excitation and 
detection of the longitudinal mode. The 
arrows indicate the direction of motion 
created by the electromagnetic force on the 
coil. 

The longitudinal motion of the vibration will cause 
the receiving transducer coil at the opposite end of the 
sample rod to move in and out of the magnetic field. The 


change in flux through the transducer coil, induced by this 


movement, generates an emf. 
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F. MATERIAL PROPERTIES AT ROOM TEMPERATURE 
1. Introduction 

The electronic instrumentation used to detect the 
induced emf can range in sophistication from an ordinary 
oscillator and voltmeter to a computer-controlled logging 
and analysis system. There are advantages and disadvantages 
associated with the chosen instrumentation. The basic 
components needed are: 

1. A suspension apparatus 

2. A driver 

3. <A receiver 

4. Magnets 

2. Electronic Instrumentation 
a. Fundamental Instrumentation 
The fundamental instrumentation scheme needed for 

taking measurements consists of an audio oscillator or 
function generator to drive the sample rod and a voltmeter, 
frequency counter, and/or an oscilloscope to measure the 
received signal. The advantages of using these instruments 
are the relatively low cost, availability, and simplicity of 
use. The disadvantages for using only these instruments are 
the limited capabilities, and the time involved in taking 
point-by-point measurements. The schematic representations 
for the various configurations of these components are shown 


in Figure 3.5. 
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Figure 3.5 Block diagram of a simple instrumentation 
setup for room temperature measurements. 


Here, this simple instrumentation scheme is further enhanced 
by adding additional electronic equipment. An Ithaco 1201 
pre-amplifier (with filter) is used to amplify the signal 
from the receiver coil. A high pass-filter with a 10 Hz 
rolloff connected to the preamplifier effectively removes 
the fluctuations on the oscilloscope caused by the soft 
suspension system or by external environmental vibration. 

On the driving side, an HP492A power amplifier improves the 
coupling of the driver's output impedance to the low- 
impedance transducer. 


b. Instrumentation Used 





More sophisticated equipment is available which 
will enhance the capabilities and decrease the time needed 


to make these measurements. An HP3562A spectrum analyzer is 
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used for this purpose. The spectrum analyzer is used as 
both the signal source and the receiver, so that the 
frequency response of the material sample is readily 
displayed on the CRT. The HP3562A also has many other 
useful features, such as pole~-zero curve fit capabilities 
and the ability to measure the frequency response (magnitude 
and the phase change) with respect to a reference input. 
With the HP3562A, one can manipulate the data with simple 
math calculations such as multiplying by i/jw to produce the 
displacement response of the system. The spectrum analyzer 
also makes it very simple to identify the various resonant 
modes. A schematic representation of the instrumentation 
actually used for the room temperature evaluation is shown 
in Figure 3.6. 


HP492A 


SOURCE cue 







| SAMPLE 


| so 


SUSPENSION 
SYSTEM 
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ITHACO 1201 
PREAMPLIFIER, 
FILTER 








Figure 3.6 Block diagram of the FFT based 
instrumentation used for room temperature 
measurements. 
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c. Lock-in Amplifier 

The lock-in amplifier is a phase sensitive 
detector and low pass filter which is capable of detecting 
small signals in the presence of noise. The lock-in 
amplifier is an essential part of the Phase Locked Loop 
(PLL) and is described in detail in section G of this 
chapter. The lock-in amplifier is also an excellent 
instrument for taking room temperature measurements. 

3. Determination of Resonance Frequency 

As shown in Figure 2.5(a) the maximum velocity 
response can always be used to find the resonance frequency. 
Another useful method for determining the resonance 
frequency is to determine where the phase crosses 0° (for 
velocity measurements) using the HP3562A spectrum analyzer 
as described in Figures 2.4(b) and 2.5(b). Probably the 
most accurate method for measuring the resonance frequency 
is to use the pole-zero curve fit from the HP3562A as 
described in Appendix B. 

4. Measurement of the Quality Factor (Q) 

The main method used for determining the quality 
factor (Q) here is that of using the pole zero capabilities 
of the HP3562A dynamic signal analyzer. The HP3562A dynamic 
signal analyzer has an internal curve fitting algorithm 
which provides a synthesized curve fit of poles and zeroes. 


The Q and resonance frequency for the fundamental resonance 








and its overtones are easily obtained from these poles and 
zeroes. This technique has been described by Brown [Ref. 
19] and is presented in further detail in Appendix B. As 
further verification of the Q for the system was also found 
by measuring the resonance frequency and the frequencies 
where the amplitude is V2 (3dB) down from the resonance 


value. This is a valid approximation when the Q is greater 


than 5. 


G. RESONANCE TRACKING 

1. Introduction 

When the temperature of the sample material is 

varied, the resonance frequency and the Q of the material 
also vary. The resonant free-free bar technique described 
earlier in Section D, can be used to measure the storage 
modulus and the loss modulus components as function of 
temperature and frequency [Ref. 4,7]. Since this is a 
resonant based method, the sample must be kept on resonance 
while the temperature is changing. This can be accomplished 
manually by readjusting the drive frequency to coincide with 
resonance every time the temperature changes. Here, a two 
channel phase-locked loop (PLL) is used to continuously 
track the resonant frequency as a function of temperature. 
As will be shown later in this chapter, the quadrature 
signal from the PLL is used to control the frequency of a 


Voltage Controlled Oscillator (VCO) and to maintain the 
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sample at resonance. The in-phase amplitude of the PLL is 
shown in Chapter III, Section G, Subsection 2 to be 
proportional to Q and is used to determine the change in the 
damping properties of the material as the temperature and 
frequency are varied. 
2. Phase Locked Loop 

The phase locked loop (PLL) is used to ensure that 
the sample will remain at resonance while the temperature is 
varied. 

a. Introduction to PLL circuits 

A PLL consists of three main components, a phase 

sensitive detector (a mixer or multiplier), a low pass 
filter, and a voltage controlled oscillator as shown in 


Figure 3.7. 











Voltage 
Controlled 
Oscillator 


Figure 3.7 Block diagram of a basic PLL circuit. 
A mixer is a circuit element whose output is the 
product of its two inputs. Typically, one input is the 


signal of interest and the other is a reference signal. 
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When the two inputs are sinusoidal, the output from the 
mixer is two sinusoids having frequencies equal to the sum 


and difference frequencies of the two inputs. [Ref. 20] 
e(t) = cos [(@,-@,) t+ (,-,) ] -5cos [(w,+0,) t+(o,+,)], (55) 


where e(t)=e,+e,, e,=sin(w,t+d,) and e,=sin(w,t+¢,). When the 
two input signals are identical in frequency, the difference 
frequency will be zero; however, there could still be a 
phase difference between the two signals. 

The output signal from the mixer is delivered to 
a low pass filter (LPF). The low pass filter removes the 
sum component from the mixer output signal leaving only the 
difference component. The output from the LPF is a dc 
voltage whose amplitude is proportional to the difference in 
phase between the two signals. A secondary effect of the 
low pass filter is to allow low frequency components of the 
mixer output signal to pass through while the higher 
frequency components are attenuated. This has the effect of 
reducing the noise associated with the mixer output signal. 
The amount of attenuation depends upon the RC time constant 
and the equivalent noise bandwidth (ENBW) of the low pass 
filter. 

Often, the reference signal used with a lock-in 
amplifier is a square wave. In this case, the output from 


the mixer will be composed of a large number of sum and 
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difference frequencies. In this case, the mixer will again 
produce a phase sensitive dc output when the signal 
frequency is equal to the fundamental or odd-harmonic of the 
reference. 

The voltage controlled oscillator (VCO) has a 
free-running frequency which may be modified by the 
application of an external control voltage [Ref. 20]. The 
change in frequency is proportional to this control voltage. 
In other words, the VCO adjusts its output frequency 
according to this feedback voltage. 

The output signal from the VCO is used as the 
input reference signal to the mixer. The circuit is 
"locked" if the frequency of the reference signal (VCO) 
equals the signal of interest and a fixed phase difference 
(0 degrees in our case) between the two signals is 
maintained. 

b. PLL utilizing a two channel lock-in analyzer 

Figure 3.8 shows a schematic diagram of a basic 
PLL circuit utilizing a two channel lock-in analyzer for the 
phase-sensitive detection. A two channel lock-in amplifier 
contains two individual phase-sensitive detector (PSD) 
circuits. Both PSD circuits use the signal driving the 
excitation coil as the reference and the received signal as 


the input. One of the channels uses a reference circuit 
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phase shifter to change the phase of the signal by 90° in 


order to produce a quadrature output. 
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Figure 3.8 Block diagram of the basic instrumentation 
used to track the resonance of the rod as a 
function of temperature. 

Accordingly, two de output signals are provided, one of 
which is in phase with the signal of interest and the other 
in quadrature. 

By utilizing a feedback loop between the 
quadrature output and the reference circuit phase shifter, 


an operational integrator can be formed. The output from 


this integrator can then be used to correct any phase shift 
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between the reference circuit in the VCO and the input 
signal. Whenever a phase shift is detected between the 
quadrature input signal and the reference signal, the 
integrator feedback loop will produce a dec error signal 
proportional to the phase difference between the two 
signals. The effect of the feedback loop is to continually 
force the output from the quadrature mixer to be zero. 

c. Application of PLL to track resonance of sample 

An HP3314A Function Geviacnbee is used to drive 
the sample rod at a predetermined resonance frequency which 
corresponds to a certain temperature. The HP3314A function 
generator also has the capability of providing the reference 
circuit discussed earlier via an internal VCO. The VCO has 
a control range of +10% to -100% of the free-running 
frequency. However, the linear operating range for the VCO 
function appears to be only +10% to -80% of the preset 
output frequency (used to drive the sample rod). Since the 
resonance frequency of the same material increases as the 
temperature decreases, the system is usually locked at a low 
temperature. 

The velocity of the motion of the receiving 
transducer coil within the magnetic field at the opposite 
end of the sample rod induces an emf voltage within the 
coil. Due to the type of rocking motion generated, the 


output signal will have a frequency and a phase related to 
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temperature and the stiffness of the sample rod. This 
voltage signal is filtered and used as the input to the two 
channel lock-in amplifier. A square wave trigger output 
from the frequency generator is used as the reference signal 
in the lock-in amplifier. The quadrature output signal of 
the lock-in is used as the feedback loop to the VCO via an 
external low pass filter. When the system is locked, the 
temperature of the sample rod can be varied and the 
resulting change in its resonance frequency will produce a 
phase error at the quadrature mixer. The feedback loop will 
force the reference frequency to change to that of the input 
signal. The quadrature signal will be forced back to zero 
and the in-phase signal will be a maximun. 

The new reference frequency, maximized in-phase 
voltage, and the temperature of the sample rod can be 
recorded. Knowing these three parameters and the physical 
dimensions of the sample rod, the dynamic modulus and the 
loss properties of the sample rod material can be determined 
as a function of both temperature and frequency. 

3. Electronic Instrumentation 
The components which were used for the resonant 
tracking measurement include: 


1. Suspension apparatus (w/ sample and attached 
transducers) 


2. Phase-locked loop consisting of 


a. voltage controlled oscillator 


52 


b. integrating circuit 
c. two channel lock-in analyzer 


3. Controllable temperature chamber 
4. Temperature gage 
5. Frequency counter 


6. Voltmeter 


Figure 3.9 is a block diagram of the system of the 
instrumentation used for the resonance tracking measurement. 
The temperature control chamber stabilizes the test sample 
at the selected temperatures so the measurements can be 
taken while the sample rod is in thermal equilibrium. A 
thermistor and a voltmeter are used to record the actual 
temperature on the sample rod. The frequency counter is 
used to record the resonance frequency of the sample rod at 
various temperatures. The voltmeter is used to record the 
in-phase output voltage which is proportional to the Q. A 
power amplifier and a pre-amp were added to improve 
performance. To facilitate the data collection process, the 
entire measurement system was automated using an HP9836 
computer. Each voltmeter and the frequency counter are 


equipped with an IEEE488 bus interface. 


53 








DIGITAL 
MULTIMETER 
HP 3456A 


ee 






















AMPLIFIER 






FREQUENCY 


COUNTER TEEE 
HP 5384A 






UOLTAGE 





OSCILLOSCOPE 
CHANNEL 1 


Figure 3.9 Block diagram of 
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temperature and 


4. Measurement of the Dynamic Moduli using a PLL 


A sample is first driven 


by adjusting the VCO manually with the integrator shorted 


and the error signal feedback pa 
from the lock-in) so that no err 


input to the vcO. The resonance 


is determined by either the maximum velocity amplitude or 


the zero crossing of the phase b 


and received signal. When the drive frequency from the VCO 
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and the resonant frequency of the sample rod are 
synchronized, the quadrature error signal will be a phase 
sensitive dc voltage as described earlier. The output from 
the quadrature signal channel of the lock-in amplifier is 
set to zero by adjusting an internal phase shifter. The 
integrator is then opened and the feedback loop reconnected. 

Once the control loop becomes locked and stable, the 
temperature of the sample rod can be altered. As the 
temperature changes, the natural resonance frequency of the 
rod shifts. The lock-in detects the frequency shift and 
alters the quadrature output voltage. This voltage is then 
used as a feedback signal to "control" the vco output 
frequency so that it continually matches the resonant 
frequency of the sample rod. The resonance frequency can 
either be recorded directly from the VCO or from an external 
frequency counter. (A thermistor connected to a digital 
multimeter is used to record the sample rod temperature.) 
An HP9000 series computer is used to record the resonant 
frequency, the sample rod temperature and the in-phase 
voltage. The control program was originally written by 
Beaton [Ref. 6], modified by Tan [Ref. 7], and modified 
again by this author [Appendix D]. 

5. Measurement of the Q (and Loss Tangent) using a PLL 
The in-phase voltage from the lock-in amplifier is 


proportional to the Q of the system when the system is 











locked onto resonance. The Q can, therefore, be obtained as 
a function of temperature by recording the in-phase voltage 
from the lock-in amplifier. The relationship between the 
in-phase voltage and the loss tangent can be derived — 
the model of a damped simple harmonic oscillator driven at 
resonance. (Shown previously by Brown and Garrett in 
Reference 5.) 

Resonance of a mechanical system is defined as the 
frequency at which the input mechanical reactance goes to 
zero [Ref. 10]. Therefore, the impedance of the mechanical 
system at resonance is the resistance of the system. Noting 


the following definition of the quality factor, [Ref. 12] 





Q= £ ’ (56) 


where w, is the angular resonance frequency, m is the mass, 
R is the resistance of the system (at resonance is equal to 


the applied force over velocity, R=F/u), the Q then equals, 


QO,mu 
so 57 
Q - . (57) 


Rearranging Equation 55 yields the following expression, 





. (58) 


If the reasonable assumption is made that the 


driving force can be held constant (independent of both 
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temperature and frequency) and that the mass of the sample 
does not change, then Equation 58 is also constant. 

The velocity of the system is proportional to the 
emf generated in the coils of wire. The velocity in the 
above equation for Q can, therefore, be replaced by the in- 
phase voltage (V,) from the lock-in amplifier for a 
particular temperature and frequency. This results ina 


proportionality constant (A) which is derived as follows, 





(59) 


Therefore, the loss tangent as a function of temperature and 


frequency is, 


1 al: 
tan6 = = = ———._.. 60 
Q A Vin®o ( } 


A disadvantage of this system, as can be seen in 
Equation 59, is that an independent measurement of the Q for 
the system must be taken in order to determine the 
proportionality constant, A. The accuracy of this 
calibration is directly dependent upon maintaining the exact 
temperature and frequency between the two measurements. 

Even so, this measurement is easily accomplished by using 
the room temperature method discussed in Section F of this 


chapter. 
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Iv. MATERIAL INVESTIGATIONS 


A. INTRODUCTION 

Data was collected for Polyurethane PR-1592, Polymethyl 
methacrylate (PMMA), and Polycarbonate. Of these materials, 
PR-1592 is the only one which has a glass transition 
temperature within the normal operating range of this 
system. Both the PMMA and the polycarbonate are relatively 
stable materials over the temperature range tested. For 
each material tested, the following sections provide 
discussion and selected data for both samples at room 


temperature and resonant tracking. 


B. POLYURETHANE PR-1592 

The PR-1592 was included in this report in order to 
demonstrate the effect of the glass transition temperature 
upon the complex elastic moduli and loss tangent of a 
material. PR-1592 is a commonly used sonar transducer 
encapsulant. An accurate characterization of its properties 
is a matter of importance. This material was discussed in 
further detail by Tan [Ref. 7]. 

1. Material Properties at Room Temperature 

Table 4.1 presents the results of the torsional mode 


(shear modulus) for the PR-1592 test sample at room 
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temperature. Figure 4.1 clearly shows the first 3 


resonances and their corresponding pole-zero data for each. 


TABLE 4.1 PR-1592 SHEAR MODULUS RESULTS (T=25.4°C). 


f/n Q Q G(Leff) 
(Hz) (3dB) (CURVE) (GPa) 





256.7 2.04 2.95 0.043 


94 4.06 
3.06 2.07 














2. Resonant Tracking Evaluation of PR-1592 

The in-phase voltage from the resonance tracking of 
the PMMA sample was calibrated by directly obtaining the Q 
at several different temperatures using the HP3562A Spectrum 
Analyzer. The proportionality constant (A) obtained from 
this calibration was then used to derive the Loss tangent as 
a function of frequency and temperature. The results of 
these measurements for the first two torsional resonances 
are shown in Tables 4.2 and 4.3. Separate plots of this 
data as a master curve for both of the first two torsional 
modes are presented in Figures 4.6 and 4.7 directly 
following the master curves plotted from the actual 
resonance tracking data as further verification of this 


measurement technique. 
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Figure 4.1 PR-1592 frequency response curves and pole 


zero information for the torsional modes at 
T=25.4°C. 
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TABLE 4.2 IN-PHASE VOLTAGE/QUALITY FACTOR CALIBRATION 
FOR THE FIRST TORSIONAL MODE OF PR-1592. 


wan 
Q/ (£V) 
(msec/V) 





-16.75 787.5 4.76 N/A N/A 
-10.16 3.05 -0.2006 | -19.31 





-5.67 579.3 3.40 


-0.2034 -28.85 


i 
ol 
oO 








427.0 2.87 -0.2288 -29.38 

2.40 -0.2859 -25.44 
3.54 -0.5149 ~30.95 
4.82 -0.8403 -30.67 
67.01 177 22 6.08 -1.0734 -31.97 
74.60 175.8 6.23 -1.2534 -28.27 
84.60 7.40 N/A N/A 
AVERAGE (A): -28.58 








Figures 4.2 and 4.3 are two of the actual data 
curves taken at the opposite extremes of the temperature 
ranges. These two figures are shown to demonstrate the 
effect of the glass transition temperature upon the dynamic 
material properties of viscoelastic materials. All of the 
curves with the data used to determine the proportionality 


constant for PR-1592 can be found in Appendix E. 
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TABLE 4.3 IN-PHASE VOLTAGE/QUALITY FACTOR CALIBRATION 
FOR THE SECOND TORSIONAL MODE OF PR-1592. 


















TEMPERATURE|FREQUENCY| QUALITY | IN-PHASE| "A" 
T Q/ (Fv) 
(°C) (msec/V) 
-i6. 75 5.25 N/A N/A 
-5.67 3.49 -0.1874 | -15.31 
4.50 898.7 2.99 =9,1962 | -18.56 
25.40 4.06 -0.3521 | -21.66 
37.02 4.61 -0.4846 | -19.20 
46.42 4.84 -0.6205 | -18.41 
55.85 5.88 -0.8041 | -18.52 
67.01 6.47 -1.0487 | -16.46 


AVERAGE (A): -18.13 
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Figure 4.2 PR-1592 direct measurement of Q, T=-16.75°C. 
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Figure 4.4 Master curve for the first torsional mode of 
PR-1592. 
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Figure 4.5 Master curve for the second torsional mode 
of PR-1592. 
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Figure 4.6 Master curve for the first torsional mode of 
PR-1592 using data obtained from direct 
measurement of Q at various temperature. 
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Figure 4.7 Master curve for the second torsional mode 
of PR-1592 data obtained from direct 
measurement of Q at various temperature. 
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The glass transition of viscoelastic materials 
explained earlier is clearly seen in Figures 4.4 and 4.5 for 
the first and second torsional modes of the PR-1592 sample. 
Figures 4.6 and 4.7, which are plotted from the data taken 
in the independent measurements for Q, provide further 
verification of this phenomena. During the transition 
between glass-like to rubber-like states the Q reaches its 
minimum value (thus, the loss tangent has its maximum). The 
plots also show that the resonance frequency has higher 
values at the lower temperatures and a lower value at the 
higher temperatures. In these four figures the shear 
modulus and loss tangent are plotted verses reduced 
frequency. The reference temperature used to derive the 
reduced frequency is 283.15 K. 

While Figures 4.6 and 4.7 essentially verify the 
measurements taken using the PLL, they are also an excellent 
example of the of the unique differences in the measurement 
techniques used to obtain the data. The data taken using 
the resonance tracking technique involving the PLL is 
continuous; providing a complete analysis of the material 
for the single resonance over a wide range of temperatures. 
The spectrum analyzer data synthesis provides complete 
information about all of the harmonic resonances for a 


particular vibrational mode at a single temperature. 
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The frequency and temperature dependence of the 
shear modulus for PR-1592 is plotted three dimensionally in 


Figure 4.8. 
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Figure 4.8 Shear modulus as a function of frequency and 
temperature for the 1" and 2™ torsional 
modes of PR-1592. 

It should be noted that the dependence of the 
frequency on the temperature is somewhat exaggerated for the 


second mode. A plot of f,/n might be more appropriate; 
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however, it then becomes very difficult to distinguish 


between the first and second modes. 


Tables 4.4 and 4.5 are a summary of data collected 


for the first and second torsional modes of the PR-1592 


sample. 


found in Append 


ix E. 


The actual and complete data for both modes can be 


TABLE 4.4 PR-1592 RESONANT TRACKING DATA FOR THE 
FIRST TORSIONAL MODE. 





TEMPERATURE 
T (°C) 


15.35 


35.19 





f (Hz) 
733.56 


580.63 


437.20 


338.70 


237.21 








FREQUENCY |SHEAR MODULUS 


G (MPa) 







14.04 


| 25.26 278.03 5.68 -0.3870 
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IN-PHASE VOLTAGE 
Vv (DCcV) 





-0.2092 


-0.2261 


-0.2951 


-0.4877 


-1.4435 





TABLE 4.5 PR~1592 RESONANT TRACKING DATA FOR THE 
SECOND TORSIONAL MODE. 


T ¢°%C) f (Hz) G (Mpa) Vv (DCV) 
sae [nse [en2 [eer 
Tere.se | s00.e [mona 
se [aoa [ones 
45.04 121.8 
se [ee [ae 


C. POLYMETHYL METHACRYLATE (PMMA) 






















-0.5984 








-1.0044 







“1.2205 








Polymethyl Methacrylate (PMMA) is a viscoelastic 
material of moderate to low molecular weight. The PMMA 
proved to be relatively stable within the normal operating 
temperatures (-15°C to 85°C) of the unmodified test chamber. 
The original PMMA rod, however, melted and the transducer 
coils drooped out of the magnetic field during a high 
temperature test at about 90°C. This is clearly an 
indication that the PMMA has a relatively rapid glass to 
viscous liquid transition, probably associated with the 


material’s low molecular weight, as explained earlier in 
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Chapter II, Section A. A second sample was then made which 
was used for the remainder of the test. 
1. Material Properties at Room Temperature 

Tables 4.6 and 4.7 show the results for the room 
temperature evaluation of the torsional and flexural modes 
of the PMMA sample. Figures 4.9 and 4.10 are the 
corresponding plots of the frequency response for the 
torsional and flexural modes, respectively. Figure 4.9 
shows the first four harmonically-related resonances and the 
corresponding pole zero data for each resonance of the 
torsional mode. Figure 4.10 shows the first three 
dispersive resonances of the flexural mode for the same 


sample. 


TABLE 4.6 PMMA SHEAR MODULUS RESULTS (T=23.3°C). 


f f/n Q Q E(Leff) 
(Hz) (Hz) (MEAS) (CURVE) (Gpa) 


2027.0 | 2027.0 | 0 17.29 17 «23 1.862 


4093.3] 2046.6 20.71 | 20.09 1.898 
c298.4| 2099.5 | 21.66 | 22.24 | 1.997 
8378.9} 2094.7 22.53 23.24 1.988 
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TABLE 4.7 PMMA YOUNG’S MODULUS RESULTS (T=21.4°C). 
E(Leff 


f f£/n2 Q Q 
(Hz) (Hz) (3ab) (CURVE) (Gpa) 


1321.9] 26.978 12.72 11.01 


The density for the PMMA sample was 1186.43 kg/m’. 









The effective length for the torsional mode was 0.309 m and 
for the flexural mode was 0.305 m. The dynamic Poisson’s 
ratio (vy) was approximated for the PMMA sample to be .330 
(v=E/2G-1) by using the shear modulus for the first 
torsional mode and the Young’s modulus for the third 
flexural mode. This is only an approximation because the 


measurement temperatures differ slightly. 
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Figure 4.9 PMMA frequency response curves and pole zero 
data for the torsional modes T=23.3°C. 
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Figure 4.10 PMMA frequency response curves and pole zero 
data for the flexural modes at T=21.4°C. 
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2. Resonance Tracking Evaluation of PMMA 

The PMMA sample showed no glass transition within 
the temperature range tested as can be seen in Figures 4.11 
and 4.12 for the first torsional and the first flexural 
modes of the PMMA sample. The reference temperature used to 
derive the reduced frequency is 211 K as shown previously in 
Equation 51. The in-phase voltage for the PMMA sample was 
calibrated by using the value for the Q obtained at room 
temperature for the first torsional and flexural modes. The 
proportionality constant (A) obtained from these 


calibrations were A=-0.0743 and A=~-0.1759, respectively. 























1.00E+10 0.1 
wo 
c kb 
¢p) 
S i 
= Oo 
2 = 
8 0.01 = 
= 
= ro) 
o ” 
< e) 
WwW a 
aks 
(7p) 
1.00E+9 0.001 

















1E-13 16-12 16-11 1610 1E9 1£E8 
REDUCED FREQUENCY (Hz) 


Figure 4.11 Master curve for the first torsional mode 
of PMMA. 
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Pigure 4.12 Master curve for the first flexural mode 
of PMMA. 


The frequency and temperature dependence of the 
shear and Young’s moduli are plotted three dimensionally in 


Figures 4.13 and 4.14, respectively. 
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Figure 4.13 Shear modulus as a function of frequency 
and temperature for PMMA. 
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Figure 4.14 Young’s modulus as a function of frequency 
and temperature for PMMA. 
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Tables 4.8 and 4.9 are a summary of the data 
collected for the first torsional and flexural modes of the 
PMMA sample. The actual and complete data collected for the 


PMMA sample is found in Appendix F. 


TABLE 4.8 PMMA RESONANT TRACKING DATA FOR THE FIRST 
TORSIONAL MODE. 





TEMPERATURE oe ne IN-PHASE VOLTAGE 

T (°C) £ (Hz) G (GPa) Vv (DCV) 

- 5.36 2.08 -0.1777 
5.39 2101.45 2.00 -0.1562 
15.30 1.92 -0.1328 
25.10 2019.47 1.84 -0.1120 
44.95 1911.01 1.65 -0.0757 
54.98 1848.78 1.55 -0.0572 
65.31 1804.99 1.47 -0.0318 
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TABLE 4.9 PMMA RESONANT TRACKING DATA FOR THE FIRST 


FLEXURAL MODE. 


TEMPERATURE ey oe ee | IN-PHASE VOLTAGE 
T (°C) £ (Hz) E (GPa) Vv (DCV) 
- 5.00 5.88 -0.4424 
15.40 254.96 5.37 -0.3062 
25.36 248.79 5.12 ~-0.2475 
35.35 241.26 4.81 -0.2104 
$4.75 228.34 4.31 -0.1054 
65.25 226.86 4.25 ~0.0612 
75.10 243.99 4.92 -0.0304 
D. POLYCARBONATE 














Polycarbonate is a material which had not previously 
been tested with free-free bar technique used in this 
thesis. Polycarbonate is a candidate material for use in 
hydrophone applications, thus the complex moduli and loss 
properties of this material are of significance. 

1. Material Properties at Room Temperature 

Tables 4.10 through 4.12 show the room temperature 
material properties for the three vibrational modes of the 
polycarbonate sample. 


Figures 4.15 through 4.17 are the 


corresponding plots of the frequency response and pole-zero 
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information. The results of a hollow polycarbonate sample 


rod tested at room temperature are in Appendix G. 


TABLE 4.10 POLYCARBONATE SHEAR MODULUS RESULTS FOR 
TORSIONAL MODES AT T=23.0°C. 











G (Leff) 
(GPa) 


1188.5 


1186.7 


1196.6 





1202.5 





1202.3 84.16 











72.81 





TABLE 4.11 POLYCARBONATE YOUNG’S MODULUS RESULTS FOR 
THE FLEXURAL MODES AT T=23.0°C. 


£ £/n2 Q Q E(Leff) 
(Hz) (Hz) (3dB) (GPa) 


123.5 13.62 40.84 





344 13.763 114.16 

















TABLE 4.12 POLYCARBONATE YOUNG’S MODULUS RESULTS FOR 
THE LONGITUDINAL MODES AT T=23.5°C. 












E(Leff) 
(GPa) 
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The density (p) of the polycarbonate sample was 
1192.93 kg/m’. The effective length for the torsional mode 


was 0.3684 m and for the flexural mode was 0.3644 m. 
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Figure 4.15 Polycarbonate frequency response curves and 
pole zero information for the torsional 
modes at T=23.0°C. 
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Figure 4.16 Polycarbonate frequency response curves and 
pole zero information for the flexural 
modes at T=23.0°C. 
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Figure 4.17 Polycarbonate frequency response curves and 
pole zero information for the flexural 
modes at T=23.5°C. 
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The dynamic Poisson’s ratio (v) was approximated for 
the polycarbonate sample to be .366 (v=E/2G-1) by using the 
shear modulus for the first torsional mode and the Young’s 
modulus for the third flexural mode. This is only an 
approximation due to the fact that the resonance frequencies 
involved differ slightly. 

2. Resonance Tracking Evaluation of Polycarbonate 

Figures 4.18 and 4.19 are plots of the shear moduli 
of the first torsional and Young’s first flexural modes of 
the polycarbonate sample verses reduced frequency. At the 
time the data was taken, a WLF equation for the 
polycarbonate sample could not be found. The closest match 
available was that which was used for the PMMA sample. 
Accordingly, the same WLf equation was used for the 
polycarbonate sample. The reference temperature used to 
derive the reduced frequency is 211 K. 

The in-phase voltage for the polycarbonate sample was 
calibrated by using the value for the Q obtained at room 
temperature for the first torsional and flexural modes. The 
proportionality constant (A) obtained from the room 
temperature evaluations of the polycarbonate sample were 
A=-0.2811 and A=-0.7532, respectively. The in-phase voltage 
data for both figures 4.18 and 4.19 has only been verified 
at one temperature. While this data is thought to be 


correct there is at this time no verification of the data. 
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Figure 4.18 Master curve for the first torsional mode 
of polycarbonate. 
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Figure 4.19 Master curve for the first flexural mode 
of polycarbonate. 


The frequency and temperature dependence of the 
shear and Young’s moduli are plotted three dimensionally in 


Figures 4.20 and 4.21, respectively. 
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Figure 4.20 Shear modulus as a function of frequency 
and temperature polycarbonate. 
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Figure 4.21 Young’s modulus as a function of frequency 
and temperature polycarbonate. 
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Tables 4.13 through 4.14 are a summary of the data 


collected for the first torsional and flexural modes of the 


Polycarbonate sample. 


The actual data recorded for the 


polycarbonate sample can be found in Appendix G. 


TABLE 4.13 


TEMPERATURE 
Tr (%¢) 
-14.32 


= 5.21. 


74.84 





83.39 


POLYCARBONATE RESONANT TRACKING DATA FOR THE 


FIRST TORSIONAL MODE. 


FREQUENCY | SHEAR MODULUS 








f (Hz) G (GPa) 
1208.87 0.964 
1195.48 0.926 
1186.29 0.911 
1167.37 0.883 
1157.42 0.868 
1138.74 0.840 
1126.85 0.822 
1120.66 0.813 


IN-PHASE VOLTAGE 
Vv (DCcV) 


-0.1338 
-0.1237 
-0.1670 
-0.1793 
-0.2325 
-0.3043 
-0.3185 
-0.3588 
-0.3769 


-0.3784 


-0.3867 








TABLE 4.14 POLYCARBONATE RESONANT TRACKING DATA FOR THE 
FIRST FLEXURAL MODE. 


TEMPERATURE | FREQUENCY | YOUNGS MODULUS|IN-PHASE VOLTAGE 
T (°C) f (Hz) E (GPa) v (DCV) 











14.86 2.50 -0.3807 
25.47 2.46 -0.4606 
35.45 2.44 -0.4874 
45.41 2.39 -0.4830 
55.37 2.38 ~0.4647 
a 
74.96 2.32 ~0.4740 





E. TRANSFER FUNCTION OF A HIGHLY DAMPED RLC CIRCUIT 

The measurement of very compliant and lossy materials 
(very low Q) poses special problems for this measurement 
method. The frequency response for this material at room 
temperature proved to be virtually flat. This is an 
interesting problem from a resonance-based measurement 
technique perspective. What is and how can the Q actually 
be measured if the frequency response at room temperature is 
essentially flat? How can the resonance be tracked if it 
can’t first be detected? These are just two of the 


questions which had to be answered. 
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The investigation of these types of materials began by 
returning to the basics. From the previous discussions in 
Chapter II, there are several ways to approach these 
problems. The first problem to be resolved was whether or 
not the resonance and the quality factor for a material with 
low Q could actually be measured accurately with the 
equipment available. This was accomplished by simulating the 
mechanical system with a simple RLC circuit. The inductance 
(L) and the capacitance (C) of the circuit are analogous to 
the compliance and inertiance of the system. The resistance 
is analogous to the resistance or real part of the 
impedance. The resonant frequency can be obtained from the 


impedance 


z aA 
Z? = R? + (joL-~—)? 
GG j = 


(61) 


The resonant frequency occurs when the impedance reaches its 
lowest value 
o,=— =. (62) 
The Q can be expressed as, 
- £ = 
pu 2" Nae” Ne ini 
R R R - 
The results of this experiment are summarized in Table 


4.15 below. The curves for both the phase and pole zero 


methods are included in Appendix H. 
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The results of this experiment are summarized in Table 
4.15 below. The curves for both the phase and pole zero 


methods are included in Appendix H. 


TABLE 4.15 RLC EXPERIMENT RESULTS. 


(2) (Hz) (Hz) ey 
4654. | 4654.9. | | 4097 | | 4072 | | ase7 | 97 
ee ee ee 
a 
ae 
ae | 
















90 





V. CONCLUSIONS AND RECOMMENDATIONS 

A. CONCLUSIONS 

The main focus of this effort was to define in detail 
the "free-free" bar technique as it pertains to the 
measurement of the dynamic moduli and loss properties. This 
research encompasses the theory, accuracy, limitations, and 
applications of the "free-free" bar measurement technique. 
This research addressed one inherent advantages and 
limitations associated with using the PLL resonance tracking 
method and the transfer function method. In the transfer 
function method, the parameters of a computer generated 
curve fit was used to determine the resonance frequency and 
Q of the sample. It was pointed out in this research that 
because the receive transducer's output is proportional to 
velocity as opposed to displacement, the maximum amplitude 
could be used to determine the resonance frequency of the 
the system. This is especially critical in the application 
of the PLL resonance tracking in that the maximized in-phase 
voltage is related to the maximum velocity of the sample rod 
and is used to determine the loss properties of the material 
as a function of frequency and temperature. An 
investigation into the ability of this technique to measure 


materials with very low Q was also initiated. 





An advantage of the transfer function technique is that 
the data synthesis performed by the curve fit actually takes 
into account the effects of overlapping modes for any given 
resonance. Therefore, the pole zero curve fit appears to be 
the more accurate way of determining the actual Q and the 
resonant frequency of a given material. In addition, the 
spectrum analyzer is capable of providing information on as 
many resonances as can be excited for each mode at the same 
time in contrast to the frequency tracking in the PLL 
method. The disadvantage of using the spectrum analyzer is 
that the data acquisition was time consuming since an 
automated data acquisition system has yet to be developed. 

The two channel PLL provides continuous tracking of the 
resonance frequency of any given mode continuously in time. 
The in-phase output is proportional to the maximum velocity 
response of the system. This ensures that the damping 
properties can also be measured. Therefore, since the 
resonance of the system occurs at the maximum amplitude of 
the velocity response, the PLL provides an accurate method 
for determining the actual temperature and frequency 
dependance of both the dynamic moduli and the loss tangent. 
Another advantage of using the resonant tracking technique 
is the data acquired using this system is easily manipulated 
by outputting the data into a tab deliminated format and 


transferring to a Macintosh computer. 
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One disadvantage with this technique is that at present 
the system is only capable of tracking one resonance at a 
time. The in-phase voltage still has to be calibrated by 
direct measurement of the Q. The set up and operation time 
required to track the resonance over a temperature range of 
-15°C to 85°C for a particular mode in this research was 
approximately six hours. 

The PR1592 sample showed a definite glass to rubber 
transition within the -15°C to 85°C temperature range of 
this Sasa The independent measurement of the © at various 
temperature utilizing the data synthesis technique involving 
the spectrum analyzer provided verification of the PLL data 
and important insights into the differences between the two 
measurement techniques. Neither the PMMA nor the 
polycarbonate samples showed any visible transitions within 
temperature range of the apparatus. Compliant materials 
that show little temperature and frequency dependence are 
candidates for hydrophone applications. 

The author was responsible for modifications to the 
environmental chamber encorporating the use of CO, as a 
coolant to decrease the low temperature capabilities down to 
-85°C. The increased temperature range will allow 
determination of the glass feneradichan temperatures of highly 
damped materials and in tracking the resonance of materials 


through the glass transition region. 
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The results of the RLC circuit simulation of highly 
damped (low Q) systems prove that it is possible to 
- electronically measure both the Q and the resonance of such 
materials with the existing equipment. The apparent 
differences between the theoretical values and the measured 
values is a direct result of the influence of the impedance 
associated with the inductor. As can be readily seen from 
Table 4.15, the theoretical values and the measured values 
match more closely as the resistance is increased and the 
importance of the inductor’s impedance decreases. In spite 
of this flaw, the LCR circuit did prove the basic point that 
the spectrum analyzer is capable of measuring very low 
values for Q at least for a single degree of freedom systen. 

These facts coupled with the new low temperature 
capabilities provide the basics needed to enable the highly 
damped materials to be cooled and characterized below their 
glass transition temperatures in the region where the 
dynamic moduli is at a maximum for that material. This is 
significant in that in this region the resonance frequency 
for that material is also a maximum. Therefore, the 
resonance frequency should be readily apparent and capable 
of being both detected and either locked onto using the 
phase locked loop technique or monitored using the spectrum 


analyzer. 
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B. RECOMMENDATIONS FOR FURTHER STUDY 

Questions surrounding the actual ability to perform 
resonance~based measurements upon low Q materials still 
remain unanswered. Investigations into automating the 
acquisition of data using the spectrum analyzer would be a 
very worthwhile endeavor. Using several PLL circuits in 
parallel to measure more than one mode at a time may also 
worth investigating but may prove to be cost prohibitive. 
Another possibility is using the spectrum analyzer in 
conjunction with the PLL in order to quantify the 
relationship between the Q and the in-phase voltage. 

The mechanical side of the resonant based measurements 
for compliant materials is still very much a challenge. 
Certainly, the question of how to support these soft 
materials is a problem to be investigated as was evidenced 
by the first PMMA rod which melted causing the coils to 
droop out of the magnetic field. Balancing this problem 
while maintaining the freely suspended condition which now 
exists could prove somewhat tricky. There are several 
possible solutions to this problem including additional 
mechanical supports. Freely and vertically suspending the 
rod vertically using a fine wire in order to retain the 
"free-free" boundary conditions is another method to be 
considered. Rigidly mounting the sample on one end and 
mounting both the transmitting and receiving coils to the 


other end of the bar is another viable alternative. ‘With 
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this alternative, however, cross talk created by the induced 
current in the transmitter and the close proximity of the 
receiver coil will have to be addressed. This effect could, 
however, probably be eliminated by using a third subtraction 
coil outside of the magnetic field. Another alternative 
worth investigating is to use a host bar with known 


properties to support the softer sample. 
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APPENDIX A. DERIVATION FOR DETERMINATION OF Q FROM SLOPE 
OF PHASE CHANGE 





The phase relationship between applied force and the 


resulting displacement of a damped harmonic oscillator is 


@ = tant—_xer_y Al 
Q(1-w7.1) 


where w,,=w/w,. Noting that the derivative of arctan is 


Otan“u _ 1 du 











nen A2 
Ox 1+u? Ox 
The derivative of equation Ai is 
Ob ea 1 1 : 71 (20,6) } 
00-61 1+[ ®r61 j 2 Q(1 -w?,1) Q? (1 -761) 2 A3 
Q(1-07.1) 
Gi) =e Q? (1-@7,;)? id Q(1-2,;) +O 61 (2Q@,.1) } A4 
00,61 Q? (1-2,7) +@261 Q7(1 -0%e1) : 
Simplification of this equation yields, 
2 2 
BH _ O(1-WFe1) +2007 01 a 
00,61 Q7(1 -@7¢1) 240% 
When w=w,, (i.e. resonance) equation A5 reduces to 
ob 7 
~ =2 . 
9@ 61 loess a 7 
Thus, 
ob = 22 ,9-1%¢ 
ag let £, ' °° 2aF*e° = 


where ¢ is in radians. This is an exact result. [Ref. 11] 
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APPENDIX B. DERIVATION OF THE POLE ZERO CURVE FIT TECHNIQUE 
FOR DETERMINING THE RESONANCE FREQUENCY AND Q 
OF A SYSTEM 


399 





The resonance frequency and Q of the modes of a 
free-free bar can be found from a pole-zero plot of the 
mechanical admittance of the oscillator as has been 
described by Brown [Ref. 19]. The theory behind this 
concept is derived. based on approximating the sample rod as 
a damped aetiver harmonic oscillator. The complex admittance 
for the oscillator is (Ref. 12] 


Y* (jw) it jJ® 


[n(jo)? + B¢jo) + X) m 


where M is the mass of the oscillator, R is the mechanical 
resistance, and K is the stiffness. By substituting s=jw 
into the equation, the admittance can be manipulated into 


pole.zero format as follows (Ref. 19] 


i) ee eee es eee ee ae 
M(s?+2g+X) M (s-y)(s-y*) sale 
M 


where y and y° are complex conjugates. The characteristic 


roots for the above equation are [Ref. 10, 19] 


= 2R . 2 R \? 
Yua * Sy * i, ols) us 


Therefore, the transfer admittance function has a pole 
at y and yy and a single zero at w=0. The pole-~zero plot 


is illustrated in Figure B.1. The poles are located on a 


circle of radius w.. 
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Figure B.1 Pole zero plot of complex admittance 
function [Ref 19]. 


The quality factor can be found, if the poles are expressed 


as y = at+tjb and y° = a-jb, as follows 


go Ve B4 


-2a 


The resonance frequency for the oscillator is found by 


2 
@, = yD? +a? = b(a+5 3) =b, B5 


where the second approximation is only valid for low or 


moderate damping. [Ref 19] 





APPENDIX C. ENVIRONMENTAL CONTROL CHAMBER SPECIFICATIONS 
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TEMPERATURE/HUMIDITY 
CHAMBERS 


ASSOCIATEO ENVIRONMENTAL SYSTEMS 


These are af slandard uns — readdy aoe — many om sock, Now you Can 
CCX Ue uw to wat your Needs and yo OUdge... Jencn wo and Flow Mogets — lour 
lemoeramure (anges io Choose Irom — choce Of 6 mormng voumnnes — 15 cferent ures on au. 


BENCH TOP MODELS—8HO Series 


The 8H0-403. 409 o 408 flast nwo argrts indicate cuds loat wortmng volume) are Dench 100 
rats that wall allow testing to meet kngwn Mu. Sid. Humaty Soecricauant These unis ave 
dasrcaity Hurmarty chamoert since they Rave 2 lente Lemoerdie range wren comoared 10 
Wne Floor Models in order 10 Grooerty ConUOl bumudity a reingeravon system ss recuWed. 
Sumole and easy (0 operate The dasic chamoers Nave a digital dry Quid and a cirect setung 
igual “AH conuotler (see canon exias lor aulomauie orogramvmng and recording! Simoly 
connect tne power, U4" water suooty ting and 1 dram line and you we reacy to operate. 
Vemoerature Range: ~ ISCO F to + STC l + OF) 

Humidity Range: 20% to 95% = S% RH wun ine range of = TOC (ES*F) 10 + ASC 18S"F) 
and g limuted by a + 4C OS"F) dew pont temoeranwe using tlandard controts, Oouonal 


MHCOOR E33OF ING TeCONdeTS are zvaHadle Lower humadiues avqilable win couomal Chemical 
dies, 


FLOOR MODELS—HO SERIES - 


This senes of Hurrecity Chama slows you lo cysiom deogn 3 cAamoer to Salisty your needs 
ana your Ovoget four wormang woumes io choose irom CR 16 13 and 64 Fui—Nee em 
perature Ranges 


Hoyemidity Range: 20% to 95% RH mites He aoe 
Smid oy a 24° TS 29°F) dew Dom lemoerature 


20°C (63°F) to +8S*C (185° anos 


With 86651 Programmacie comroier Comic aomy od 2350 AM, 





HOW TO OROER SY MODEL NO. 


Examole: MOOEL A 254 
HO gid | a 
Floor model tie 
00 Senes. wm stanard 
(-70°C w s177*C) paca 4 
t 
64 cudtc lee Gece 
fon standard 











Select the temperature range: 

The 400 teres nas a temperate range of — 1ST o + TTC ITF wo + IF) 

The SOO senes nas a temoerature range of — AC ia + 17TC1 — ATF to + ISO) 
The 200 senes has a iamperanie range d 70°C to +177°C (-100°F to -250°F) 


Select the working rotume? 
Inceeaie Cas the last hwo digits Of tne model number - Sives: C8. 16. 3 o 64 FU, 
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FEATURES Intustrarea on page 4) 


FAILSAFE= An aarusiaole Argh temoerature tadsaie 
1S Bamgaed on ail Associated Emmrommencat System 
Chamoers Wis iactory set at the ultimaie mign tem. 
Detairme Ca0aONNY OF INA CA|MOE Loon receror of 
tne chamoer, the customer should set the dry ould 
conuotier (9 ic d0QvE IME Mighest les! Cemoeraiore 
ang Jonst we <ansaie This well qvewent tne test 
cAamoer Hom a emoerature overrun, 


FUSIBLE LINKS fin tne neater cucusth - for over 
semoerarure ororecuon, 


MULTIPANE VIEWING WINDOWS = 19 Ine coor of 
ail chamoers (aoniona: on EHO Senesi 


INTERNAL UGHT = wan external swuen stancard 
on ail Soon models cationg! On dench tog UNIS 


PILOT LIGHTS = to monntor praoes lunction of wan 
oud systems, 


INTERIORS — All stauvess sree! toe 204 helarc 
werded 


EXTERIORS ~ Heavy gauge cold colled steel win 
two coats of texted enamel! hewsn 


INSULATION — Kign censty, low K cox non semmng 
Adergass, 


VAPOR GENERLIORS = Wren humidity i required 
the controller ngnais ine vaow generator lo wduce 
water vaoor ing ine chamae’, Tas meinog a! creat- 
1g Fegan MumGHIES 1s suCeMOr 12 tne atder methods 
OF SOWING water of Gassing warm air aver wacer 


VAPOR GENERATOR FAILSAFE = aU 3 £5. huma- 
iy chamoers are emnooed wrin a laciory set larisaie 
COMIra! (6 Shuk down the vaoor generaion 4 ine event 
of a na water srtuaiton, 


FREEZE PROTECTION = Since all ALS. Chamoers 
have the caoadsliry of going Celow ine freering paint 
Of water, £acn Mumagity Camoer Aas a Ireete orotec- 
lor wrwcn rs Jactory set at LS7°C + 2SFL Ths com 
Vol acivaces a sovenand varve Io dram water Jrom the 
ehamoer and prevent aamage due to Ireszing of 
water. 


FORCED AIR CIRCULATION = High voiume aie cite 
culation is stangard on ail units (9 achieve unilorruty 
in the wort saace The intenow tan is connected ta an 
eatemal moor wun a stamiess steel sae 


WIRING = Coded readily accestvbie. In accordance 
wun NEC 


CONTROLS ~ Ad located in 2 hinged access gane! 
wan Ml relays aed uses readily accessiole. 


ACCESS PORT - aT crameter access cont aod giug 
ta len ude wall on Soor moaels ana nqnt side wall 
On bench loos 1 ort-QUncAed mm ihe metal — addi- 
bona! ports and plugs avarice. 


DIGITAL TEMPEAATURE & HUMIDITY CONTROL. 
LERS 





ANO READILY ACCESSIBLE. 





FEATURES 


PAILSAFE mands 

CONTROL PANEL a] - VIEWING WINDOW 
FREEZE PROTECTION i 

VAPOR GENERATOR Ss name ts. EXTERIOR HAS 
FAILSAFE a TEXTURED 


: ENAMEL FINISH 
: : INTERIOR 1S 
ACCESS OR —>_ _—_—_——_—__— * STAINLESS STEEL. 


ANO PLUG 
WIRING 1S COOED 

















AIR COOLED 
REFRIGERATION 
SYSTEM 
WATER COOLED 
(optional) 
BHD-L00 Senes HO-400, 500 er 200 Series 
43 | S| oa we | Bl: w& 











Dimensions: invem | invem | inser | indem | insem |} iasem | in/em 
Quistae_Heraar 2es8t_i_ 28/71 | 32/81 44/21} 


42/107 | Sost77? | $$7340 | 48/117 Sasya7 | 70/178 





2aye1 | 0/78 | 36/91 | S2i3e { soz | 22/18) | Bov7CI 
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SENCH TOP MOOELS 
3HG-4c. 405, 408 


UNIT SPSTIFICATIONS 
! 10x TOF wOoHLS 
Teme sooge: <1 Cm - Oe 
tom nee ee 
Teno | 


Towa Asege: 2 IFT ie = ITS 
won Ac A 








Fger mooets 
HO-400, 96. 200 





Teme Aenge, «OTe UTE | amm Aamee: = Fe 2S 


AGORA MOONS | 
som, Aiaae: STR 8 Pe t mom Ames SR 




















F 90, i mm con 6 ome comer | 
eee me 











eat away Maleate 
astae | sresag | Sireladt 


Siar ta ala 


Wacamg Yeiome: | 
4220 cm | 


















CSTICNAL EXTRAS 
1. Srogrammacie Conroier S645! (sandard on af floor modes) 
i 36549! 488 wwertace ( recwed /1 oosont 
2 Proguct Prorector lena concrens 
4 Moneweed 2 Sen 12° Crouar Chan Recoraer 
WL a% 1 4° Winco Door lor 310400 moaes 
5 Yarua Wino Moe: 
7 Srtemas gn san Exemat Swucn (stancard on af floor modes) 
Water Storage Gan Externe 
2 Water Recyorg Sysiem jimcuces onerna moraga Gnx) 
"Q) Water Oermeneraazer Syren 
Esta Oermmeranzer Caanage 
SL Futy Acasa Denes 
‘1 Acceona Srewer 











7 ; 
Cemens | Opn | Onsen | mass! PROGAAMUABLE COMTAQLLEA STAMGAKO | 


MaMa) Badan 
Sabie es 





3 <a An Oa s0-0a AMIS 
ae a oe ee | a a = ao 7 a 
rn a) ow ? aa ni? 
TS, 02) 13, 9] 0, | oe, |, yo, (ra, 2a! 
: ‘ redyt ' | en) 
Td: f Leet | ed ee | 
lets wiwlwiwtes {esa} 
Lies ws Laisxiaeatae : 3 2 sit x tural 
= mt Ti at om i : E = mm an 
5 tay a a! a ee ee ee : i __ a; 
: 1 i : 1 : P | : 1 : 
ee ea ta nin karat av acl vig h ae mh mpage qe Gant earl sar are wel 


P28 Ol mm weet omen eet ee eee j 





+ TE RL cewees semeeee waemeeer 


"Q-400, $00, 300 SEAIES . 


Pee east 
Naa 


aeae ! 


a I, 4 | a | “ | 
lige ns 


“Sebald 
























6 Chamoer Cant (lor EHO“) 

17. Aleor Stand (BHO-405 6 408) 

1@ Uxaed CO} Cooeng Cscicety~ mecty ss 

12 Udwd N2 Coceng Cicaonty 

M1. Nose Regucvon Arcage lor hoor mode t 

21. Water Cooed Concenser 

iter Power Source Aecumemenc: 

TA Descran Oenurmoter lor low 44 AH reauwemems (const ‘aciory 
noe = 





APPENDIX D. RESONANT TRACKING COMPUTER PROGRAM 
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10! RT20DLB : Resonance Tracking System 


20! Written by —: Brian Beaton, Beng Hock Tan, David A. Brown, 
30! Date : 18 Mar 90 

40! Last modified : 21 FEB 92 

50! Features added : 1. End mass correction using effective 
60! length. In order to calculate 

70! corrected moduli. D. Bartlett, 

80! 21 Feb 92 

90! 

100! 2. Input and display total mass. 

110! D. Bartlett, 21 Feb 92 

120! : 

130! 3. Change ruuning time for program 
140! to 5 hours with measurements taken 
150! every 2 minutes. D. Bartlett, 

160! 21 Feb 92 

170! 

180! 

190! 


9 00 atasatahaaaniaaaaaneanailaaaioaraaieiaiaaiasaiaiaaiailaseaiaeediaeiaiaseaaaiasesesaaia 


210! Program initialization 
900) PERE EAE CEASA IAR EITC TATA TAT TAA TATA IA AA AAA 


230! 

240 MASS STORAGE IS ":INTERNAL,4,0" 

250 RE-STORE "RT20B:INTERNAL,4,1"” 

260! 

270 OPTION BASE 0 

280 DIM Thermistor(0:240), Resfreq(0:240), Time(0:240), 
Inphmag(0:240) 

290 DIM Mod(0:240), X(0:240), Y(0:240), Alphat(0:240), 
Redfreq(0:240) 

300 DIM Temp(0:240) 

310 DIM Gtor(0:240),Elong(0:240),Eflex(0:240) 

320 DIM Label$(0:25)[50] 

330 DIM Main_title$[50],Sub_title$[50] 

340 DIM X_axis_name$[50], Y_axis_name$[50] 

350 DIM Block1$[50],Block2$(50],Block3$[50],Block4$[50] 

360 DIM C1$[25],C2$[25],C3$[25],C4a$[25],C4b$[25],C5$[25] 

370 DIM C6$[25],C7$[25],C8$[25],C9$[25],C10$[25] 
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DIM L1$[50],L2$[50],L3$[50],L4$[50] 


Clear$=CHR$(255)&CHR$(75) ! CLEAR SCRN key 
Home$=CHR$(255)&CHRS(75) ! HOME key 
Flag=0 'Flag=1 to indicate editing data 
Tmode=1 ! FIRST torsional mode 
Lmode=1 ! FIRST longitudinal mode 
Fmode=3.0112 ! FIRST flexural mode 
Sample_rate=2 ! 1 sample every 2 minutes 
Arg1=60/Sample_rate ! samples per hour 
Pi=PI 1 3.14159265 
Default_grid$="Partial” ! Full, Partial or No 


Label$(1)="Resonance Frequency vs Time"! Graph Main Titles 
Label$(2)="Temperature was varied"! Sub Titles, Axes names 
Label$(3)="Time, hours" 
Label$(4)="Resonance frequency, Hz" 
Label$(5)="Temperature, C” 
Label$(6)="Resonance Frequency vs Temperature" 
Label$(7)="Temperature vs Time” 
Label$(8)="Frequency variation was observed” 
Label$(9)="Young's Modulus vs Temperature" 
Label$(10)="Shear Modulus vs Temperature” 
Label$(11)="Shear modulus, Pa, x104” 
Label$(12)="Young's modulus, Pa, x104" 
Label$(13)="In Phase Magnitude, DCV" 
Label$(14)="In Phase Mag vs Time” 
Label$(15)="In Phase Mag vs Temp” 
Label$(16)="In Phase Mag vs Freq” 


Penl=1 ! White Default colors 
Pen2=2 ! Red 
Pen3=3 ! Yellow 

Pen4=4 ! Green 

Pen5=5 ! Cyan (greenish blue) 
Pen6=5 ' Blue 
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TIDE RHE RAR HAHA HAHAH HAHA AIHA AAA AAAI AHHH IH ASIA HAHA IA AHA HSA A SIR 


780 Options: ! Allow user to process existing data or collect new data 
TOD) a I HA AI HH IH AAA HAIR AA AHA A RAHA AE AAR HAAS IIASA HAIR IAAT A A. 


800! 


810 OUTPUT KBD;Clear$; ! Clear the CRT 

820 OUTPUT KBD;Home$; ! Home display 

830 STATUS 1,9;Screen ! Get screen width 

840 Center=(Screen-28)/2 ! Leading spaces for centering 
850 GCLEAR 

860 PRINTER IS CRT ! Use CRT for displaying menu 
870 PRINT TABXY(1,1) ! Start at top with blank line 

880 PRINT TAB(Center);"Key Purpose” 

890 PRINT TAB(Center);"---------------------------- 4 

900 PRINT TAB(Center);" 0 Process existing data" 

910 PRINT TAB(Center);"4 Collect new data” 

920 PRINT TAB(Center);" 5 Editing existing data” 

930 FOR Keynumber=0 TO 9 ! Off all keys 

940 ON KEY Keynumber LABEL "" GOSUB Comment] 

950 NEXT Keynumber 

960 ON KEY 0 LABEL "Existing Data" GOTO Existing_data ! Key 0 
970 ON KEY 4 LABEL "New Data" GOTO New_data  ! Turn on Key 4 
980 ON KEY 5 LABEL "Edit Data" GOTO Editing_data !Tum on Key 5 
990 BlinkO: WAIT 1 

1000 DISP "Select an option” 

1010 WAIT 1 

1020 DISP 


1030 GOTO Blink0 
1040 RETURN 


LOS 0)! 2 a ee ee hk He IA IA AAA AAA HH IH A AAA AAA HAA HAA HEIR H AHH HATA IAA K 


1060 New_data: ! Main program driver for collecting new data 
AK O71 0) Gaialutaaaataeaaiaeaiaaeaeiaaaaeaaeiaaieaiateaitaieaieaiia aaa eibiediat deat 


1080 GOSUB Clear_keys 

1090 GOSUB Init_param 
1100 GOSUB Init_periph 
1110 GOSUB Input_run_data 
1120 GOSUB Create_new_file 
1130! GOSUB Draw_freq_time 
1140 GOSUB Draw_mag_time 
1150 GOSUB Take_plot_data 
1160 GOSUB Close_file 
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1170 GOSUB Free_periph 
1180 GOSUB Post_process 
1190 STOP 


LQ AAA AAA AAAI AIA A SAA AAA HAA AAA AHH A SAS AAA AAA AE AEH AEH TAA 


1210 Existing_data: ! Main program driver for process existing data 
120) | A AAI AIA HR IIHT ITAA AAAI IAA TAHA AA AAA AA ITAA AIA AAS AIEEE AAA 


1230 GOSUB Clear_keys 
1240 GOSUB Open_old_file 
1250 GOSUB Post_process 
1260 STOP 


JQ TDN RRR A AAA RASA AA AAS AAA AAAS AAA EAH IAI AAI HI AIA A TRIAS ITAA 


1280 Editing_data: | Main program driver to compute Redfreq for old data 
12000 Fe HE RABE ca EEE AROSE EAA ECAC AAI A A TIT TIE AACR A AOTC 


1300 Flag=1 

1310 GOSUB Clear_keys 

1320 GOSUB Open_old_file 
1330 GOSUB Calc_redfreq 
1340 GOSUB Create_new_file 
1350 GOSUB Close_file 

1360 STOP 


13770 FRA REE EEE EEE EEC AG EEE AAT TATE TTT ASAT AA AAA 


1380 Clear_keys: ! "Turns off" all softkeys and clears the screen 
1300 eA AE ICE ECAC AAAI AI AA AR AAA AAA AAA TA AAA AAA TA ATA TAA. 


1400 FOR Keynumber=0 TO 9 
1410 ON KEY Keynumber LABEL "" GOSUB Comment 
1420 NEXT Keynumber 


1430 OUTPUT KBD;Clear$; 
1440 RETURN 
1450 pH ee se doe dab aa EE EASA Ga ESO CE BARGES AATCC TA AAA 


1460 Init_param: ! User sets the minimum and maximum values for initial plot 
ATO RH RR AR I I RIA IC III TCI FR TCR ATOR FIR IR AC IIIA A AIA ATR ARICA HCA ACE 


1480 Center=(Screen-46)/2 ! Leading spaces for centering 

1490 PRINT TABXY(1,1) ! Start at top with blank line 

1500 PRINT TAB(Center);"Key Purpose” 

1510 PRINT TAB(Center);"-------------------------------------- 
1520 PRINT TAB(Center);" 0 Set minimum resonance frequency” 
1530 PRINT TAB(Center);” 1 Set maximum resonance frequency” 
1540 PRINT TAB(Center);" 2 Set minimum in phase magnitude” 
1550 PRINT TAB(Center);” 3 Set maximum in phase magnitude” 
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1560 PRINT TAB(Center);" 5 Set minimum temperature in , C” 
1570 PRINT TAB(Center);" 6 Set maximum temperature in , C" 
1580 PRINT TAB(Center);” 7 Set the default values” 

1590 PRINT TAB(Center);" 8 Review the assigned values” 
1600 PRINT TAB(Center);" 9 To proceed” 

1610! 

1620 ON KEY 0 LABEL "Min Freg, Hz" GOSUB Min_frequency 
1630 ON KEY 1 LABEL "Max Freq, Hz" GOSUB Max_frequency 
1640 ON KEY 2 LABEL "Min Mag, DCV" GOSUB Min_mag 

1650 ON KEY 3 LABEL "Max Mag, DCV" GOSUB Max_mag 

1660 ON KEY 4 LABEL "" GOSUB Comment! 

1670 ON KEY 5 LABEL "Min Temp, C" GOSUB Min_temperature 
1680 ON KEY 6 LABEL "Max Temp, C” GOSUB Max_temperature 
1690 ON KEY 7 LABEL "Default Values” GOSUB Default_values 
1700 ON KEY 8 LABEL "Review values" GOSUB Review_values 
1710 ON KEY 9 LABEL "Proceed" GOTO Moveon 

1720! 

1730 Blinkl: WAIT 1 ! Wait for user selection and 

1740 DISP "Select an option" =! then take appropriate action 
1750 WAIT 1 

1760 DISP 

1770 GOTO Blink1 

1780 Moveon: GOSUB Clear_keys 

1790 DISP "Proceeding .....” 

1800 WAIT 1 

1810 DISP 

1820 RETURN 

1 B30 te Ae eee eke ek te eh i A A AH AAA RH AHA HA AH HAA HR AR HAAR ARAL HAE ATA I AAA. 
1840 Commentl: ! Alerts user when an unassigned soft key is selected 
J BS 01 HH a ee AA I HI A AIH AAA AR AAA AAA HAAR AAA HAA A AAA A AAA HE TA AIA. 
1860 BEEP 300,.1 

1870 DISP "This soft key is unassigned" 

1880 WAIT 1 

1890 DISP 


1900 RETURN 


19] Qt AAA RA IA AA AAA A HAH A ARIA AAA HAA HAI A AA AAA A AHHH AA IAA AAA AK. 


1920 Default_values:! Assigns default values for minimum and maximum T,F 
O30 A A HR AIH HA AH HAAR AHH A ARH AIHA AAAI HATA AAA T EAA ITAA AAA 


1940 User_freq_min=1650 
1950 User_freq_max=1900 


neske 





1960 User_temp_min=0 


1970 User_temp_max=25 

1980 DISP "The Default Values are Set” 
1990 WAIT 1 

2000 GOSUB Review_values 


2010 RETURN 


2020 [2 oC CH A IH ACR I AR I AE CAI HE RR AR HAR EE HE RHE HE IH I HE IR I IIH IE ICR I I 


2030 Min_frequency: ! Accepts user input for minimum frequency 
2DA aR ee edcaidoicciooiide dd odo tag oa aaa TAR TATRA ATTA AA AAA 
2050 INPUT "The min freq for the plot, in Hz 2", User_freq_min 
2060 DISP "The min frequency is set at: ";User_freq_min;" Hz" 
2070 WAIT 1 

2080 DISP 


2090 RETURN 


vl 00! SE HAC ee ee he i HC He CK eH A IH HR HR AH EH AE A AI I I A IC 


2110 Max_frequency: ! Accepts user input for maximum frequency 
5 19) 9 acacia aati siaaaadatiiataatatdiaiiadaadadiaiidinediadiiaiaiaindadaiiitanatl 


2130 INPUT "The max freq for the plot, in Hz ?",User_freq_max 
2140 DISP "The max freq is set at: ";User_freq_max;"Hz” 

2150 WAIT 1 

2160 DISP 


2170 RETURN 


2. 80 Res eee real eoidciar dca aE EAHA EERO ETA CAA EI EEE ATTA TATA AAA 


2190 Min_mag: ! Accepts user input for minimum magnitude 

D2 A AA AA HAH A AI ACR A ATOR AR AR AR AI A RRA AA AH A AAR AAA A TAA ATA AH. 
2210 INPUT "The min mag for the plot, in DCV ?”,User_mag_min 
2220 DISP "The minimum magnitude is set at: ";User_mag_min 
2230 WAIT 1 

2240 DISP 


2250 RETURN 


2260 I HH AHA A HAAR AHA A AMAR AHA HAA K ERLE KARA AAI IK TERIA. HIRE EAE 


2270 Max_mag: ! Accepts user input for maximum magnitude 

QI ROI BAHAR AA AA AAA AHHH HR AAA AH FTIR AIH AIH HAIR ATR AA HH AIH HAHAH HAHA IH ALARA K. 
2290 INPUT "The max mag for the plot, in DCV ?",User_mag_max 
2300 DISP "The maximum magnitude is set at: ";User_mag_max 
2310 WAIT 1 

2320 DISP 


2330 RETURN 


2340! HORAK HK AK HAR HR AH HK HHH AH IKK HI KR KICK HK KKK KKK RRR KKK HAIR KK IK ICE KE 
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2350 Review_values: ! Presents currently assigned values for review 
ATs 0) kaetaaiaaiaieaiaaeieetaa eae aaa 


2370 DISP "Minimum frequency: ";User_freq_min;" Hz" 
2380 WAIT 1 

2390 DISP " Maximum frequency: ";User_freq_max;" Hz" 
2400 WAIT 1 

2410 DISP "Minimum magnitude: ";User_mag_min;" DCV" 
2420 WAIT 1 

2430 DISP " Maximum magnitude: ";User_mag_max;" DCV" 
2440 WAIT 1 

2450 DISP "Minimum temperature: ";User_temp_min;" C" 
2460 WAIT 1 

2470 DISP " Maximum temperature: ";User_temp_max;” C" 
2480 WAIT 1 

2490 DISP 


2500 RETURN 


Pls$I 0 estseahatshaasababaeletiaaeieiadaiaiilaiaiaediaieieiaiaiaiaiaiaiadaiaieieeeiaiainiaiadeaininiedbainseaeaiaseaiaes 


2520 Min_temperature: ! Accepts user input of minimum temperature 
75.3.0) RAAF A a a I a CITE IR AA A AIR A A A AA A A I 
2540 INPUT "Minimum temperature, degrees C ?",User_temp_min 
2550 DISP "The min temp is set at: ";User_temp_min;"Deg C" 

2560 WAIT 1 

2570 DISP 


2580 RETURN 


D5 OO Re HHA B SBSH Ea ABBE EEE CAAA AAA BSAA AAO TAT TA AIR 


2600 Max_temperature: ! Accepts user input of maximum temperature 
261.0) Ea A I I IR IA A TE I I AR I I I RR RR FR TR IR IR FAIR AA IH I I 
2620 INPUT "Maximum temperature, degrees C ?" User_temp_max 
2630 DISP "The max temp is set at: ";User_temp_max;"Deg C" 

2640 WAIT 1 

2650 DISP 


2660 RETURN 


22.670 2A AREER EGE E ACEO ITER ACAI TAT a TIT IR A TTT TTI TA A 


2680 Init_periph: ! Initializes the voltmeter and the frequency counter 
2690)! IH EE ER AIH AH RI AHH IH I KHIR HR AHH HH AH HAH FAIA HHH AA RH HHA HAA HH HK 
2700! 720 - HP5316A Universal Counter 

2710! 722 - HP3456A Digital Voltmeter 

2720! 724 - HP3478A Digital Multimeter (added on 23 Jan 91) 
2730! 

2740 OUTPUT 720;"IN" ! Default state 


a3 





2750 OUTPUT 722;"HF4R1M66STGIOOSTI” ! 2 wire ohms, THMS deg C 
2760 ! 100 line cycles integration 

2770 OUTPUT 724;"FIRANS5T1Z1D1" 

2780 RETURN 

QT HRA AAA AAA AAR ERIKA ASI EAHA FASEB EEE IIR IH HH HEHE EH EATER, 


2800 Free_periph:! Frees the voltmeter and the counter from the HPIB bus 
DB HAHAH AAA HAHAH AI AAA AH ARE RE AHH KIER IAAI REAR EITHER HK 


2820 LOCAL 720 
2830 LOCAL 722 
2840 LOCAL 724 


2850 RETURN 


DRED HEA RHA EAH AAR HARA EAHA SAK KARA AAA HE IK HAHA EHR EAA HE FAH HAHAH TR, 


2870 Input_run_data: ! Accepts user input of selected run data 

23.80 pea Red r Aabaonddibidaiodieiinidaiiccoincicadiiacinaderaicc isaciciisiceice ick ie 
2890 Center=(Screen-42)/2 ! Leading spaces for centering 
2900 PRINT TABXY(i,1) ! Start at top with blank line 

2910 PRINT TAB(Center);"Key Purpose” 

2920 PRINT TAB(Center);"------------------------------------+-+ P 
2930 PRINT TAB(Center);" 0 Enter rod identification” 

2940 PRINT TAB(Center);” 1 Enter the run number for this mode” 
2950 PRINT TAB(Center);" 2 Enter the mode for this run” 
2960 PRINT TAB(Center);" 3 Enter the date for this run” 

2961 PRINT TAB(Center);" 4 Enter the total mass of rod & coils” 
2970 PRINT TAB(Center);" 5 Enter the mass of the rod only" 
2980 PRINT TAB(Center);" 6 Enter the length for this rod” 
2990 PRINT TAB(Center);” 7 Enter the diameter for this rod” 
3000 PRINT TAB(Center);" 8 Review the information entered” 
3010 PRINT TAB(Center);" 9 To proceed” 

3020! 

3030 ON KEY 0 LABEL "Filename (Rod ID)" GOSUB Rod_id 

3040 ON KEY 1 LABEL "Run No.” GOSUB. Run_number 

3050 ON KEY 2 LABEL "Mode" GOSUB Mode 

3060 ON KEY 3 LABEL “Date” GOSUB Date 

3070 ON KEY 4 LABEL "Total mass" GOSUB Total_mass 

3080 ON KEY 5 LABEL "Rod Mass” GOSUB Rod_mass 

3090 ON KEY 6 LABEL "Rod Length" GOSUB Rod_length 

3100 ON KEY 7 LABEL "Rod Diameter" GOSUB Rod_diameter 
3110 ON KEY 8 LABEL "Review entries” GOSUB Review_entrys 
3120 ON KEY 9 LABEL "Proceed" GOTO Onward 
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3130 Blink2: WAIT 1 


3140 DISP "Select an option or proceed” 
3150 WAIT 1 
3160 DISP 


3170 GOTO Blink2 
3180 Onward: Volume=(Length/100.0)*Pi*.25*(Diameter/100.0)42 


3190 Density=(Mass/1000)/Volume 
3191 M=Totalmass-Mass 

3200 GOSUB Clear_keys 

3210 DISP "Proceeding ..... " 

3220 WAIT 1 

3230 DISP 


3240 RETURN 


3250! HAH AHA HHA HE HH AC A HEA AC CH AE EA I a A A IC a a ae a CK i a eo kei ak lei kk 


3260 Rod_id: ! Accepts user input of filename & rod identification 
332770) Pe RRR ORR AEA ECE IAEA AE EAGER ATCC TAA AA AAA AAA 


3280 INPUT "Filename (Rod ID) block label (i.e. ECP4) ?" Block1$ 

3290 WHILE LEN(Block1$)>12 

3300 DISP "Limit Rod Identification to 12 characters" 

3310 BEEP 300,.1 

3320 WATT 1 

3330 DISP 

3340 INPUT "Rod identification block label (i.e. ECP4) 
2" Block1$ 

3350 END WHILE 

3360 DISP "The Rod ID is set at: ";Block1$ 

3370 WAIT 1 

3380 DISP 


3390 RETURN 


IA OO PH eSB EAE AEG A RATCHET TIAA AAA TO ATTA TAR AK 


3410 Run_number: ! Accepts user input of run number 

30.20 pepe sa edadadccedcicaiccioc adda de do ea ag aaa RAGA CE AE ACACIA 
3430 INPUT "Run number for this mode (i.e. 4) ?",Block2$ 

3440 DISP "The Run number is set at: ";Block2$ 

3450 WAIT 1 

3460 DISP 


3470 RETURN 


30.9 0 Je RHE HSER SSSA ESSEC AAA ESE E EGE ASAT AA TAA 


3490 Mode: ! Accepts user input of the mode 


11.5) 





BS OD RH IIH A AR HH AA IH AR HHI AAA A AAA TAA TAHA AEA HASH TIAA HARA. 


3510 INPUT "Mode (Torsional,Longitudinal,Flexural) ?” Block3$ 
3520 SELECT UPC$(Block3$) 

3530 CASE "FLEXURAL","TORSIONAL","LONGITUDINAL" 
3540 DISP "The mode is set as: ":Block3$ 

3550 WAIT 1 

3560 DISP 

3570 CASE ELSE 

3580 BEEP 300,.1 

3590 DISP "Choices are: Torsional,Longitudinal, Flexural” 
3600 WAIT 1 

3610 DISP 

3620 GOTO Mode 


3630 END SELECT 
3640 RETURN 
F650 Pes eae dodrddddcddecdcicaniccicecd dace aaa gaa ASS IAAT AAA 


3660 Date: ! Accepts user input of run date 
GTQ RH HH AR AAAI HH AAA AAAI AA AAA SAA AAAS AA IA HAHA AEA HAHA ISIAH ASA AA SAHA 


3680 INPUT "The date for this sample (i.e. 10 SEP 89) ?",Block4$ 
3690 WHILE LEN(Block4$)>9 

3700 DISP “Limit date entry to 9 characters” 

3710 BEEP 300,.1 

3720 WATT 1 

3730 DISP 

3740 INPUT “Enter the date (i.e.10 SEP 89) ?",Block4$ 
3750 END WHILE 

3760 DISP "The Date is set at: ";Block4$ 

3770 WAIT 1 

3780 DISP 


3790 RETURN 


BTQ |b A HH AIK AA AH AAA AAA AA IH AAA HAHA AAR AAA AAA AA IA AA AAA AA SA AAA A TH AA AE 


3792 Total_mass: ! Accepts user input of total sample mass 

793 1 HH AA AA AH IH AH AA AA AAA AA AA AA SAA AIA H A IAAT AA AHS A SHA AK AAA AAA HER. 
3794 INPUT "The total mass for sample rod (grams) ?”,Totalmass 
3795 DISP “The total mass is set at: ";Totalmass;" grams” 

3796 WAIT 1 

3797 DISP 


3798 RETURN 
39.0010) 32 3 2 3 2 I I I IGA IATA ICICI OCCA IACI CAA ICA A 
3810 Rod_mass: ! Accepts user input of the rod mass 
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3820! HAH AAA HA AHH HH HAH AA HAH HA AA HA AH HAH HHA HAHA HHA HHH HA AA HK AA AAA HAHAHA HACK 


3830 INPUT "The mass for this rod (units: grams) ?",Mass 
3840 DISP "The mass is set at: ";Mass;" grams" 

3850 WAIT 1 

3860 DISP 


3870 RETURN 


RRO) 3 ee eee i ek HH HHH A IIE A HAHAH AHR AAR AAR AERA AEA KEI AAA EAHA AAI 


3890 Rod_length: ! Accepts user input of the rod length 


BOD RAK AAA AAR ARH ARH FRI H AHIR HAE AHIR AREA A AEA A AA TAREE TAKE HATA AEA 


3910 INPUT "The length for this rod (units: cm) ?",Length 
3920 DISP "The length is set at: ";Length;"centimeters” 
3930 WAIT 1 

3940 DISP 


3950 RETURN 


3960! 2 AR IAI IOC FAIA IG I a 3 a ik a oe keke ak 0 2 dk 3 2 i i oe 


3970 Rod_diameter: ! Accepts user input of the rod diameter 

BQ RY HAH AAA AAA AAA HAIR ARR A AH RA AAA AAA HAAR EAHA A AH HAA IAAT AA AIK 
3990 INPUT "The diameter for this rod (units: cm) ?",Diameter 
4000 DISP "The diameter is set at: ";Diameter;" centimeters” 

4010 WAIT 1 

4020 DISP 


4030 RETURN 


0A0 3B ERE EEA OEIC ESSA G AA TOCA A TA TATA IAAT AT ATTA AA AIA 


4050 Review_entrys: ! Presents currently assigned values for review 
AQ EL A AR AAA IA AAA AAAI HA AAA ASIA AHH AEE HERI RAIA AEH TAA 
4070 DISP "Rod ID: ";Block1$ 

4080 WAIT 1 

4090 DISP " Run: ";Block2$ 

4100 WAIT 1 

4110 DISP "Mode: ";Block3$ 

4120 WAIT 1 

4130 DISP " Date: ";Block4$ 

4140 WAIT 1 

4141 DISP "Total Mass";Totalmass;”"grams” 

4142 WAIT 1 

4150 DISP “Rod Mass: ";Mass;"grams” 

4160 WAIT 1 

4170 DISP " Length: ";Length;” centimeters” 

4180 WAIT 1 
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4190 DISP “Diameter: ";Diameter;" centimeters” 
4200 WAIT 1 

4210 DISP 

4220 RETURN 


0730 FOR RARER EARS IAEA ISAT T TATA T ATA AAT TATA TISAI ATA 


4240 Draw_mag_time: ! Produces Magnitude vs Time graph w/o curve 
QDS (a a A a A AA IAHR II A AA AAA AA HAHA AAA AAAI AA AAA AA HAHA AAA TA A AHH RAT 
4260 Main_title$=Label$(14) 

4270 Sub_title$=Label$(2) 

4280 X_axis_name$=Label$(3) 

4290 Y_axis_name$=Label$(13) 

4300 Xmin=0 

4310 Xmax=5 

4320 Ymin=Actual_mag_min 

4330 Ymax=Actual_mag_max 

4340 GOSUB Generic_plot 


4350 RETURN 


13.60 FARRAR AEE dE AE GEE EAE E CISA a ATA IT TAIT TAT AT AT AT TATA 


4370 Draw_freq_time: ! Produces Frequency vs Time graph w/o curve 
03 80 asada dace aad a gaggia a AEA AERA R AECC AGAR ICCC AGAR 

4390 Main_title$=Label$(1) 

4400 Sub_title$=Label$(2) 

4410 X_axis_name$=Label$(3) 

4420 Y_axis_name$=Label$(4) 

4430 Xmin=0 

4440 Xmax=5 

4450 Ymin=Actual_freq_min ! Actual minimum frequency 
4460 Ymax=Actual_freq_max ! Actual maximum frequency 
4470 GOSUB Generic_plot 


4480 RETURN 


A. Q() | 962 9 a A CR AIA A RI AIH AAAI IH A ITH IH RA AHA AIHA AA AHH HAA AAA AT AH HAA HA 


4500 Young_mod: ! Compute Young's modulus in flexure modes 
45 10 2a BR Ee Re EEA ECE SCAG EAE SEI ATTA AAS TTI TA IAAI A 


4520 D=Diameter/100 

4530 L=Length/100 

4531 Leff=L*(1+(M/Mass)) 

4540 Arg6=Density*((32*Leff*2)/(Pi* D*Fmode’2))42 
4550 Eflex(I)=Arg6*Resfreq(I)*2 


4560 RETURN 


IST RRA RA AAA AR AAA AAA AAA A AAA ATA TAR AAA TATA TRH TA TAR AAR TR ATR TOT HA AR IH. 
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4580 Shear_mod: ! Compute Shear Modulus 


LES OA 1 ee ee tte ek hk kk Hk A AA HR AA HR A AHHH AAA HA AAA AHH ARH AR AA AAA H. 


4600 D=Diameter/100 

4610 L=Length/100 

4611 Lefs=L*(1+(2*M/Mass)) 

4620 Arg5=(Density*4*Lefs*2)/(Tmode42) 
4630 Gtor(J)=(Arg5*Resfreq(J)*2) 


4640 RETURN 


PU) azat tazatatasaae aca eleabaitiediuiiababaonaaabainiiaiaiataiinmaiaaaaaaaaia 


4660 Lyoung_mod: | Compute Young's modulus in longitudinal mode 
GTQ) RAE HA RAH HAAR AAAI HHH IAI HAHEI AHHH HH HAHA SAE AHHH AAS AEA AAS H. 


4680 D=Diameter/100 


4690 L=Length/100 

4691 Lefl=L*(1+(M/Mass)) 

4700 Arg7=(Density*4*Lefl42)/(Lmode2) 
4710 Elong(K)=Arg7*Resfreq(K)*2 


4720 RETURN 


LT BQ RHA AA AA AAA AAA AIHA SA IRA H HARIKA HAHEI HA HAI H AAA AIS ATA AH ATA TH 


4740 Compute_redfreq: ! Compute redfreq of new data 
11750 Feb RR EERE AAAS E ECC E IAC AAT ATTA A AAA TAA AT AAA AA AAA 


4760 Temp(I)=Thermistor(I)+273.15 

4770! X(D=-12.9*(Temp(1)-283.15)/(107+Temp(I)-283.15) ! PR 1592 
4780 X(1=-21.5*(Temp(I)-211)/(43.1+Temp(I)-211) ! Plexi-glass o 
4790 Y(I)=X(I)*LOG(10) 

4800 Alphat(I)=EXP(Y(1)) 

4810 Redfreq(I)=Resfreq(I)*Alphat(I) 


4820 RETURN 


4.330) 3A AERA A IAAT AAT ATA TATA TATA TTA TA TAT ITALIA IAS IATA 


4840 Calc_redfreq: ! Convert RT20B to RT20C format 


48.50 ;eR HRA Ere Aaa AEE AEE ACCA AKA EAT AA TAA AT AAAS A TAIRA 


4860 FOR I=0 TO 5*Argl 

4870 Temp(I)=Thermistor(I)+273.15 

4880 X(I)=-12.9*(Temp(I)-283.15)/(107+Temp(I)-283.15)!PR 1592 
4890 Y(D=X(1D*LOG(10) 

4900 Alphat(I)=EXP(Y(1)) 

4910 Redfreq(I)=Resfreq(I)*Alphat(I) 

4920 NEXT I 

4930 RETURN 


QA) a A eR HH AHH AIH AIH HAH HH AIA AHHH HAAR HRA HRI AA IER TAHT TARAS A HRA. 


pi he) 





4950 Compute_modulus: ! Computes appropriate modulus based on mode 
40 60 Bebe paseo a EAE ACRE ASAE CSCC A OAR AAA IA ARIA ATCT TC TATA TAA AK 





4970 DISP "Computing appropriate modulus ........ is 

4980 D=Diameter/100 ! convert to meters 
4990 L=Length/100 ! convert to meters 
4991 Leff=L*(1+(M/Mass)) 

4992 Lefs=L*(1+(2*M/Mass)) 

4993 Lefl=L*(1+(M/Mass)) 

5000 SELECT UPC$(Block3$) 

5010 CASE "FLEXURAL" 

5020 Arg2=Density*((32*Leff*2)/(Pi*D*Fmode’2))2 
5030 Eflex_max=Arg2*Actual_freq_max/2 

5040 Power=LOG(Eflex_max) DIV LOG(10) 

5050 Scale_factor=104Power 

5060 Scale_eflex_max=Eflex_max/Scale_factor 

5070 Scale_eflex_min=(Arg2*Actual_freq_min‘2)/Scale_factor 
5080 FOR I=0 TO 5*Argl 

5090 Eflex(I)=(Arg2*Resfreq(1)42)/Scale_factor 
5100 NEXT I 

5110 CASE "TORSIONAL" 

5120 Arg3=(Density*4*Lefs42)/(Tmode2) 

5130 Gtor_max=Arg3*Actual_freq_max42 

5140 Power=LOG(Gtor_max) DIV LOG(10) 

5150 Scale_factor=104Power 

5160 Scale_gtor_max=Gtor_max/Scale_factor 

5170 Scale_gtor_min=(Arg3*Actual_freq_min‘2)/Scale_factor 
5180 FOR J=0 TO 5*Argl 

5190 Gtor(J)=(Arg3*Resfreq(J)42)/Scale_factor 
5200 NEXT J 

5210 CASE "LONGITUDINAL" 

5220 Arg4=(Density*4*Lefl42)/(Lmode‘2) 

5230 Elong_max=Arg4*Actual_freq_max42 

5240 Power=LOG(Elong_max) DIV LOG(10) 

5250 Scale_factor=10“Power 

5260 Scale_elong_max=Elong_max/Scale_factor 

$270 Scale_elong_min=(Arg4*Actual_freq_min‘2)/Scale_factor 
5280 FOR K=O TO 5*Argl1 

5290 Elong(K)=(Arg4*Resfreq(K)42)/Scale_factor 
5300 NEXT K 


5310 END SELECT 
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5320 
5330 


DISP 
RETURN 


53.40 SEH EEE EERO E EEC EAC T ISA AIA TIA AA TAA ASIA STAI IIIA 


5350 Percent_moduli:! Determines the percent change in modulus per deg C 
S36 HIRI HHH AAA A IAAI HRA A AAA ARAMA AA AAA RAIA HER AAT FAA HER, 


5370 
5380 
5390 
5400 
5410 
5420 
5430 
5440 
5450 
5460 


Delta_temp=Actual_temp_max-Actual_temp_min 
Fsum2pt=Actual_freq_max+Actual_freq_min 
Fdiff2pt=Actual_freq_max-Actual_freq_min 
Fminls=Actual_temp_max*Slope+Intercept 
Fmaxls=Actual_temp_min*Slope+tIntercept 

Fsumls=Fmaxls+Fminls 
Fdiffls=Fmaxls-Fminls 

Prent_per_c_2pt=100*4*(Fdiff2pt/Fsum2pt)/Delta_temp ! 2 pt 

Prent_per_c_ls=100*4*(Fdiffls/Fsumls)/Delta_temp !L squares 
RETURN 


SAT 0 AIHA RIAA AAA AAR HAHA HAA AA AH HA AH HHA HHH KATHI HAMAR HHA IAA 


5480 Mag_temp: ! Produces In phase mag vs Temperature graph w/ curve 
5400 bebe seedaeaadeiedocdecdcecciccae dada aac rca aca AAAI ACCA AAO AAA 


5500 
5510 
5520 
5530 
5540 
5550 
5560 
5570 
5580 
5590 
5600 


Main_title$=Label$(15) 
Sub_title$=Label$(2) 
X_axis_name$=Label$(5) 
Y_axis_name$=Label$(13) 
Xmin=Actual_temp_min 
Xmax=Actual_temp_max 
Ymin=Actual_mag_min 
Ymax=Actual_mag_max 
GOSUB Generic_plot 
CALL Generic_curve(Pen3,Thermistor(*),Inphmag(*),5*Arg1) 
RETURN 


5610 a ea a ee ak HA A HAH AAA AA IH AA AHH AA AH AAR HA ARH AAA AAA AHHH AAA. 


5620 Mag freq: |!Produces In phase mag vs Resonant frequencies 
SE 301 RR AH AIH AAAI AA AHHH RA AAA A HA AIHA A HRA HARSHA HAAR AAAI HAHA 


5640 
5650 
5660 
5670 
5680 
5690 
5700 


Main_title$=Label$(16) 
Sub_title$="_" 
X_axis_name$=Label$(4) 
Y_axis_name$=Label$(13) 
Xmin=Actual_freq_min 
Xmax=Actual_freq_max 
Ymin=Actual_mag_min 





5710 Ymax=Actual_mag_max 

5720 GOSUB Generic_plot 

5730 CALL Generic_curve(Pen3,Resfreq(*),Inphmag(*),5*Arg1) 
5740 RETURN 


EW bat) Estates ateaaeaaalaaaabeataiaiaiaaiaaeateataateaeaeaeaaieaeatleatdaataiateseataiatiaadiaiiaaiaial 


5760 Shear_temp: ! Produces Shear modulus vs Temperature graph w/ curve 
ST TQ ERE A AA AH AAAI A AAA AAA AAA AAAH RASA SAHA AA AAAS AAA HAS AHA A 
5780 Main_title$=Label$(10) 

5790 Sub_ttieS="_" 

5800 X_axis_name$=Label$(5) 

5810 Y_axis_name$=Label$(11)&VAL$(Power) 

5820 Xmin=Actual_temp_min ! Actual minimum temperature 
5830 Xmax=Actual_temp_max ! Actual maximum temperature 
5840 Ymin=Scale_gtor_min ! Scaled minimum shear modulus 
5850 Ymax=Scale_gtor_max ! Scaled maximum shear modulus 
5860 GOSUB Generic_plot 

5870 CALL Generic_curve(Pen3 ,Thermistor(*),Gtor(*),5*Arg1) 


5880 RETURN 


SRO | aa a a a a A HAI A HA HR AA HHA HAH FH A AAA HAHA AAA AAA A HAHA AHA ATA AAI AAA. 


5900 Young_temp:! Produces Young's modulus vs Temperature graph w/ curve 
50 ] Ope aaa aaa ea ara ECCS ICICI IAA CAAA AAA TATA ATH AAAI AA AH A 


5920 Main_title$=Label$(9) 

5930 Sub_title$="_" 

5940 X_axis_name$=Label$(5) 

5950 Y_axis_name$=Label$(12)& VAL$(Power) 

5960 Xmin=Actual_temp_min ! Actual minimum temperature 
5970 Xmax=Actual_temp_max ! Actual maximum temperature 
5980 SELECT UPC$(Block3$) 

5990 CASE "FLEXURAL" 

6000 Ymin=Scale_eflex_min -! Scaled minimum Y (flexural) 

6010 Ymax=Scale_eflex_max ! Scaled maximum Y (flexural) 
6020 GOSUB Generic_plot 

6030 CALL Generic_curve(Pen3,Thermistor(*),Eflex(*),5*Arg1) 
6040 CASE "LONGITUDINAL" 

6050 Ymin=Scale_elong_min ! Scaled min Y (longitudinal) 

6060 Ymax=Scale_elong_max ! Scaled max Y (longitudinal) 

6070 GOSUB Generic_plot 

6080 CALL Generic_curve(Pen3, Thermistor(*),Elong(*),5*Arg1) 


6090 END SELECT 
6100 RETURN 
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OL LQ) a AA AAA AA AA AA AAA AHHH AAA IAA IIASA ASEH AA AAAS HAA TI AAT 


6120 Generic_plot: ! Produces design and layout of line graph 


6] SERRA HAHAHAHA AHH AAR HAA KAA AHHH AA IH EAHA HAHA AHH RIAA HAHA AIA RA 


6140 GINIT ! Initialize various graphics parameters 

6150 GCLEAR ! Clear the graphics display 

6160 GRAPHICS ON ! Turn the graphics display on 

6170 OUTPUT KBD;Clear$; ! Clear the CRT 

6180! 

6190 IF Plot_device$="Plotter" THEN |! Route to HP7475A plotter 

6200 GRAPHICS OFF 

6210 BEEP 1000,.1 

6220 DISP "Press CONTINUE when the plotter is ready" 

6230 PAUSE : 
6240 OUTPUT 705;"IP 1300,1000,9000,6750;"! Set scaling points 

6250 PLOTTER IS 705,"HPGL” 

6260 END IF 

6270! 

6280 Xgdumax=100*MAX(1,RATIO) ! How many gdu’s wide screen is 
6290 Ygdumax=100*MAX(1,1/RATIO) ! How many gdu’s high screen is 


6300 CALL Scale(.02,.98,.10,.98,Top,Left,Xcntr, Ycntr,X gdumax, Ygdumax) 
6310 PEN Pen2 
6320 FRAME 


6330! 

6340 IF LEN(Main_title$)>30 THEN ! Size main title 
6350 CALL Label(4,.6,0,6,Pen2,Xcntr,Top-2,Main_title$) 
6360 ELSE 

6370 CALL Label(6,.6,0,6,Pen2,Xcntr, Top-2,Main_titleS) 
6380 END IF 

6390! 

6400 CALL Label(4,.6,0,6,Pen2,Xcntr,.9*Top,Sub_title$) 


6410 CALL Block_info(Block1$,Block2$,Block3$,Block4$,X gdumax,Pen5,Pen2) 
6420 CALL Scale(.15,.90,.25,.8,Top,Left,Xcntr, Ycntr,Xgdumax, Ygdumax) 
6430 CLIP OFF 


6440 CALL Label(4,.6,0,4,Pen5,Xcntr,-18,X_axis_name$) 
6450 CALL Label(4,.6,90,4,Pen5,-12, Ycntr,Y_axis_name$) 
6460 CLIP ON 

6470! 

6480 SELECT Xscale$ 

6490 CASE “Auto scale X" ! Auto scale X axis 


123) 





6500 CALL Xscale(Xmin,Xmax,Xminor,Xmajor) 


6510 CASE "User scale X" ! User scale X axis 

6520 SELECT Pass$ 

6530 CASE "Initial" ! Initial pass 

6540 Xmin=0 ! Minimum time, Hours 

6550 Xmax=5 ! Maximum tme, Hours 

6560 Xmajor=.5 

6570 Xminor=.1 

6580 CASE "Follow on" ! Use user modified values 
6590 Xmin=Xmin_manual 

6600 Xmax=Xmax_manual 

6610 Xmajor=Xmajor_manual 

6620 Xminor=Xminor_manual 

6630 END SELECT 

6640 END SELECT 

6650! 

6660 SELECT Yscale$ 

6670 CASE "Auto scale Y" ! Auto scale Y axis 

6680 CALL Yscale(Ymin,Ymax, Yminor, Ymajor) 

6690 CASE "User scale Y" ! User scale Y axis 

6700 SELECT Pass$ 

6710 CASE "Initial" ! Initial pass 

6720! Ymin=User_freq_min ! Use user entered min freq 
6730! Ymax=User_freq_max ! Use user entered max freq 
6740 Ymin=User_mag_min ! Use user entered min mag 
6750 Ymax=User_mag_max |! Use user entered max mag 
6760 Ymajor=(Ymax-Ymin)/5 

6770 Yminor=Ymajor/5 

6780 Pass$="Follow on" 

6790 Xscale$="Auto scale X" 

6800 Yscale$="Auto scale Y" 

6810 CASE "Follow on” ! Use user modified values 
6820 Ymin=Ymin_manual 

6830 Ymax=Ymax_manual 

6840 Ymajor=Y major_manual 

6850 Yminor=Ym.> >r_manual 

6860 END SELECT 

6870 END SELECT 

6880 WINDOW Xmin,Xmax,Ymin, Ymax 
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6890! 


6900 CALL Lbl_axes(2,.6,Pen4,Xmin,Xmax,Xmajor, Ymin, Ymax, Ymajor) 
6910 PEN Pen2 

6920 AXES Xminor, Yminor,Xmin, Ymin,Xmajor/Xminor, Ymajor/Yminor,3 
6930 AXES Xminor, Yminor,Xmax, Ymax,Xmajor/Xminor, Ymajor/Yminor,3 
6940! 

6950 SELECT Grid_type$ ! Grid 

6960 CASE "Partial" 

6970 GRID Xmajor, Ymajor,Xmin, Ymin 

6980 CASE "Full” 

6990 GRID Xminor, Yminor,Xmin, Ymin,Xmajor/Xminor, Ymajor/Yminor, | 
7000 END SELECT 

7010! 

7020 PENUP 


7030 GRAPHICS ON 
7040 RETURN 


7050! Se a aac aaa ac aiciack Jaci dak cok ica ccc dtc ck detec 


7060 Take_plot_data: ! Collects freq vs temperature data during the run 
1091 0 faseaataaaoaaaaeeaeaaeaiaiaiaaadaaiaaahailaiaeiaaiaiaaiaiaeseaia 
7080 ON KEY 0 LABEL "Take Data" GOTO Take 

7090 FOR Keynumber=1 TO 9 

7100 ON KEY Keynumber LABEL "" GOSUB Commentl 

7110 NEXT Keynumber 


7120 Wait: GOTO Wait 
7130 Take: ON KEY 0 LABEL "" GOSUB Comment 


7140 PEN Pen3 

7150 OUTPUT KBD;"Taking data ....."; 
7160 TO=TIMEDATE 

7170 SELECT UPCS$(Block3$) 

7180 CASE "FLEXURAL" 

7190 FOR I=0 TO 5*Arg1 

7200 ENTER 720;Resfreq(I) 
7210 ENTER 722;Thermistor(I) 
7220 ENTER 724;Inphmag(1) 
7230 Time(1)=I/Arg1 

7240 PLOT Time(I),Inphmag(I) 
7250 GOSUB Compute_redfreq 
7260 GOSUB Young_mod 

7270 OUTPUT @Path_1;Eflex(I),Thermistor(I), 


Resfreq(I), Inphmag(I),Redfreq(I) 
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7280 
7290 
7300 
7310 
7320 
7330 
7340 
7350 
7360 
7370 
7380 
7390 
7400 


7410 
7420 
7430 
7440 
7450 
7460 
7470 
7480 
7490 
7500 
7510 
7520 
7530 


7540 
7550 
7560 
7570 
7580 
7590 
7600 
7610 
7620 
7630 


DISP I 
WAIT 118.315 !Adjust for 1 sample every 2 min 
NEXT I 
CASE "TORSIONAL" 
FOR J=0 TO 5*Argl 
ENTER 720;Resfreq(J) 
ENTER 722;Themmistor(J) 
ENTER 724;Inphmag(J) 
Time(J)=J/Argl 
PLOT Time(J),Inphmag(J) 
GOSUB Compute_redfreq 
GOSUB-Shear_mod 
OUTPUT @Path_1;Gtor(J), Thermistor(J),Resfreq(J), 
Inphmag(J), Redfreq(J) 
DISP J 
WAIT 118.315 
NEXT J 
CASE "LONGITUDINAL" 
FOR K=0 TO 5*Argl 
ENTER 720;Resfreq(K) 
ENTER 722;Thermistor(K) 
ENTER 724;Inphmag(K) 
Time(K)=K/Arg1 
PLOT Time(K),Inphmag(K) 
GOSUB Compute_redfreq 
GOSUB Lyoung_mod 
OUTPUT @Path_1;Elong(K),Thermistor(K),Resfreq(K), 
Inphmag(K),Redfreq(K) 
DISP K 
WAIT 118.315 
NEXT K 
END SELECT 
T1=TIMEDATE 
DISP "IT TOOK ";DROUND(T1-T0,4);"SECONDS" 
WAIT 10 
OUTPUT KBD;Clear$; ! Clears the CRT 
BEEP 
ON KEY 0 LABEL “Post Process" GOTO 7650 


7640 Waitl: GOTO Waitl 
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7650 


RETURN 


TEE) EAA AAA AAA H SARA AI AHHH AAA AISA IAA AHA SIA ASIA AAS AIA ASE TA 


7670 Post_process: ! Permits user to extract plots, tables & other info 
TORQ RH AAA RA HAH AA HAH AAAI SA AAA AA LAH IH HR AIA IAA IAA HHA SEA INAH IH HII 


7690 
7700 


7710 
7720 
7730 
7740 
7750 
7760 
7770 
7780! 
7790 
7800 
7810 
7820 
7830 
7840 
7850! 
7860 
7870 
7880 
7890! 
7900 
7910 
7920 
7930 
7940 
7950 
7960 
7970 
7980 
7990 
8000 
8010 
8020 
8030 


GCLEAR 
CALL Least_squares(5*Arg1 ,Thermistor(*),Resfreq(*), 
AB,CDEFG) 

Slope=A _! Keeps call to Least_squares to one program line 
Intercept=B 
Correlation=C 

Slope_error=D 

Intrepterr=E 

Tmean=F 
Fmean=G 


Actual_temp_min=MIN(Thermmistor(*)) ! Find max and min 
Actual_temp_max=MAX(Thermistor(*)) ' temp and frequency 
Actual_freq_min=MIN(Resfreq(*)) 
Actual_freq_max=MAX(Resfreq(*)) 
Actual_mag_min=MIN(Inphmag(*)) 
Actual_mag_max=MAX(Inphmag(*)) 


Volume=.25*(Length/100)*Pi*(Diameter/100)42 
GOSUB Percent_moduli 
GOSUB Compute_modulus 


OUTPUT KBD;Clear$; ! Clear the CRT 
OUTPUT KBD;Home$; ! Home display 
GCLEAR 
PRINT TABXY(1,1) ! Start at top with blank line 
Center=(Screen-42)/2 ! Leading spaces for centering 
PRINT TAB(Center);"Key Purpose” 


PRINT TAB(Center);"--------------------------------------- " 
PRINT TAB(Center);" 0 Plot Frequency vs Time” 
PRINT TAB(Center);" 1 Plot Temperature vs Time” 
PRINT TAB(Center);" 2 Plot Frequency vs Temperature” 
SELECT UPC$(Block3$) 
CASE ="TORSIONAL" 
PRINT TAB(Center);" 3 Plot Shear modulus vs Temp” 
CASE ="FLEXURAL" 
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8040 
8050 
8060 
8070 
8080 
8090 
8100 
8110 
8120 
8130 
8140 
8150! 
8160 
8170 
8180 
8190 
8200 
8210 
8220 
8230 
8240 
8250 
8260 
8270 
8280 
8290 
8300 
8310! 
8320 
8330! 
8340 
8350 
8360 
8370! 


PRINT TAB(Center);" 3 Plot Young's modulus vs Temp” 
CASE ="LONGITUDINAL" 
PRINT TAB(Center);" 3 Plot Young's modulus vs Temp" 
END SELECT 
PRINT TAB(Center);" 4 Send Information to the Printer” 
PRINT TAB(Center);" 5 Select the graph output device” 
PRINT TAB(Center);" 6 Plot In phase mag vs Temperature” 
PRINT TAB(Center);" 7 Plot In phase mag vs Temperature" 
PRINT TAB(Center);" 7 Select the type of X axis scaling” 
PRINT TAB(Center);" 8 Select the type of Y axis scaling” 
PRINT TAB(Center);" 9 Exit this program” 


ON KEY 0 LABEL "Freq vs Time" GOSUB Freq_time 
ON KEY 1 LABEL "Temp vs Time” GOSUB Temp_time 
ON KEY 2 LABEL "Freq vs Temp" GOSUB Freq_temp 
SELECT UPC$(Block3$) 
CASE ="TORSIONAL" 
ON KEY 3 LABEL "G vs Temp" GOSUB Shear_temp 
CASE ="FLEXURAL" 
ON KEY 3 LABEL "E vs Temp" GOSUB Young_temp 
CASE ="LONGITUDINAL" 
ON KEY 3 LABEL "E vs Temp” GOSUB Young_temp 
CASE ELSE 
ON KEY 3 LABEL "" GOSUB Commentl 
END SELECT 
ON KEY 4 LABEL "Print Info" GOSUB Dump_info 
ON KEY 5 LABEL "Output Device" GOSUB Output_device 
ON KEY 6 LABEL "Grid Option” GOSUB Grid_option 
ON KEY 6 LABEL "Imag vs Temp" GOSUB Mag_temp 
ON KEY 7 LABEL "X scale option” GOSUB Xscale_option 
ON KEY 7 LABEL "Imag vs Freq" GOSUB Mag_freq 
ON KEY 8 LABEL "Y scale option” GOSUB Yscale_option 
ON KEY 9 LABEL "Exit Program” GOSUB Program_end 


8380 Blink3: WAIT 1 


8390 
8400 
8410 
8420 


DISP "Make a DECISION” 
WAIT 1 

DISP 

GOTO Blink3 
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8430 RETURN 


910 PHBH ededddadcdcdanced edad aad aaa ag aaa AE CAAA AAAI Ag Aa aaa 


8450 Grid_option: ! Accepts the user's choice for the plot grid 


8460! HAH HCA HAH A HA HH HH AH A HH AHH HH HHA A HAA HAA HK HH AH HH HC HC HAH A HH IK ek a CK 


8470 GRAPHICS OFF ! Turns off graphics display 

8480 OUTPUT KBD;Clear$; ! Clears the CRT 

8490 DISP Grid_type$;" grid is currently selected” 

8500 INPUT "Enter F - full; P - partial; N - No grid 
2" Response$ 

8510 SELECT Response$ 

8520 CASE "F" 

8530 Grid_type$="Full” 

8540 CASE "P" 

8550 Grid_type$="Partial” 

8560 CASE "N" 

8570 Grid_type$="No" 

8580 CASE ELSE 

8590 DISP "No change" 

8600 WAIT 1 

8610 DISP 

8620 WAIT 1 

8630 END SELECT 

8640 DISP Grid_type$;” grid is selected" 

8650 WAIT 1 

8660 DISP 


8670 RETURN 


165.10 Galassi 


8690 Dump_info: ! Sends selected data and table info to the printer 


8700 a A IE II TIAA TAR AAA ATTA A AA TTA STAI TIT T AT TIAA SAA AT AISA A 


8710 PRINTER IS 701 

8720 Perfskip$=CHR$(27)&CHR$(38)&CHR$(108)&CHR$(49)&CHR$(76) 
8730 Formfeed$=CHR$(12) 

8740 PRINT Perfskip$ ! Skip on Perforation 

8750 PRINT Formfeed$ 

8760 PRINT USING "3/,#" ! Three line feeds 

8770 PRINT “Filename (Rod ID): "&Block1$;TAB(70);"Page 1 of 2" 
8780 PRINT "Run: "&Block2$ 

8790 PRINT "Mode: "&Block3$ 

8800 PRINT "Date: "&Block4$ 

8810 PRINT USING "3/,7#" ! Three more line feeds 
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8820! 
8830 
8840 
8850 
8860 
8870 
8880 
8890 
8900 
8910 
8920 
8930 
8940! 
8950 
8960 
8970 
8980 
8990 
9000 
9010 
9020 
9030 
9040 


9050 


9060 
9070 
9080 
9090 
9100 
9110 
9120 
9130 


9140 
9150 
9160 
9170 


C1$="Time" 
C2$="Temperature” 
C3$="Frequency" 
C4a$="Shear modulus” 
C4b$="Young's modulus” 
C5$="In Phase Component" 
C6$="minutes” 


C7$="deg C" 
C8$=" Hz" 
C9$="Pa*10*" 
C10$="DCV" 


SELECT UPCS$(Block3$) 


CASE "TORSIONAL" 
PRINT USING 8980;C1$,C2$,C3$,C4a$,C5$ 
IMAGE 4A,6X,11A,5X,9A,5X,13A,5X,18A 
PRINT 
PRINT USING 9010;C6$,C7$,C8$,C9$,Power,C10$ 
IMAGE 5A,8X,5A,11X,2A,12X,6A,ZZ,14X,3A,12X 
PRINT 
FOR I=0 TO 5*Arg] STEP Arg1/2 
PRINT USING 9050;Time(I), Thermistor(I),Resfreq(), 
Gtor(I),Inphmag(I) 
IMAGE X,3D,8X,SDD.DDD,4X, DDDD.DDDE,9X,D.DDD,13X, 
SDD.DDDDDE 
NEXT I 
CASE "FLEXURAL" 
PRINT USING 9140;C1$,C2$,C3$,C4b$,C5$ 
PRINT 
PRINT USING 9010;C6$,C7$,C8$,C9$,Power,C10$ 
PRINT 
FOR J=0 TO 5*Arg] STEP Arg1/2 
PRINT USING 9050;Time(J),Thermistor(J),Resfreq(J), 
Eflex(J),Inphmag(1) 
IMAGE 4A,6X,11A,5X,9A,5X,15A,3X,18A 
NEXT J 
CASE "LONGITUDINAL" 
PRINT USING 9140;C1$,C2$,C3$,C4b$,C5$ 
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9180 
9190 
9200 
9210 
9220 


9230 
9240 
9250! 
9260 
9270 
9280 
9290 
9300 
9310 
9320! 
9330 
9340 
9350 
9351 


9360 
9361 
9370 


9380 
9390 


9400 
9410! 
9420 
9430 
9440 
9450 
9460 


9470 


9480 


PRINT 


PRINT USING 9010;C6$,C7$,C8$,C9$,Power,C10$ 
PRINT 


FOR K=0 TO 5*Argl STEP Argl/2 
PRINT USING 9050;Time(K),Thermistor(K),Resfreq(K), 
Elong(K),Inphmag(1) 


NEXT K 
END SELECT 
PRINT Formfeed$ ! Advance to top of next page 
PRINT USING "3/,#” ! Three more line feeds 


PRINT "Filename (Rod ID): "&Block1$;TAB(70);"Page 2 of 2” 
PRINT "Run: "&Block2$ 

PRINT "Mode: "&Block3$ 

PRINT "Date: "&Block4$ 


PRINT USING "3/,#" ! Three more line feeds 
PRINT "Physical properties:" 
PRINT 
PRINT USING "3X,12A,14X,DDD.DDD";"Total mass, grams:", 
Totalmass 


PRINT USING "3X,12A,14X,DDD.DDD";"Mass, grams:",Mass 

PRINT USING "3X,20A,06X,DDD.DDD";"Length, cm:",Length 

PRINT USING "3X,20A,06X,DDD.DDD";"Effective Length, cm:", 
Leff 

PRINT USING "3X,22A,04X,DDD.DDD”";"Diameter, cm:" Diameter 

PRINT USING "3X,21A,07X,D.DDDE";"Volume, cubic meters:”, 
Volume 

PRINT USING "3X,16A,09X,DDDD.D";"Density, kg/m43:” Density 


PRINT USING "3/,#" ! Three more line feeds 
PRINT “Least-squares fit results {freq vs temp]:” 
PRINT 


PRINT USING "3X,12A,13X,DDDD.DDD";"Slope, Hz/C:" Slope 

PRINT USING "3X,18A,07X,DDDD.DDD";"Slope error, Hz/C:", 
Slope_error 

PRINT USING "3X,14A,11X,DDDD.DDD";"Intercept, Hz:", 
Intercept 

PRINT USING "3X,20A,05X,DDDD.DDD";" Intercept error, Hz:", 


13a. 





9490 


9500 
9510 
9520 
9530 
9540 
9550 
9560 
9570 
9580 
9590 
9600 
9610 
9620 
9630 
9640 
9650 
9660 
9670 
9680 
9690 
9700 
9710 
9720 
9730 
9740 
9750 
9760 
9770 
9780 


9790 


9800 
9810 
9820 
9830 
9840 


Intrcpterr 
PRINT USING "3X,12A,13X,DDDD.DDDDD";"Correlation:”, 
Correlation 
PRINT USING "3X,20A,05X,DDDD.DDD";"Mean temp, C:",Tmean 
PRINT USING "3X,19A,06X,DDDD.DDD";"Mean freq, Hz:",Fmean 
PRINT 
PRINT 
PRINT "Other statistics:" 
PRINT 
Tmin=Actual_temp_min 
Tmax=Actual_temp_max 
Tave=(Tmax+Tmin)/2 
Fmin=Actual_freq_min 
Fmax=Actual_freq_max 
Fave=(Fmax+Fmin)/2 
PRINT USING "3X,23A,02X,DDDD.DDD";"Min temp, C:",Tmin 
PRINT USING "3X,23A,02X,DDDD.DDD";" Average temp, C:",Tave 
PRINT USING "3X,23A,02X,DDDD.DDD";"Maximum temp, C:",Tmax 
PRINT 
PRINT USING "3X,22A,03X,DDDD.DDD";"Min frequency, Hz:",Fmin 
PRINT USING "3X,22A,03X,DDDD.DDD";"Ave frequency, Hz:",Fave 
PRINT USING "3X,22A,03X,DDDD.DDD";"Max frequency, Hz:",Fmax 
PRINT 
L1$="Percent change in shear modulus per degree C: 
L2$="Percent change in Young's modulus per degree C:" 
L3$="Two point max-min approach:” 
L4$="Multi-point least-squares approach:" 
SELECT UPCS$(Block3$) 
CASE "TORSIONAL" 
Gmin=Scale_gtor_min 
Gmax=Scale_gtor_max 
PRINT USING "3X,30A,ZZ,5X,D.DDD";"Min "&C4a$, 
Power,Gmin 
PRINT USING "3X,30A,ZZ,5X,D.DDD";"Maximum "&C4a$, 
Power,Gmax 
PRINT USING "3/,#" ! Three more line feeds 
PRINT L1$ 
PRINT 
PRINT USING "3X,27A,8X,DDD.DDD";L3$,Prent_per_c_2pt 
PRINT USING "3X,35A,DDD.DDD";L4$,Prent_per_c_ls 
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9850 CASE "FLEXURAL" 


9860 Emin=Scale_eflex_min 

9870 Emax=Scale_eflex_max 

9880 PRINT USING "3X,30A,ZZ,5X,D.DDD";"Minimum "&C4b$, 
Power,Emin 

9890 PRINT USING "3X,30A,ZZ,5X,D.DDD";"Maximum "&C4b$, 
Power,Emax 

9900 PRINT USING "3/,#" ! Three more line feeds 

9910 PRINT L2$ 

9920 PRINT 

9930 PRINT USING "3X,27A,8X,DDD.DDD";L3$,Prent_per_c_2pt 

9940 PRINT USING "3X,35A,DDD.DDD";L4$,Prent_per_c_ls 

9950 CASE "LONGITUDINAL" 

9960 Emin=Scale_elong_min 

9970 Emax=Scale_elong_max 

9980 PRINT USING "3X,32A,ZZ,5X,D.DDD";"Min "&C4b$, 

Power,Emin 

9990 PRINT USING "3X,32A,ZZ,5X,D.DDD";"Maximum "&C4b$, 
Power,Emax 

10000 PRINT USING "3/,#" ! Three more line feeds 

10010 PRINT L2$ 

10020 PRINT 

10030 PRINT USING "3X,27A,8X,DDD.DDD";L3$,Prent_per_c_2pt 

10040 PRINT USING "3X,35A,DDD.DDD";L4$,Prcnt_per_c_ls 

10050 END SELECT 

10060 PRINT Formfeed$ ! Advance to top of next page 


10070 PRINTER IS CRT 

10080 RETURN 

LOCO QE I HH I AAA HAA AA AH IAA HAR AA HAA AAA AA HA ASHI IH AHI 
10100 Freq_time:! Produces frequency vs time graph w/curve 

101 1) a a a I A IIA AAR AH AAA IRA AAA A IR AIH AR HA AAAI IA IRA 
10120 GOSUB Draw_freq_time 

10130 CALL Generic_curve(Pen3,Time(*),Resfreq(*),5*Arg1) 
10140 RETURN 


STO) ot) bataaiohatatatattaaatatateataabatatahtceaeaeaeaaiiaaeaeaiinaiaaaaeedeeaiaiaaiacediaaseatataeaia 


10160 Output_device: ! Permits user to route graphs to screen or plotter 
STON G10) lasatneailstataealataaaaaataiaaiaeaiaieieaaaalateaieataaaaoataiaaiataiailasiaedeaiediidacbaaataeatsiaa 


10180 OUTPUT KBD;Clear$; ! Clear the CRT 
10190 OUTPUT KBD;Home$; '! Home display 
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10200 
10210 
10220 
10230 


10240 
10250 
10260 
10270 
10280 
10290 
10300 
10310 
10320 
10330 
10340 


10350 
10360 
10370 
10380 
10390 
10400 
10410 
10420 
10430 
10440 


GRAPHICS OFF ! Tur off the graphics display 
SELECT Plot_device$ 
CASE "Plotter" 
INPUT "Graphs appear on the screen, OK? (Y/N) ", 
Response$ 


IF UPC$(Response$)<>"Y" THEN 
Plot_device$="Plotter” 
DISP "Graphs will be sent to the Plotter” 
ELSE 
DISP "Graphs will remain on the screen" 
WAIT 1 
DISP 
GRAPHICS ON 
END IF 
CASE "Plotter" 
INPUT "Graphs are sent to the plotter, OK? (Y/N) ", 
Response$ 
IF UPC$(Response$)<>"Y" THEN 
Plot_device$="Screen" 
DISP "Plots will be sent to the Screen” 
ELSE 
DISP "Plots will stay routed to the plotter” 
END IF 
END SELECT 
WAIT 1 
DISP 
RETURN 


10450! HOA ACHC ACH AC HAH AC A ACA He he HH A A He AC ACA HAC HE HC ACH AAC a HCI A i a A i i a RI AC RK AC 


10460 Xscale_option: ! Permits user to auto or manually the X axis 
10470! We AC HR A A 9 9 A A RO II IO II I II RO ACR IC IO a I I te of i fake ke ie fe aR aK eke ie koe 


10480 
10490 
10500 
10510 
10520 
10530 


10540 
10550 
10560 


OUTPUT KBD;Clear$; ! Clear the CRT 
OUTPUT KBD;Home$; ! Home display 
GRAPHICS OFF ! Turn off the graphics display 
SELECT Xscale$ 


CASE “Auto scale X” 
INPUT "X axis is automatically scaled, OK? (Y/N) ", 
Response$ 
IF Response$<>"Y" THEN 
Xscale$="User scale X” 
INPUT "Minimum X axis value ?”,Xmin_manual 
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10570 INPUT "Maximum. X axis value ?",Xmax_manual 


10580 INPUT "Major X axis increment ?”,Xmajor_manual 
10590 INPUT "Minor X axis increment ?”,Xminor_manual 
10600 DISP "The X axis will be scaled manually" 

10610 ELSE 

10620 DISP "The X axis will remain auto scaled" 

10630 WATT 1 

10640 DISP 

10650 GRAPHICS ON 

10660 END IF 

10670 CASE "User scale X" 

10680 INPUT "X axis is man scaled, OK? (Y/N) ", Response$ 
10690 IF Response$<>"Y" THEN 

10700 Xscale$="Auto scale X" 

10710 DISP "The X axis will be scaled automatically" 
10720 ELSE : 

10730 INPUT "Change manual limits ? (Y/N) ",Response$ 
10740 IF Response$="Y" THEN 

10750 INPUT "Minimum X axis value ?",Xmin_manual 
10760 INPUT "Maximum X axis value ?”,Xmax_manual 
10770 INPUT "Major X axis increment ?",Xmajor_manual 
10780 INPUT “Minor X axis increment ?",Xminor_manual 
10790 DISP "X axis will be scaled with new values" 
10800 ELSE 

10810 DISP "The X axis will remain manually scaled" 
10820 WATT 1 

10830 DISP 

10840 GRAPHICS ON 

10850 END IF 

10860 END IF 


10870 END SELECT 
10880 WAIT 1 
10890 DISP 

10900 RETURN 


N15) 0S Sasasatdsedaaeelaaiaaiaalaiaiealadeaaadaiadiaalniadiiaiedealadaaiaiadia 


10920 Yscale_option: ! Permits user to auto or manually scale the Y axis 
10930! HAHA HEHE HE HHA HAH HHH HACIA HE HAA IH HEH HEH HHH HEH HH HH A HH HH I IK III 


10940 OUTPUT KBD;Clear$; ! Clear the CRT 
10950 OUTPUT KBD;HomeS; '! Home display 
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10960 
10970 
10980 
10990 
11000 
11010 
11020 
11030 
11040 
11050 
11060 
11070 
11°80 
11090 
11100 
11110 
11120 
11130 
11140 
11150 
11160 
11170 
11180 
11190 
11200 
11210 
11220 
11230 
11240 
11250 
11260 
11270 
11280 
11290 
11300 
11310 
11320 
11330 
11340 
11350 





GRAPHICS OFF ! Turn off the graphics display 
SELECT Yscale$ 
CASE "Auto scale Y" 
INPUT "Y axis is auto scaled, OK? (Y/N) ", Response$ 
IF Response$<>"Y" THEN 
Yscale$="User scale Y" 
INPUT "Minimum Y axis value ?",Ymin_manual 
INPUT "Maximum Y axis value ?",Ymax_manual 
INPUT "Major Y axis increment ?”,Ymajor_manual 
INPUT "Minor Y axis increment ?",Yminor_manual 
DISP "The Y axis will be scaled manually” 
ELSE 
DISP "The Y axis will remain auto scaled” 
WATT 1 
DISP 
GRAPHICS ON 
END IF 
CASE "User scale Y" 
INPUT "Y axis is man scaled, OK? (Y/N) ", Response$ 
IF Response$<>"Y" THEN 
Yscale$="Auto scale Y" 
DISP "The Y axis will be scaled automatically” 
ELSE 
INPUT "Change manual limits ? (Y/N) ",Response$ 
IF Response$="Y" THEN 
INPUT "Minimum Y axis value ?",Ymin_manual 
INPUT "Maximum Y axis value ?",Ymax_manual 
INPUT "Major Y axis increment ?",Ymajor_manual 
INPUT "Minor Y axis increment ?",Yminor_manual 
DISP "Y axis will be scaled with new values" 
ELSE 
DISP "The Y axis will remain manually scaled" 
WATT 1 
DISP 
GRAPHICS ON 
END IF 
END IF 
END SELECT 
WAIT 1 
DISP 
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11360 RETURN 


L137 DER AAA AAAI AAA AA AAA AAA AAAI AAAI AAAS SSI AAAS AAS ASR IMA AEA TER 


11380 Temp_time: ! Produces Temperature vs Time graph w/curve 
LLB Q QI I AAAI AA AAA AAA AIA IAA AA AAA AAA AAAI AA AAA ISHII AAA AIA HAT 


11400 Main_title$=Label$(7) 

11410 Sub_title$=Label$(8) 

11420 X_axis_name$=Label$(3) 

11430 Y_axis_name$=Label$(5) 

11440 Xmin=0 ! Minimum time is 0 hours 

11450 Xmax=5 ! Maximum time is 5 hours 

11460 Ymin=Actual_temp_min ! Actual minimum temperature 
11470 Ymax=Actual_temp_max ! Actual maximum temperature 
11480 GOSUB Generic_plot 

11490 CALL Generic_curve(Pen3,Time(*), Thermistor(*),5*Arg1) 


11500 RETURN 


TLS 1 01 9 9 2 ee te ek tk ok A A AA AA IAA AA AAA AAA A AAA AIH AA HAHAH AAA AA AAA 


11520 Freq_temp: ! Produces Frequency versus temperature graph w/curve 
L153 01 a AA I A AI AA AAA AAA A AAA AAAI AAA IAAT AHI I HH IAI IAA AHHH 
11540 Main_title$=Label$(6) 

11550 Sub_titleS="" 

11560 X_axis_name$=Label$(5) 

11570 Y_axis_name$=Label$(4) 

11580 Xmin=Actual_temp_min ! Actual minimum temperature 
11590 Xmax=Actual_temp_max ! Actual maximum temperature 
11600 Ymin=Actual_freq_min ! Actual minimum frequency 
11610 Ymax=Actual_freq_max ! Actual maximum frequency 
11620 GOSUB Generic_plot 

11630 CALL Generic_curve(Pen3,Thermistor(*),Resfreq(*),5*Arg1) 
11640! PEN Pen4 

11650! MOVE Actual_temp_min,Actual_temp_min*Slope+tIntercept 
11660! DRAW Actual_temp_max,Actual_temp_max*Slope+Intercept 


11670! PEN 0 

11680 RETURN 

LL 6 QQ A I AAA AIA A AAA AAA SATA AA AAA TAA AAA SHA ASA IH AAA TAA H HERA 
11700 Create_new_file: ! Creates a file for a new data mm 

LT 1 QA I a I I I IAA AAA AAA AAA AAA SAA AAR SEAR AES AERA AAA HR, 
11720 Pass$="Initial” 

11730 Grid_type$=Default_grid$ 

11740 Xscale$="User scale X” 
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11750 
11760 
11770 
11780 
11790 
11800 
11810 
11820 
11830 
11840 
11850 
11860 
11870 
11880 
11890 
11900 
11910 


11920 
11930 
11940 


Yscale$="User scale Y" 
Plot_device$="Screen" 
BEEP 1000,.1 
ON ERROR GOTO Fix2 
INPUT "Filename to store the data ?" FilenameS 
CREATE BDAT Filename$,50 
ASSIGN @Path_1 TO Filename$ 
OFF ERROR 
OUTPUT @Path_1;Block1$,Block2$,Block3$,Block4$ 
OUTPUT @Path_1;Mass,Length,Diameter,Density 
BEEP 1000,.1 
DISP "Output will be stored under Filename: ";Filename$ 
WAIT 1 
DISP 
IF Flag=1 THEN 
FOR I=0 TO 5*Argl 
OUTPUT @Path_1;Mod(I),Thermistor(I) ,Resfreq(1), 
Inphmag(I),Redfreq(1) 


11 OS 01 3 9 2 2 AR IAA IIH A IH AA ATA ATA ASIA A IAAL AI AAA IANA IAAI AK 


11960 Open_old_file: ! Opens and retrieves data from an existing file 
1197 Qt a ee a IC A AAA AA AAA AA AAA AA AIH AH AAA AAAI AAAS AAAS IAAT AAA 


11980 
11990 
12000 
12010 
12020 
12030 
12040 
12050 
12060 
12070 
12080 
12090 
12100 
12110 
12120 


Pass$="Follow on” 
Grid_type$=Default_grid$ 
Xscale$="Auto scale X" 
Yscale$="Auto scale Y" 
Plot_device$="Screen" 
BEEP 1000,.1 
ON ERROR GOTO Fix1 
INPUT "Filename to retrieve the data ?",Filename$ 
ASSIGN @Path_1 TO Filename$ 
ENTER @Path_1;Block1$,Block2$,Block3$,Block4$ 
ENTER @Path_1;Mass,Length,Diameter,Density 
BEEP 1000,.1 
DISP “Retrieving data stored under Filename: ";Filename$ 
FOR I=0 TO 5*Argl 
ENTER @Path_1;Mod(I),Thermistor(I),Resfreq(I), 
Inphmag(I),Redfreq(1) 
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12130 
12140 
12150 
12160 
12170 


Time(I)=I/Argl 
NEXT I] 
OFF ERROR 
ASSIGN @Path_1 TO * 
RETURN 


SUS :1 0) aiattahacaieieaddaiaeeiaslaiaiadiaiadeaiedieraceciaaininiesaiedadeiiaiia 


12190 Fixl: ! 


12200! HEA ACH Re AACA AH AH HA AC A A HCH He HH A AA HA HICH HH I IIH A I ACI A II I I I AI I IC 


12210 
12220 
12230 
12240 
12250 
12260 
12270 
12280 
12290 
12300 
12310 
12320 
12330 
12340 
12350 
12360 


SELECT ERRN 
CASE 53 
DISP "Limit file names to 10 characters. No punctuation” 
CASE 56 . 
DISP "This file doesn't exist on the data disk” 
CASE 58 
DISP "This file is not a BDAT file” 
CASE ELSE 
DISP " This file can not be processed by RT20B” 
END SELECT 
PRINTER IS CRT 
BEEP 300,.1 
WAIT 1 
DISP 
WAIT 1 
GOTO Open_old_file 


123.70 Pa R REAR RR ERGO Ra TOIT AIT Aa ATT ITA TAT AT TT TT IATA TTT 


12380 Fix2: ! 


12390 Fe eee eaddadddddad eed a daa aaa ROR R ARAKI AAA AI AT ATTRA ATI AIA 


12400 
12410 
12420 
12430 
12440 
12450 
12460 
12470 
12480 
12490 
12500 


SELECT ERRN 
CASE 53 
DISP "Limit file names to 10 characters. No punctuation” 
CASE 54 
DISP “Duplicate filename! Try another name.” 
END SELECT 
BEEP 300,.1 
WAIT 1 
DISP 
WAIT 1 
GOTO Create_new_file 


125.0 Fee HaHa dd dee dd aa aaa CaCO RACE TERA AC CAAT AAR T IAA AK 


1.39 





12520 Close_file: 'Closes new data file after data collection is over 
1253 Qa BB AR pepe die ioc iad dog ga aaa Go GRA AGETESA AA A AEA ATT RAE FER Fok 


12540 GCLEAR 


12550 ASSIGN @Path_1 TO * 

12560 BEEP 1000,.1 

12570 DISP "Data is stored under Filename: ";Filename$ 
12580 WATT 1 

12590 DISP 


12600 RETURN 
12.6] 0 BERRA R ae iad idaidei acerca ea daa ia ci ai aa I IA a ISA I AIA AAR AR AAA AK 
12620 Program_end: ! Shuts down shop, plays a little melody 


1.26.3 008 2 9 2 tee eae hh A A AAA AA AAI AA IHR AAA A AIHA AA AHA AA AA AIHA TAT 


12640 OUTPUT KBD;Clear$; 
12650 GRAPHICS OFF 


12660 CONTROL 1,12;0 
12670 BEEP 157,.1 
12680 BEEP 201,.1 
12690 BEEP 178,.1 
12700 BEEP 272,.1 
12710 WAIT .5 

12720 BEEP 272,.1 
12730 BEEP 178,.1 
12740 BEEP 157,.1 
12750 BEEP 201,.1 
12760 DISP 

12770 DISP “Press RUN when you are ready for another try ....” 
12780 END 


1079 0 pesados a icdidcniddeciescdddde td acdc cg go ic Ga c ag ACG CaCI 


12800 SUB Label(Csize,Asp_ratio,Ldir,Lorg,Pen,Xpos, Ypos,Text$) 
1281.0) 989 2 9 2 kk 2 a A ACAI HR RATA AIA AR HAA AA AAA AAA AAA AAR 
12820! This sub defines systems variables (Csize, LDIR, etc.), 
12830! and labels the text (if any) accordingly. 

12840 DEG 

12850 CSIZE Csize,Asp_ratio 

12860 LDIR Ldir 

12870 LORG Lorg 

12880 PEN Pen 

12890 MOVE Xpos, Ypos 

12900 IF Text$<>"" THEN LABEL USING "#,K";Text$ 
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12910 PENUP 
12920 SUBEND 


1293.0 2 A BR a a EG A IIIT TA Ta a a aT Ta TR TAI ATA ATA 


12940 SUB Lbi_axes(Csize,Asp_ratio,Pen,Xmin,Xmax,Xstep, Ymin, Ymax, Ystep) 
{DOS Qe a a ee A EA AA AHHH AA HHA AHH RAHA HAHAHAHA AHF ALAA IIA 
12960 DEG 

12970 CSIZE Csize,Asp_ratio 

12980 PEN Pen 

12990 CLIP OFF 

13000 LDIR 0 

13010 LORG 6 

13020 Yrange=Ymax-Ymin 

13030 Yoffset=.02*Yrange 

13040 FOR L=Xmin TO Xmax STEP Xstep 


13050 MOVE L,Ymin-Yoffset 
13060 IF ABS(L)<.001 THEN L=0 
13070 LABEL USING "#,K";L 
13080 NEXTL 


13090 LORG 8 
13100 Xoffset=.02*(Xmax-Xmin) 
13110 ~~‘ IF Yrange<=1 THEN 


13120 Mmax=DROUND(Yrange/Ystep, 1) 
13130 Yval=Ymin 

13140 FOR M=0 TO Mmax 

13150 IF ABS(Yval)<=.001 THEN Yval=0 
13160 MOVE Xmin-Xoffset, Yval 

13170 LABEL USING "#,SD.DD"; Yval 
13180 Yval=Yval+Ystep 

13190 NEXT M 

13200 ELSE 

13210 FOR M=Ymin TO Ymax STEP Ystep 
13220 IF ABS(M)<=.001 THEN M=0 
13230 MOVE Xmin-Xoffset,M 

13240 LABEL USING "#,K";M 


13250 NEXT M 
13260 ENDIF 
13270 CLIP ON 
13280 PENUP 
13290 SUBEND 


BS OQ | a a a ee eh a A AA I FI IH AIA AR I IR AAA HR AO AE IIR TT TE IO 
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13310 SUB Block_info(Block1$,Block2$,Block3$,Block4$,Xgdumax,Pen,Pen2) 

13320 9a A I I AA IA IIA AH AI RIAA AAA AA AA AAA AAA AA AAI AAAI ATA TR TA 

13330 CALL Label(3,.6,0,2,Pen,2,2,"Rod: "&Block1$) 

13340 CALL Label(3,.6,0,5,Pen,Xgdumax/3,2,"Run: "&Block2$) 

13350 CALL Label(3,.6,0,5,Pen,Xgdumax*2/3,2,"Mode: "&Block3$) 

13360 CALL Label(3,.6,0,8,Pen,.97*X gdumax,2,Block4$) 

13370 MOVE 0,4 

13380 PEN Pen2 

13390 DRAW 133,4 

13400 PENUP 

13410 SUBEND 

34.2.0)! a A AA IA A AA AAA AHA HAH AR AA AAA AAA A AH AA ASIA ASIA IAA ASAE 

13430 SUB Least_squares(Imax,X(*), Y(*),Slp,Int,Cor,SIp_er,Int_er, 
Xmean, Ymean) 

[BGAN ea 2 ea 2 ek ke I I I AIR AAI HAHA AR A HAHA AA AA HAHA AHA AAA HA IIA AAA. 

13450 DISP "Computing least-squares fit ......... ’ 

13460 Sumx=0 

13470 Sumy=0 

13480 Sumxx=0 

13490 Sumxy=0 

13500 FOR I=0 TO Imax 


13510 Sumx=Sumx+X(I) 

13520 Sumxx=Sumxx+X(I)42 

13530 Sumy=Sumy+Y(I) 

13540 Sumxy=Sumxy+X(I)*¥(1) 

13550 NEXTI 

13560 Delta=(Imax+1)*Sumxx-Sumx42 

13570 Int=(Sumxx*Sumy-Sumx*Sumxy)/Delta 
13580 Slp=((Imax+1)*Sumxy-Sumx*Sumy)/Delta 


13590 Sumerrerr=0 
13600 FOR J=0 TO Imax 


13610 Sumerrerr=Sumerrerr+(Y (J)-Int-Slp*X(J))42 
13620. NEXT J 

13630 Sigmayy=Sumerrerr/(Imax+1-2) 

13640 Sigmay=SQR(Sigmayy) 

13650 Int_er=Sigmay*SQR(Sumxx/Delta) 


13660 Slp_er=Sigmay*S QR((Imax+1)/Delta) 
13670 Xbarsum=0 
13680 Ybarsum=0 
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13690 FOR K=0 TO Imax 


13700 Xbarsum=Xbarsum+X(K) 
13710 Ybarsum=Ybarsum+ Y(K) 
13720 NEXT K 

13730 Xmean=Xbarsum/(Imax+1) 
13740 Ymean=Ybarsum/(Imax+1) 


13750 Sigmaxxsum=0 
13760 Sigmayysum=0 
13770 Sigmaxysum=0 
13780 FOR L=O TO Imax 


13790 Sigmaxxsum=Sigmaxxsum+(X(L)-Xmean)42 

13800 Sigmayysum=Sigmayysum+(Y(L)-Ymean)*2 

13810 Sigmaxysum=Sigmaxysum+(X(L)-Xmean)*(Y(L)-Ymean) 
13820 NEXTL 

13830 Sigmax=S QR (Sigmaxxsum/(Imax+1)) 

13840 Sigmay1=SQR(Sigmayysum/{Imax+1)) 

13850 Sigmaxy=Sigmaxysum/(Imax+1) 

13860 Cor=Sigmaxy/(Sigmax*Sigmay 1) 

13870 = -DISP 


13880 SUBEND 


13 BOC A ee ee eke kk I a A AI I AA II IH IIA AAAI A AA HAHA AK AAA HAA AA AAA AAA IAA A. 
13900 SUB Scale(L,R,B,T,Top,Left,Xcenter, Ycenter,X gdumax, Y gdumax) 
1391.0) 9 9a ee ek ee kk eh a HH HH A AAA AA AIA AR AAA AHH HARTA A AH AA 
13920 Top=T*Ygdumax 

13930 Bottom=B* Ygdumax 

13940 Left=L*Xgdumax 

13950 Right=R*Xgdumax 

13960 Xcenter=(Right+Left)/2 

13970 Ycenter=(Top+Bottom)/2 

13980 VIEWPORT Left,Right,Bottom,Top 

13990 SUBEND 

LAI 1 9 a a ah a a A AE I IR AA A AIH A AAA AH A IER A IA A IRA AIA IAA I IH 
14010 SUB Yscale(Ymin, Ymax, Yminor, Ymajor) 

JQ.QD 01 A a HH HAI A AA AH I IA A AAA AH AH HI A AHA RAH AAA AA AA TAA TA ATH AA ACA H TR. 
14030 DIM Diff(36),Minor(36) 

14040 DATA .1,.005,.2,.01,.25,.01,.3,.02,.4,.02,.5,.02,.75,.05,1,.05 


14050 DATA §,.2,10,.5,15,1,20,1,25,1,30,2,40,2,50,2,75,5,100,5,125,5 

14060 DATA 150,10,200,10,250,10,300,20,400,20,500,20,750,50,1000,50 
14070 DATA 1250,50,1500,100,2000,100,2500,100,3000,200,4000,200,5000 
14080 DATA 200,7500,500,10000,500 
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14090 FOR J=1 TO 36 

14100 READ Diff(I),Minor(1) 
14110 NEXTI 

14120 Yrange=Ymax-Ymin 
14130 Index=1 

14140 WHILE Yrange>Diff(Index) 


14150 Index=Index+1 
14160 END WHILE 
14170 Yminor=Minor(Index) 


14180 Ymajor=Diff(Index)/5 
14190 IF Ymin<0 THEN 


14200 Newmin=Ymin-Ymajor+ABS(Ymin MOD Ymajor) 
14210 ELSE 

14220 Newmin=Ymin-(Ymin MOD Ymajor) 

14230 ENDIF 

14240 Newmax=Newmin+Diff(Index) 

14250 WHILE Ymax>Newmax 

14260 Index=Index+1 

14270 Yminor=Minor(Index) 

14280 Ymajor=Diff(Index)/5 

14290 IF Ymin<0 THEN 

14300 Newmin=Ymin-Ymajor+ABS(Ymin MOD Ymajor) 
14310 ELSE 

14320 Newmin=Ymin-(Ymin MOD Ymajor) 

14330 END IF 

14340 Newmax=Newmin+Diff(Index) 


14350 END WHILE 

14360 Ymax=Newmax 

14370 Ymin=Newmin 

14380 | SUBEND 

GSO 1 He a AIR AR AH A AAAI A HH AA HAA AHA AAA AA HAHAH AA AAAI AAA AR, 
14400 SUB Xscale(Xmin,Xmax,Xminor,Xmajor) 

1G] (1 a a a I I A A A A AH AR AAA RAHA AHR AAA A IAAI AAAI IHA AHA IA AAA A 
14420 DIM Diffx(40),Minorx(40),Majorx(40) 

14430 DATA S5.,0.2,1.,6.,.2,1.,7,.2,1.0,8.0,.2,1,10,.5,2.0,12,.5 

14440 DATA 2.,14.,.5,2.0,16.,.5,2.,20.,1.0,4.0,25.0,1.0,5.0,30.0,1.0 
14450 DATA 35.0,35.0,1.0,5.0,40.,1.,5.,50.,2.0,10.,60.,2.0,10.,70.,2. 
14460 DATA 10.0,80.0,2.0,10.0,100.0,5.0,20.0,120.,5.,20.,140.,5.,20. 
14470 DATA 200.0,10.0,20.0,400.0,10.,50.,600.,20.,100.,800.,50.,100. 
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14480 DATA 1000.0,50.0,100.0,1200.0,50.0,200.0,1400.0,50.,200.,1600. 


14490 DATA 50.0,200.,1800.,100.,300.,2000.,100.,500.,2500.,100.,500. 
14500 DATA 3000.,100.,500.0,3500.,100.,500.0,4000.,200.,400.0 
14510 ~— FOR I=1 TO 34 

14520 READ Diffx(I),Minorx(!),Majorx(1) 

14530 NEXT I 


14540 Xrange=Xmax-Xmin 
14550 Index=1 
14560 WHILE Xrange>Diffx(Index) 


14570 Index=Index+1 
14580 END WHILE 
14590 Xminor=Minorx(Index) 


14600 Xmajor=Majorx(Index) 
14610 IF Xmin<0 THEN 


14620 Newmin=Xmin-Xmajor+ABS(Xmin MOD Xmajor) 
14630 ELSE 

14640 Newmin=Xmin-(Xmin MOD Xmajor) 

14650 ENDIF 

14660 Newmax=Newmin+Diffx(Index) 

14670 WHILE Xmax>Newmax 

14680 Index=Index+1 

14690 Xminor=Minorx(Index) 

14700 Xmajor=Majorx(Index) 

14710 IF Xmin<0 THEN 

14720 Newmin=Xmin-Xmajor+ABS(Xmin MOD Xmajor) 
14730 ELSE 

14740 Newmin=Xmin-(Xmin MOD Xmajor) 

14750 END IF 

14760 Newmax=Newmin+Diffx(Index) 


14770 END WHILE 

14780 Xmax=Newmax 

14790 Xmin=Newmin 

14800 SUBEND 

148] 0 pesee eee cacaaadeddcoiacciacacde a dae de agg gaara CAAA AIA A 
14820 SUB Generic_curve(Pen,X(*), Y(*),Max_point) 

1G G3 01 a ee a A A A AC A AA AHA AAA AAA AA AAA TH AAA AHA HHA AT TAHA AAA HE 
14840 PEN Pen 

14850 FOR Point=0 TO Max_point 

14860 PLOT X(Point),Y(Point) 

14870 NEXT Point 
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14880 PENUP 

14890 SUBEND 

AQ A I AA A AAA AAI ARIA IIIA III HII EH AIRSET AAA AA TE HE 
14910! Nothing follows. 

AQ) a He de Ae A AAA HAHAH AAA HH HHA AHA AHIR HH AAA HE HAA ARH HK AH KAMA A HHH FR. 
14930 Pessoa iaciootcectciac oe ag aaa ag AACE AIA IAA AAA AI IH TH Ik 
AQMD A HAHA IA AHA AA AHHH AIM HAHAHAHAHA AAA ALIAS HAHAH IIH AH REA IA HAR FO 
14950 SUB Calc_redfreq(Imax,Thermistor(*),Resfreq(*),Redfreq(*)) 


14960! HEA eH HHH He He HICH AIA HAHA AHH KAA HHA HAA HAA HAH HH HAHA HH HEH HHH HHA HHH A HH AK 


14970 FOR I[=0 TO Imax 


14980 Temp(I)=Thenmistor(I)+273.15 

14990 X(I)=-129*(Temp(I)-283.15)/(107+Temp(I)-283.15) 
15000 Y(I)=X(1)*LOG(10) 

15010 Alphat(I)=EXP(Y(I)) 

15020 Redfreq(I)=Resfreq(I)* Alphat(1) 

15030 NEXT I 


15040 SUBEND 
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APPENDIX E. PR-1592 DATA 
CONTENTS: 
- PR1592 FIRST TORSIONAL MODE DATA (TABLE E.1) 
- PR1592 SECOND TORSIONAL MODE DATA (TABLE E.2) 


- PR-1592 DIRECT Q MEASUREMENT DATA AT VARIOUS TEMPERATURES 
FOR CALIBRATION OF THE IN-PHASE VOLTAGE (FIGURES E.1-E.12) 
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TABLE E.1 PR1592 FIRST TORSIONAL MODE DATA 
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|i Label A | B | c D i 
Label SHEAR (CORREGTEMP (C) IRES FREQ INPHASE i 
fl Gee 2558-al -145711 73355994) -0.204615} 
2\<e 3.440E-81 -13.338| 721.3546i] -0.206158 
3[<rsar ,a% 3.340E-8] -12.741] _710.72646| _-0.203729 
at 5. 249e-8) = 12.262] 700.97543]__ -0.201988 
s{[ TS 7SE-8] = 11.877] 692.9811] -0.201168 
6 3.11S5E-8[ -11.553] 686.58914| -0.201192} 
7 
8 2.971E-8 -10.518|  670.38667| -0.200194 
9 
10 651.01918| _-0.2004Sél 
11 2.71 8E-8] -8.76}  641.19772 -0.201 166} 
12 2.634E-8 -8.145] 631.22634| _-0.198666! 
13 2.55 4E-81 -7.596] 62154814] _-0.196348 
i4f sd 2.479E-8! -7.095] 612.32536| -0.195655} 
i el 2.399E-3i -6.493|  602.37463| -0.196473' 
16 2.30 i£-a] -5.901| _589.97692| _-0.200562I’ 
17{ 2.2298 +8 -5.179| $80.62669| -0 209166 
18 2.158E-8l -4.545] 572.65712| -0.211506| 
19 2.11 2E-8) -3.963] $65.16869| -0.207602} 
20 2.013€-8i -3.243] 551.78365| __-0.204488 
Di 1.944E-a! -2.588! $42.26711|. -0.206024} 
22 1. B75E-8! -1.99] _532.57848 -0.20645] 
233. —— — = 1,803E-8i =1.232) 522.23082| _-0.207881 
24 1.734E- 38) -0.618) 512.10938| -0.209376I 
25 1.668E-a! -0.009] _502.29582| _-0.213031 
26 1.61 4E-8/ 0.627; 49409934] _-0.216686 
27 1.558E-8l 1.423] 485.46042] __-0.218909 
28 1.499E-8l 2.082} 476.17158| -0.220244| 
29 1. 445E-81 2.734| _46756027| _-9.222705 
30 1.395E-8l 3.37] 459.405| _-0.225057 
Zr 1.350E-6] 4.032} 451.83432| _-0.227289} 
32 1.305E-al 469] 444.33496| —-0.229338| 
33 1.264E-3 5.356] 437.19854 -0.22611 
34 1.216E-8] 6.01 428.83433 -0.227894 
35 1.172E-8l 6.673) _420.96193| - -0.230282} 
36 1.1 30E-81 7.339]  413.45877| _ -0.234065) 
7. — -7 — 1.089E-8} 8.002 405.9057 -0.22729 
33, 1,049E-al 8.697| 39838314] -0.242718 
39] 1.014E-al 9.3841 39158948] -0.247857 
4o[ | 9. 799-7] 10.024 384.983| _-0.252229 
aif 9.479E- 71 10.702 378.6403| _-0.2584321 
a2p | __9.156€-7| 11.364] _372.14826| _-0.262851| 
43|_ {i886 2E-7} 12.028| _366.12717| -0.268638 
44] 5.5835E-71 12.692] 360.31887| -0.2771391 
45] 3.21 4E-7} 13.333) 35460743| _-0.278625} 
46, sd 2 OS7E-7! 14.012! 749.0958] _-0.287906/ 


TABLE E.1 CONTINUED 


|S Label A | 8 | c | p) | 
47 7.613€-7| 14669] 343.75973{ _-0.290056} 
45{ sd 7.534€-7! 15.354] 338.7007} _-0.295088 
49] sid 7 SSE-7 16.01} 33353144! -0.300249] 
50 7.138E-7 16.678 328.573 -0.30658} 
5] 6.943E-71 17.339} _324.05237 -0,21346 
52 6.76 1E-7I 18] 319.78518 -0.31969 
53] 8 5 85E-7| 18.672] _715.59445 -0.32568 
S4 -0.32927 
ss) 6.256E-7! 20.003! 307.50354 -0.32995 
56 6.049E-7} 20.651] 302.47952 -0.33045) 


S7ip. ee 5.870E-7 21.312! _297.97926 -0.3426} 
sats 5.721E-7! 21.99] 294.1769 -0.35217! 


59 5.586E-7/ 22.652] _290.66802 -0.36078| 
60 23.301] _287.30684 -0.26804 
61] sd S33 4E-7} 23.959] _284,04529 =0.37485 
62 5.220€-71 24.623 280.9944 -0.38154 
63 =0.38702I 
64 5.007E-7! 25.925] _275.19598 -0.3925] 
65 4.889E-7/ 26.595] _271.94384 -0.78424} 
66 4.751E-7{ 27.248] _268.07988 -0.39235} 


67|. —s—ss—s—idE:sC “i 6 22E-7) 27.908] _264.40884] -0.4035} 
a 4.515E-7] 28.58] 261.72992 -0.41436} 
69, sd 4.420€-71 29.225} _253.57502 -0 42416} 





70 4.332E-7| 29.873| 255.98542 -0.43227 
71 4.241E-7] 30.532!  253.27751) -0.47964} 
72 4.138E-7! 31.179 250.1834) -0.45577} 


73. | 4.065E-7! 31.882| 24795498] _-0.46195} 
| nea 3.988E-7| 32.534] _245.60855| _-0.46908] 


75 3.91 7E-7 33.178]  243.41238 -0 474171 
76 3.856E-7| 33.855|_ 241.51547 ~0.47927 
77 |. 3.790E-7/ 34.508] 239.41813 -0.48234) 
78; sd 3.720E+7] 35.193 237.2095 -0.48768 
79, sd] 3.6 49-7] 35.858] _234.93002 -0.49373! 
so} | 3.571E-7! 36.521]  232.42163|  -0.50474| 


31 3.508E-7/ 37.188] _230.33902 -0.51896 
62, ssid 3.450E-7] 37.866| _228.438611 -9.53036] 
esti iQ E-7I 33534] 226.48174 0.545331 


Ec 3.336e-7} 39.194| 22463599 -0.55661! 
as{ sd] 3288-77) 39.862] 223.00513 0.567881 
a6] | 26 E-7I 40.535) 221.24825 -0.58337 
a7] | S31 94-7} 41.218] _219.79415 -0.591781 
ea 3.1S1E-7] 41.876] _213.32857 -0.60326 
go, sd NOE“? 42535] 216.8855 -0.61481} 


90 3.070E-7I 43.236] _215.47867 -9.62977} 
gif YT E-71 43.699! 21421148 -0.63979} 
92 2.998E- 7] 44562) 212.9569 -0.650661 





or] | 296 5E- 7 45.277! 21'.77474 -0.66305) 
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210.6433 
20950511 


208.3323) 
2.840E-7| 47.958] — 207.27333 
2.783E-7{ 49.284] 205.16444 
204.16496 

50.637]  203.16787 

2.705E-7! 202.27369 
2.682E-7} 51.964| 201.40724 
2.635E+7} 55.2711 199.64121 
2.513E-7 53.937)  198.79755 
















2.571£-7! 55.277| 197.1887 
109 2.549E-7} 55.947] _196.25229 -0.8.4034} 
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19413874 -0.87835 












113 2.47 4E-7! 58.626] 193.42544) -0.89325i 
1 2.45 7E-7{ 59.296] 192.76385 -9.90713] 
115 192.0869 -0.9198} 
116 2.424E-7} 50.622 191,47189 -0.97 4221 
117 2.407E+7! 61.291] 190.8043 -9.95009} 
118 2.392671 61.973] 190.21615 -0.96 199} 
119 2.376E-7} 62.65] 189.59101 -0.98015i 
120 188.9780] -0.99487 
121 2.346E-7I 64.4538 188.38921 1.01129 
122 
123 
124 2.304E-7/ 66.417| 186.69455: -1.75918] 
125 2.293E-7/ 67.073| 186.2!917 -1.07342} 
126 2.280E-7! 67.736]  185.68471 -1.09107! 
127{ | 2.26 8E- 7) 68.377] _185.21015 -1t1157i 
128 2.258E-7| 69.0311  184.78674 -1.12665 
129 2.247E-7| 69.694 184.76136 -1.143 al 
130 2.236E-7/ 70.345}  182.88961 -1.15939] 
if 225e-7f 71.007] 183.44837| =1.1749i 
32 2.21S5E-7} 71.671 183.04204 -1.1895| 
133 2.206E+7 72.344]  182.54761 -1.19659 
134 2.195E-7| 73.028] 182.21211 -1.2236 
135 - 1.24404) 
126 2.130E-7{ 78.393 179.47319 -90.64647! 
37, SY 2.121€-71 79.063] 179.11 125] -0.65615 
at 2.115E-7! 79.704) 178.84327 -0.662691 
179 | 2.107E-71 30.379) 173.51625 -9.67139 
140 2 100E-7i 31.046] 175.2167) -9.56125} 
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TABLE E.2 PR1592 SECOND TORSIONAL MODE DATA 



























Label A I B | c D | 
Label [corrected [SHEAR MODULUITEMPERATURE RESONANT FREGIN-PHASE VOLi 
1 F.5816-8! -13.765| 1479.4419| -0.204575] 
2 -13.05] _1455.3742| _-0.203569 
3 3.369E-8] -12.491] 1434.9878] -0.201524 
af 3.278E+8] -11.962| 1415.4678| -0.199589 
st 3.18 16-8 -11,335) 1394265] -0.197678l 
i | 3.093E-8] -10.843) 1374.8478 -0.1956} 
7/2. 990€-8 =10.132| 1351.8373]  -0.193.46| 
af sd 2903-8} -9.649] 1331.9333| -0.192172! 
of 2 BO5E- 8} -8.967| 1309.3382| -0.192172! 
tof | 2724-8} -8.453| 1290.23116| _-0.190769} 
es ee en | -7.74| 1268.17618] -0.189372/ 
12 2.552E-3} -7.265| 1248.96705] -0.1875931 
7 2.465E-8 -6.59] 122732184] -0.186087! 
14 2.376E-8) -5.896] 1205.01801| -0.187699} 
15 2.294E- 8} -5.238| 118412597! -0.187179} 
16 2.21S5E-3i ~4.611] 1163.43138]  -0.186778 
17 2.140E-ai -4.012] 1143.558561 _-0.186475| 

2.067E-8! 294j 1123.99536 

1.999E-3i 1105.33972 
20 1.936E-8f -2.243|  1067.7452| -0.186059 
21 1.872E+8} -1.602) 1069.73184] -0.1856711 
22 1.812E-5} -1.076] 1052.48279| -0.1856211 
22 1.75 1£-al -0.468} 1034.36234] -0.185805I 
24f si 1.688E-3} 0.179} 1015.82959| -0.187108i 
| ae 1.626E-8) 0.817] 996.8548]  -0.187935| 
26], 1.56 7E-8I 1.494] _978.62057| _ -0.189358 
27, 1.511€-8l 2.118] 961.i2765| -0.190594I 
28; id 1.45 7E-8 2.745] 94363414] -0.191595 
og]. 1.405E-8} 3.294| 926.732723| -9.192635! 
30 1.354E-3i 4.044] 909.76509| _-0.193.422) 
1.302E-8] 4.692] _892.05993 197496] 
1.260E-8] 5.342} 877.50896; -0.200118] 
1.216E-8} 861.99127 -0.201633} 
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8.880E- 7} 
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8.038E-7! 
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8 46.46.4472 
831.82277 
817.17401 
802.76867 


-0.204303} 
-9. 206897 
-0.209815] 
-0.213018} 












8.654] 











9.313} 768.80904|  -0.216215} 
9.981]  775.30585 -0.21979) 
10.628] 762.93628] -0.223594l 
11.268] 749.23633| -0.227535/ 
11,938]  736.72315{ -0.231776) 
12.605} 724.4171)  -0.276012I 
13.275] 712.49196]  -9.240624] 
13.925) 700.90944}  -90.245246! 
145691 569.7SE56! -0.250075I 
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5.01 5E-7| 
4.901E-7} 


25.174 
25.825| 


553.5055 2) 
$47.33195 















MS abel A t 8 | c ») | 
47 7.542E-7] 15.261 678.93951| -0.255227) 
48 7.307E-7] 15.931] 668.3011] -0.26059a! 
49 =0.265975 
50 -0.271756| 
51 6.669E-7} 17.918 638.4309} -0.27791 4] 
52 ~0.284022 
53 =0.289932 
34 6.124E-7 19.905] _611.80764| _-0.295973 
55 ~0.302276} 
56 5.805E-7! 21.222] 595.65898 -0.30879 
57 0.31541 
53 5.513€-7| 22.532] _580.50025 -0.32226i 
$9 
60 5.250E-7| 23.839] 566.4781 -0.33583 
61 5.123E-71 24.511] 559.86453 -0.342721 


64 4.794E- 7] 26.503] 541.28514| —_-0.3635 11] 
65 4.691E-7} 27.166] _535,45388 -0.37065] 
66 4.594E-7; 27.855] _529.90476 -0,.37776 
67 4.503E-7} 28.489} 524.59373 0.38436 
68 4.419E-7! 29.156! 519.15705 -0.291 33] 
69 4.224E-7} 29.816] 51408658 -0.39856| 
70 4.242E-7] 30.475] _509.19762| _-0,40537] 
71 | 4.163E-7} 31.15| _504.43773 -0.41212 
72 4.088E-7/ 31.8] _499.85922] _-0.41809| 
73] 4012-7} 32.4441 495.20548| —_-0.42263} 
74 3.942E-7| 33.113] _490.82007 -0.42835} 
75 3.873E-7 33.793| _486.56068 0.43654 
76 3.804E-7| 34.463] —_-482.19882| -0.44428} 
77 3.73.4E-7] 35.137] _477.74448 -0.45441 
78 3.671€-7| 35.794] 473.65325 -0.46507 
i) as ae 36.446] 469.8187|  -0.a7448 
80 2.55 4€-7| 37.106] 466.0964 -0,48454{ 
ai] ——s—s—sd~Ssi3.SO0E-7] 37.778] _462.52808 -0.49781 
82] sd Ssédz*«CAAE-7{ 38.439] . 459.02312 -0.50374 
asf 5596E-7] 39.11] 455.58934]  -0.513503) 
saf 3.348E-7| 39.76) _452.32786 -0.52216| 
as| sd i302 - 7} 40.427] _449.21158 -0.53152 
36] 7.256E-7/ 41.063] _446.12375 -0.54071 
a7]. 321-7} 41.722] 443.21414 -0.5497 
saj 3.172€-7| 42.394| 440.29983 -0.55932} 
ao] sid iS IE 43.029] 43 7.46565| 0.569401 
gof CT 091 E- 73.715] 43.4.65158]_ -0.57941] 
gif sd LOS3E-7 44.387] 431.9919 0.58942 
92; sid 45.043] _429.37615 -0.59841 
ee ee 45.693] 42659713 -9.60969 
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TABLE E.2 CONTINUED 


| sm Label A i 5 | ow a) i 
ga] 2.943E-7i 46.359] 424.13577 -0.62046] 
o[ 2.910E-71 47.024] 421.71895}  -0.63076} 
96{ 2.877E-7 47.696] _419.32634 -0.64117] 


97 2.844E-7| 48.34]  416.94878 -0.65181 
98 2.813E-7| 49.011 414.62369 -0.66187 


99] 2.782E-7| 49.655 412.3848 -0.67366 
100] —i—“‘ié‘wY:*Ci(itsi SET 50.315] 410.20276 -0.68586 


101 50.984] 408.10427 -0.69694 
1o2] sd 2.697E-7I 51.661] 406.03963 -0.71121 
103 2.671E-7i 52.32 404.0668 -0.72541} 
104 2.545E-71 52.992! _402.10598 -0.74147} 
105 2.621E-7| 53.636} 400.26866 -0.75571 
1o6f sd 2.598E-7! 54.295] 398.47694 -0.77046} 
107 2.575E-7| $4.969| _396.74903 -0.785385} 
108ft—‘“<té‘idYSs*si SS ET= 7} 55.6421 395.08846] _-0.80051} 
109 2.534E-7] 56.286] _393.54082 -0.81395} 
110 2.51 4E-7| 56.96] _391.99084 -0.8301f 
Mt 2.49 4E-7| 57.628] _390.44766 -0.84369} 
112 2.476E-7 58.277|  388.99518 -0.85936 
113 7 2.453£-7 58.949] 38758625 -0,87417 
114 2.440E-7| 59.628]  386.17226 -0.88907 
WS -0.90301 
116 2.407E-7} 60.949] 383.52576 -0.91894I 
It? 2.391€-7i 61.618] 382.25039 -0.93359 
118 2.375E-71 62.264} _381.00867 -0.94842| 
119 2.362E-7} 62.898] 379.91389) -0.96294| 








2.548E- 7} 64.088 










378.79252 -0.97671| 






































121 2.23 4E-7} 64.758] 377.68567 -0.99095| 
2.3216-7! 65.416] _376.61856 

123 66.033| 375.60847 -1.01885} 
124 2.297E-7I 66.68] _374.66032 - 1.03378] 
125 373.69562 - 1.04873] 
126 2.274E-7} 67.965 372.79373 - 1.06219} 
127 2.267€-7| 68.623} 371.8829 -1.08003} 
128 2.252E-7 69.3021 _370.96992 -1.09622 
129[.—“‘~*:SsiCi a E-F} 69.932] _370.09256 = 1.112421 
130 2.23 1E-7| 70.579] _369.24057 -1.12564] 
134 2.221€-7} 71.229} 368.41917 -1.14162] 
132 2.21 1E-7| 71.892) 367.63375 -1.153805} 
133 2.202E-7I 72.53] _366.83531 -1.17491 
134 2.192E-7 73.15]  366.05546 -1,18774] 
135 2.183E-7I 73.816] 365.28087 -1.20351| 
136 2.175E-7] 74.431] 36458518 - 1.22052 
137 2.166E-71 75.193 363.37027 -1.23791 
2.157E-71 75.855) 362.07272| - 1.254891 

76.486) 552.46716 -1.27178l 

77.147) 261.7849 - 1.28826) 
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Figure E.1 PR-1592 direct measurement of Q, T=-16.75°C. 
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Figure E.2 PR-1592 direct measurement of Q, T=-10.16°C. 
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Figure E.3 PR-1592 direct measurement of Q, T=-5.67°C. 
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Figure E.4 PR-1592 direct measurement of Q, T=4.50°C. 


158 


FREQ RESP 






























































CURVE _FIT 
=3 0 =i} 
=19 
[=j=) 
dB 
-90. 90 
-90.0 | 
aes | 
FxdxXY 200 H2 Ak 
Curve Fit 
Poles And Zeros 
POLES Oo ZEROS Oo 
pi “87 le 752 235. 404 
2 218. 909 437.171 
3B -68.7345 +)j 322.643 -207. 0S2 &)j S95.. 228 
4 -100. 198 +Jj 673.981 663.713 +) 1. 06869k 
S 285. 256 =) 1« 10642 650. 187 +) 1. 88526k 
is} -408. 192 +j 1. 70371k S20. 447 +) 2.73424k 
7 =457.86Ge +] 2. 35k BOG. 245 +j 3. 56935k 
8 —-448.82 +J 3. O341kK 141.169 + 3. S0O45k 
9 -313.949 +) 3.72105k —-S63. 038 +Jj 4.14441k 
10 183.125 +Jj 3. 80795k 











Time dalay= G.0Q S Gatin=-e2. OO3K Scalea= iG } 


nel 








Figure E.5 PR-1592 direct measurement of Q, T=14.23°C. 
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Figure E.6 PR-1592 direct measurement of Q, T=25.40°C. 
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Figure E.7 PR-1592 direct measurement of Q, T=37.00°C 
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Figure E.8 PR-1592 direct measurement of Q, T=46.42°C. 
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11 19.2138 +j 2. 81961ik 20.9659 +) 2. 82133kK 
L2 123.506 +Jj 3. 1524ik 

Tima dealay= 0.0 S Gain=-19. 014 Scala= i. 'O 





Figure E.9 PR-1592 direct measurement of Q, T=55.85°C. 
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X=179. 6 Hz 
Yor-12.513 dB 
FREQ RESP 
Yb=-12.095 dB 


CURVE FIT 
o.0 











dB 





dB 











-64.0 











-64.0 
FxdxyY 100 Hz 2k 
































Curve Fit 
Poles And Zeros 
POLES o ZEROS Oo 

1. -3. O3679k -S16. 39S 
2 -14.5652 +) 176.623 168. 584 
3 -28.939S8 +] 373.534 -145. 55S «Jj «A72. 82 
4 -—38. BSe22 2.) SSS. 168 287.662 +) 498.02 
S -70.4523 +) 784.575 B26.113 +) 942.796 
6 -53.6018 +] 982. 986 -115.065 +) 955.716 
7 -160. 787 +3} 1.1801Sk S62. 852 =) 1. BOiL24& 
8 -164.017 +) 1.47345k 266.773 +) 1. S6879k 
9 -46.2129 +) 1.6785Sk -60.6101 +3 1. 68025k 
10 -148.596 +J 1. 8797K 208. 658 +.) &. O33SSK 

Tima dalay= 0.0 S Gain=—-4.193 Scole= leo QD 





Figure E.10 PR-1592 direct measurement of Q, T=67.01°C. 
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xe174.6 Hz 


FREQ RESP 

Yb=-10. 887 dB 
CURVE_FIT 
o.0 









































0.0 
dB 
dB 
-64.0 
-64.0 
FxdxXY 100 H2 y 2k 


Curva Pit 
Poles And Zeros 


POLES ie} ZEROS a) 


1. 20563 154. 08s 
-14.1148 +) 175.266 -105.933 +) 346.865 
-17.6794 +) 358.454 224.096 +) 438.092 
-42.4391 +) 549. 382 -105.149 +) 754. 852 
-37.4546 +) 753.913 254.498 +) 856.805 
-86.5696 +} 962.877 -74, 2487 +} 1.15128 
-S2. 3113 +J 1.16796 261.346 +Jj 1.25141k 
-1.19598k+] 1.2689 120.194 +) 1.64988 
-220.569 +} 1.3634k 226.447 +} 1.69923 
10 -297.246 +) 1.61223 201.572 +j 2.01947k 
11 112.463 +) 1.66329 
12 -215.929 +) 1.8707 









ODONOUSKWNY 





| 27 uaS dalay= 0.0 S Gain=27. 4E+9 Scala= 1.0 








Figure E.11 PR-1592 direct measurement of Q, T=74.60°C. 


165 





X=169.6 Hz 


FREQ RESP 

Yb=-10.549 dB 
CURVE FIT 
o.0 





o.0 


~ FAIA 


o aie a 





—— 


-64.0 





FxdXY 100 Hz Sk 













































Curve Fit 
Poles And Zeros 
POLES Oo ZEROS o 

1 -1. 1S5396k 182. 6938 
2 ~L1,. SBB4 & | 17D. 758 -149. S67 + | 366. 207 
3 -17.9099 +) 350. 331 275.096 +) 470.433 
4 -28.9193 +) S33. 241 -194.339 +) 745.371 
S -39.4818 +) 732.259 302.64 +) 874.078 
6 -56. 6217 +) 934. 359 -126.903 +J 1.14217k 
7 —-53. 4814 +) 1.14138k 281.275 +Jj 1. 2615ek 
8 -118.151 +J 1.31283 260.384 +) 1.6608 
9 -126.233 +) 1.SS83ik ~-le22.202 +) 1.77631 
10 -76.709 + J 1. 76231k 185.811 +) 2. 0O0981k 
11 -116.356 +} 1. 92914k 

Time delay= 0.0 S Gain= —2. 48K Scala= 1.0 











Figure E.12 PR-1592 direct measurement of Q, T=84.60°C. 
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APPENDIX F. PMMA DATA 


CONTENTS: 
-PMMA FIRST TORSIONAL MODE DATA (TABLE F.1) 


-PMMA FIRST FLEXURAL MODE DATA (TABLE F.2) 
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MS tae! A I B | € D i 
Label] SHEAR MODULU{TEMP (C) [RES FREQ |INPHASE | 
i 2.13E+9} -13.584| 2170.2295] __-0.186175 
2{ eof -12.428] 2167273] -0.186765 
sf 2e~of = 11 962] 2 165.391] 0.186045] 
at 2.12 +9 11.586] 2163.9711|  -0.185724 
sf Eo} = 11.153] 21624858] -0.18587| 
6] E-9f 10.682] 21.60.8468] -0.186084) 
7] te-9f 10.161] 2159.026 1] -0.185518) 
a] 2 0E-9f = 9.579] 257.0361] 0.185038) 
of 0E-9f -9.01 7] 2155.0013]  -0.184178) 
ro[ C2 10 +9} =8.395| _2152.9106] _-0.183241 
nf 09e-9f = 7.799] 2150.7247|-0.1823.44) 
rat 2 09E-9f 7.305] 21.48.7597] -0.181346] 
i3{ Eso} 6.611] 2146.365]  -0.179982) » 
ia[ 2.08E-9} -5.958 2144.058| _-0.179076 
is{f ss sdT 2. 08 E+ 9} -5.363] 21.41.8177] __-0.177692 
tel. 2 07E+9} -4.788] _2139.6827| _-0.176412 
i7[ 2.07E+9} -4.273| 2137.7027| _-0.175547] 
1af 2.06 E =] -3.593| 2135.4514 -0.17404 
1iof sd 2.06E-9 -2.972] _2133.1573| __-0.172888 
2c} | 2.05E*9f 2.54] 21313394] -0.17181) 
7h en 2.05E+9i -1.679| 2128.6048] — -0.17015} 
22, sd 2.04F-9 -1.01] 2126.0178| -0.168617 
23, | _2.04e-9| _-0.461| 21237212) _-0,167364 
24, id 2.04 - 91 0.206] 2121.2944] -0.16598! 


TABLE F.1 PMMA FIRST TORSIONAL MODE DATA 


7) ee 2.03E-9! 0.772] 2119.13.45 -0.16653 


26, sd 2.03E-9} 1.527] 2116.5967 
27, 2.02E-9 2.223} 2113.8699 


-0.16317 
-0.16154 


: 


28; sd 2.02E-9} 2.844] 2111.2689| -0.160346 


29[ TO E-Gf 3.487) 2108.8021]  -0.158949) 
sof TE 21] 2106.33.44] -0.157669 
3si[ 2 Eso! 7717 2103.8559|  -0.156204 
32]; sd 2.00E-9 5.385] 2301.4502| -0.154686 
33 1.99E+9 6.054] 2098.8933| -0.153449 
gat Cd Eso 718] 2096. 51491 -0.151969] 
35{ 1.98E+9 7.3671  2093.7681| -0.150656 
sof CT 9E-9f 8.052] 209.2058] -0.149201 
37] E+ 9f 8.703] 2088.6161)  -0.147574 
sof OTE} 9.567] 2085.9914]  -0.146263) 


399 E+ 9) 10.028] 2083.3422] _-0.144589 


ao] Cd EO 10.686]  2080.6743 


-0.143087 


7 1.95E+9 11.349] _2077.9868| _-0.141561 





a2{ SEO] 12.013] -2075.329'  -0.140178 
az[ Eso 2.665] 2072.6075]  -0.138532] 
44f 1.94E+9 13.322] _2069.8903| -0.137134 
a5[ Tt 93e-9f 15.984] 20671384 -0.13566) 
7 eel L.O3E-9f 14.644) 2064.3723| __-0.134306 
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TABLE F.1 CONTINUED 































= Label A i B ¢ D | 
47 1.92E-9} 15.3] 20615514] -0.132768! 
4[ id 1.92-9 15.948] 2056.7444] _-0.13157 
ag{_ CT EO} 16.615) 2055.9833} -0.130222 
50 2053.1806|  -0.128708 
st] 90E-9f 17.93] 20504226] -0.127139 
52 
sxs[ CT 9-9] 19.249] 204489761 -0.124718 
54 
55 : 1.88E-9 20.557] _2039.3161 
56 = 2036.5754)  -0.120365 
57 1.B7E-9 21.856] _2033.7549| -0.118723) 
58 1.87E-9 22.495) 2035.0002| _-0.117053 
59 1,.86E+9 23.161], 2028.0434] -0.115529 
60] CidEC( (CC BGEHD 23.808] 2025.1252|  -0.1 14124 
61 1.85E-9] 24.451 2022.308| _-0.112963! 
62 w= -1,.84E+9 25.104] _2019.4712| -0.112031 
63 1.84E+9 25.762! _2016.4934] _-0.110583] 
64 1.83€+9! 26.419] 2013.4772| -0.109642 
65 1.83E+9 27.08] _2010.4914| _-0.108044 
66 -1.82E+9 27.74] 20076932] _-0.105481 
67 1.82E*9 28.392 2004.729| -0.104024) 
68 1.B1E+9 29.045] 2001.7857| _-0.102619 
60, TB TE*F 29.697] 1998.8531)  -0.101044; 
70 1.80E+9 30.368 1995.884 -0.09968 
71 1.80E+9 31.015 1992.903| _-0.098511 
72 1.79E+9| 31.669, 1989.9731| -0.097088 
-0.095585] 
1.78E+9! 33.007] 1983.9916| _-0.094185 





-0.092622 


1.78E+9 33.659} _1980.9691 
34.328] _1977.9937| _-0.091032 
1.76€-9! 34.994] 1975.0318 -0.08946 


35.665 1972.0582 -0.087813 










79 L75E*9 36.32] 19691111] -0.086077 
gol 7 2E-of 36.974 1948.335]  -0.092582 
81 L.21e-9| 37.644] 1944.5879)  -0.091 188) 
82 172E-9 38.308] 19473617] _-0.088995 
83 L71E*9} 38.98] _1943.7924| _-0.087663 








84 1.70E+9t 39.642 1940.21] -0.086483 


as] 1.70E+9 40.317] _1936.5939| -0.085198 
66, sd 1.69E+9) 40.964] _1933.0179] _-0.083899 








87 1.68E-9 41.63] _1929.4515| _-0.082371 
gay C=C “C*‘CBE-OL 42.282] 925.8672] -0.081149 
89 
90 43.636] 19185008] -0.078538 
gif sd 1666-9] 44.284) 1914.7348] -0.077124} 
92 1.65E+9 44.953] _1911.0058| _-0.075723 





a ee 1.65E-9 45.611) 1907.1974)  -0.074261 
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3 | 
46.27 2| 
46.941| 
47.506} 






























48.269] 
48.9.4] 
49.607] 
50.251] 
50.895} 




























105 6 
106 SOE Ot 
lo7] 1.SSE-9I 
108 1.5 4E-9! 

Q 153+ 
110 1.S3e 
V4 1.S2ZE 
112 LSIE 
ae 1 "S1E-9i 
It4 SOE-9i 
11s — 1.50E-9 
116 1.49E+9) 
117 1.49E+91 
118 1.486+91 
119 1. 46E-Sj 
129 i. 46E-9: 
Tea 1 47E-S! 
22| | 47E-9 
123 eda | a7E-Q 
124 1,43€-9! 
ie i 


























53.616] « 
54.317} 
5 4.96| 


943 


60.861 
61.552] 
52.22 
62.879] 
63.906 
64.57 
55.313! 
65.979} 
66.646} 


















1804.9657 


TABLE F.1 CONTINUED 











1903.1464 
1899.0685 
1694.6504 

1820.032! 
1885.6676 
1881.5396] 
1877.5627 
1872.4947 
1869.3214 
13865.1634 
1861.0747 
1856.9836 
















18.48.7837 
1844.7561 
1840.9334 
1827.1454 
1333.4599 
1829.993 
1826.447 
LS2S. 01S 
1819.6695 
1816.8279 
1813.6706 
1$11.0656 
1806.6168 
1806.9164 
1805.5463 


1$05.4492 
1807.7255 


1813.45] 

















B) i 
-0.072794} 
-9.971298} 
-0.079089! 
-9.069288! 
-0.968618} 
-0.067665; 
-0.066566! 
-0.065518! 
-9.064159 
-9.0629031 
-0.061554! 
-0.06010S! 




















i 1352.9529 -9.058664l 


-9.057202| 

-0.05578! 
-0.0543 42) 
-0.052878! 
-90.0515 14 
-0.050028!} 
-0.04853 11 
-0.046935i 
-0.045358} 
-0.043866 4] 
-9.041856{ 
-0,040127| 
-0.036251| 
-0.076304] 
-9.034112{ 











-0.031776| 


-3.029136} 
-0.0260628} 


-0.0222559I 











1826.1035 -9.917304i 
1851.6973| -9.9114399} 
1894.1251| -0,0063043} 


1936.0094 





1964.325 
































19€1.0419 
1991.0099 
1996.4142 
2094.3418 

2909.291 
2915.7246 
2015.9262 
2921 4291 


2024.77293 






-9.0033287i 














-0.0020376) 
-0.0014411 
-0.0010852! 
-9.0008779i 
-0.000683 |} 
-0.000539 11 
-9.9094066i 
~9.9002596: 











-9.9900873} 


TABLE F.1 CONTINUED 


1.89E-9} 83.231] 2045.2737| _0.00072.43} 
isof 90-9 83.886] 2047.5409] 0.00081 98] 
TCC 90-9 82.851] 20476868] 0.000788 





[i A | 8 
Label [PMMA ——_—«[YOUNG'S MODUUTEMP (C) 


TABLE F.2 PMMA FIRST FLEXURAL MODE DATA 


Label 


| 


iC >) 
RES FREQ (Hz) |INPHASE VOLTA 


T['ST FLEXURAL | _6.001€-9 =14.252| _271.49919| __-0.47969 
2 [MODE DATA 6.073E-9| __-13.115| _271.09085| _-0.48556! 
| ae 


6.059E-9 12.561 270.77732 - 0.48373 


af 6048E-9} 12.163] 27054238] -0.48119 
s[ 6 058E-9f 11.794] 270.31047|  -0.47936) 
6f 6 028e-9f = 1.39] 270.09154l-0 47806] . 
7] 60 bE-9f = 10.83] 269.8718] -0.47493] 
rr ae -10.282| 269.54729 -0.47379 
of 5 990E-9f 9.6 45] 269.25162] -0.469717 


10 
VW 
12 
13 
14 
1S 
16 
is 
18 





5.979E-9] 9.102] 268.98 782 -0.4659 


a -8.509| 268.64069 -0.462 
| 5 9 49E +9 -7.876| 26830832 ~0.45769 
; 5.9336+9 -7.281| 267.95162 -0.45504 


5.919E-9} -6.703] _267.63059 -0.45082 
5.907E-9 -6.192| _267.37337 -0.44861 
5.892E-9} --5.615} — 267.02394 -0.44611 


| ststs—é~dT:=Ci‘“‘S;«C@#7QE-Of 4.996] 266.71982 -0.44244 


5.865E+9 4.387 266.40057 -0.43772 





i9f 5.848E+9 -3.702| _266.01735 -0.43413} 
20 pe ~0.42969 
21 =0.42585) 
22 5.800E-9 -1.674| 264.92506 -0.42164 
23 -1.092} 26456703 -0.41785 
24 
25 5.753E-9} 0.204| _263.86577 -0.4082) 
26 263.5001 -0.40427} 
27[ 5. 723E- 9) 1.506] _263.16049 -0.39807 


2a] 5 705E+9] 2 76 262.7626 -0.3942 
71 | 5.689E+9 2.848] 26238114)  -0.3899] 
so] 5.6 728+9 3.479] _ 261.98618 -0.38508 
Zif 55 7E-9f 17] 26164662] -0.38003 









PCT SS 908-9) 6.75] 26009807] 0.36183 


= 5.523E-9} 9.434| 258.53054 


0.34297 
5.507E+9 10.095] 258.14459 
a 10.76 


| 5.423£-9 13.407! _256.17371 -0.31646 
















32] sti‘(isé‘*Y:Sti<‘<it«‘S = 4.767} 261.25799 -0.37568 
33 5.623—-9 5.419, 260.86032 -0.37125} 
34 5.607E-9 6.08] _260.48888 





9 
5.574E-9 7.395] _259.72091 ~0.3567 
S.557E-9 8.102] 259.31066| __-0.35256 
5.540E-9] 8.766] —258.93155 -0.34804 



















257.373.49 
256.97245) 
256.56169 









-0.32122: 


5.406E-9f 14.087| 255.77406 -0.31469 


172 


TABLE F.2 CONTINUED 











MS abe A | 8 | c D | 
47, 5.389E-9} 14.74] 255.37819 -0.31019 
48|SC«dT:SsC«CS3 7 ZED 15.4] 254.95939 -0.30622} 
2 9-7) 16.088 25454521 -9.30175 
50. 5.55 7E-9 16.742] 254.1472) ~0.2983 
51 —— 5.320E+9 17.399} 253.73427 -0.29445 
52 5.303E-9 18.058] _253.31341] -0.29004) 
53, 5 285E-9 18.736) 252.88746| _-0.28631 
54 
ss[ dT 5.25089] 20.067] 252.06478| -0.2793 
56 5.234E+9 20.732 _251.66251 -0.27589 
57 5.216E+9 21.379| _251.25042 -0.27208 
58 5.200E-9 22.045} 250.65408| _-0.268172 
sof CT 5 183E+9] 22.692) 250.4315] -0.264384 
60 5.165E+9 23.379] _250.00432| _-0.260238 
61 5.1 486-9 24.032) _249.58872) -0.256662 






5.131€-9} 24.693 


26.017] 


249.1881! 
248.78806 
248.35937 
247.9417) 
247.51388 
247.00949 


246.313519 


~0.25 1666 


-0.224496 




















69 4.983E-9 29.3581 245.56749| _-0.227018) 
70 4.962E-9 30.006! 245.04262| _-0.227396 
“TI 4.942E-9| 30.68} 2445485| _-0.224925) 
72 -0.223304l 
73 24360181]  -0.221045 
_74[ dT 885E-9f 32.687] 243.15515] 0.219335 
75 
76 -0.21464) 
77| dT 829-9 34684) 241.7266] _-0.212782/ 
78 =0.210438} 
7ol——s—~dL*SCs=Ci«QOE* 9 36.035! _240.74994] _-0.210961 
60, sssid(ESC si BED 36.702} _240.41202| -0.197931 
81 4.76 4E-91 37.361] 240.09571! _-0.199636 
rs ee ee | 38.027] _239.26357| _-0.195996 
83 4.713E-9| 38.697;  238.82677| -0.196073 
ga[ | _4.605e-9] 39.388] _236.35605] _-0.195201] 
Es ee e471) 40.044] _237.86973|__-0.192627! 
86 4.656£-9 40.702| _237.36183{ -0.190512 
B7 41.375 236.879] _-0.187734) 
6af | S61 7E-OF 42.035 236.3735, 0.185289 
so[si*dL:C A SOGE-9} 42.708] _-235.89405]__-0.183189! 
oof —~dYSC“‘sK SS 7QE=Of 43.56] 235.40172]  -0.180755 
a1 4.560E-9 44.045] 23490971 -0 17835) 
92s 4.541E+9i 44.698} 234.427! -0.175398 
33 4.521E-9! 45.348] 233.915i12) -0.174064; 
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TABLE F.2 CONTINUED 


|_oM Label A | B z D 
o4{ sd 4.503E-9 46.018] _233.43427| -0.171791 
: en 4.482E-9 46.691} 232.69825| -0.169821 
96[ Cd a 6 SE-9} 47.367] 232.4034] -0.168203) 
g7] a aaae-9f 48.024] 231.90219[ -0.165549) 
gaff 48.685] 23137904] -0.163268) 
got 405-9} 49.356] 230.88204|  -0.161692 
yoo[ | 585E-9f 50.036] 230.36055] 0.159443] 
roif 36 5E-9f 50.701] 229.83951)  -0.157363] 
vo2f Cd 4 3a6E-9f 51.578] 22932451] -0.155281] 
ro3[ a 327-9] 52.063] 228.81621]  -0.153116 


104 4.354E-9] 52.744] 229.54222| _-0.112945 
105 4.339E-9 53.399] 229.15715] -0.110144) 


roof 4 323E-9f 54.053] 228.72323| -0.107914 
ro7[ Cd a 08E-9] 5.4.75] 228.53618] -0.1053811] 
108] i dT 2928-9 55.405] 227.89321| -0.103455] 
1 4.278E+9} 56.088] 22753613] -0.101111 
wo] CT 4 265E-9f 56.764] 227.17343|_ -0.09865 
mit 2s3e-of 57.431] 226.85622|  -0.096147) 
iat are-of 58.113] 22653992) -0.093714 
m3], 4 230€-9f 58.773] 226.2575] -0.091179) 
114 L_ 4.221E+9 59.443) 225.99864] _-0.088552 
115 4.212E-9| 60.137] _225.77141| _-0.085786 
nef dT 4 205E+9f 60.789] 225.59249) -0.083125 


ie 4.200E-9 61.465] _225.44929} -0.080141 
118 4.200E+9 62.142) 225.44588] -0.076889) 





11s 4.201E-+9 62.806! 225.48359 -0.073502 
120 4.210E-9} 65.037 225.70207 -0.069912} 


121) 4.224E+91 64.58] _226.09955| _-0.066015/ 
122] 4. 253-9 65.248] 226.85468] -0.061175, 
a3 297E-9f 65.885] 228.0447] -0.056175] 
124 -0.050956 
125[ sd 4 429-9) 67.1791 231.51122]  -0.046428 
iz6{ 4. 495E-9] 67.796] 233.21795] -0.042856) 
127, 4.567E-9 68.434] 235.10142) -0.039086 





128 4.652E-9} 69.0771 —237.26624) —-0.035826 
129] sd 4.735€-9| 69.736] _239.38204 -0.03804 





isof 4 7608-9} 70.421 239.99487] -0.0235041 
i31| sds 8 80-91 71.095| _243.01838] _-0.029325. 
i32[ a 96 6E-9f 71.803] 245,142.48] -0.027936] 
i33[ 4 998-9] 72.454] 2.45.93592]  -0.052924) 
isa] 87e-9f 73.101] 23815062] -0.036712 
3st ae e-ol 73.781] 237.67969] _ -0.035818) 
i3ef ae e-9f 74.45] 238.7571] -0.056.457 
i37{ 4 920E+9f 75.102] 243.99396]  -0.030403 
i38f_ 5.502E+9 75.801! _258.04527| -0.0250032 
139] TS SQNE-9 76.468! 260.1136] -0.0240697 
i4of sd 5.610E-9} 77-1661 260.5575| -0.0239374 
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TABLE F.2 CONTINUED 





) 









| iw Label A B Cc | 

i41[ 5621 E- 9} 77.742| _260.82056 
142], 563 0E 9 76.386) _261.02645] _-0.0233475 
rax[ 56 40e-9] 79.077] 261.240597_-0.0229913} 
14a 5.6508 +9} 79.784] 261.48127| -0.0229279) 
145[ TS 659€-9] 80.501] _261.68881] -0.0224509} 
ta6[ 5.66 96-9| 81.273) 261.93298]  -0.0223389) 
ta7[ 5679-9] 81.969) 262.14919] -0.0220627) 
i4e[ dT S56 88E+9) 82.66] 262.36215/ -0.0216973 
ragf 5 697E+9) 83.321] 262.55957]_ -0.02158.42) 
sof TS. 706E+9] 83.999] 262.77687/ -0.0211899 
[LS ISE*9] 82.698] _ 262.93856] -0.0213681) 
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APPENDIX G. POLYCARBONATE DATA 
CONTENTS: 
~ POLYCARBONATE FIRST TORSIONAL MODE DATA (TABLE G.1) 
- POLYCARBONATE FIRST FLEXURAL MODE DATA (TABLE G.2) 
- POLYCARBONATE FREQUENCY RESPONSE CURVES AND POLE ZERO 


INFORMATION FOR BOTH THE TORSIONAL AND FLEXURAL MODES OF 
THE HOLLOW SAMPLE ROD (FIGURES G.1 AND G.2) 
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TABLE G.1 POLYCARBONATE FIRST TORSIONAL MODE DATA 


z= Label | A | 8 | Cc ny) 
Label |POLYCT! SHEAR MODULUITEMP {C} IRES FREQ INPHASE | 


1ST TORSIONAL ; -14.516] 1225.37934j -0.133809) 






























3.695E-al -13.232| 122356535) -0.137884 
3-19-93 9.575E-al -12.671) _1222.2908| -0.140726 
9.661£-3| -12.221] 1221.36923} -0.141794] 
9.651£-3} -11.909] 1220.75875{ -0.143.439} 





9.640E-al -11.47{  1220.0566) 
9.626E-al -10.983} 1219.;5391 
9.609E-8l -10.435|  1213.12276 
9.S91E-8) -3.772| 1216.94905) 
9 S7SE-8) -9.2691 1215.95711 
9.558E-5] -8.672! 121484478 
9.54i£-4 -6.024) 1213.76545| -0.121538] 
9.525E-a} -7.472| 1212.78635| -0.121718} 
14 9.508E-8} -5 359] 121170271] -0.122435| 
9.497£-3i -5.291] 1210.72757 
16 9.47 76-3} ~5.741| 1209.70394)  -0.124491 
9.46 4E- a} -5.208] 1206.86628] -0.1236571 
9.4452-3) -4.4251 1207.67921 -0.123676 
19 9. 430E-a} -7.772| 1206.67449} -0.123072! 
9.41SE-3} -3.i41} 1205.71682| -0.135792} 
9.400E-8} -2.541] 1204.78427| -0.152431 
9.786E-a} -1.572! 1203.86879) -0.136115 
9.372E-a! -1.247| 1203.01206| -0.132107 
9.359E-al -0.614) 1202.17127 j 
25 9.247E-8) 0.03! 1201.24901 
26 9. 324E-8} 0.694] = 1200.5611 
9.322E- 8} 1.369] 1199.76499 
9.210£-al 2.036| 1199.01288 
9.298E-a} 2.667) 1198.24817 
9.288E- aj 1197.54846 






















-0.118362I 
-0.123588 




















































































-0.136224] 







































2.277E-af 1196.85019|  -0.1457031 
22 9.2562-6) 4645} 1196.14473|  -0.17946.4} 

9 255E-3I 1195.47763  -0.167006| 
34 9.2462-ai 5.964} 1194840591 -0.174313} 
35 0.276E-3/ 6.6351 119418978!  -9.164825} 
3 225€-3i 1193.52875| -0.154226} 

3.215E-i 1192.92552] 

o.296E-3l 1192.30088 

Ai 


1191.6868 
1191.08726 
10.671! 11904996] = -0.1597231 
11.306) 1189.55073{ -0.171644) 
11.9563 1139.2206 -0.172017] 
"1838590261 = -0.174337/ 
1P38607755}  -0.1731011 
1187 4640S| — -0 178579} 
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Label 





















TABLE G.1 CONTINUED 














A | B c ) | 
9.127€-af 14.623] 1186.89273} -0.1771981 
9.114€-8i 15.302] 1186.28628 -0.179307} 
9.103E-8I 15.961] 1185.600961 -0.181255} 
9.094E-al 16.6341 118499556] _-0.183068 
9.065E- aj 17.294] 1184.42276]  -0.189827} 
9.076E-at 17.957] 1182.3258i| -0.193584} 
9.066E-8} 18.615} 1183.20417 -0.198164 
9.0S7E-8} 19.286] 1182.59745} -0.2020471 
9.047E-3} 19.953] 1181.95599] -0.207766] 
9.027E-6i 20.617] = 1181.2102| -0.209218} 
9.023E-3ai 21.266] _1160.6884) -0.213297 
9.019E-3i 21.928] 1180.09464} -0.2162991 
9,000E-3} 23.284) 1178.3623| _-0.223699} 
8.991E-51 23.9171 1178.24662| -0.225883} 
3.981E-3) 24.572} 117764259] _-0.229738l 

7023248 4| 

6.9S6iE-Sf 25.915] 1176.21068 
8.95 1E-al 26.586] 1175.662305| -0.238062 
3.342£-3] 27.252] 1175.02057} -0.242211{ 
1174,.39662! -0.30239] 


1173.75079 
1173,09212 
1172.49072} 
1171.38196) 
1171.22765} 
1170.62065 
1169.86199 
1169.31279 
1168.61764 
1168.02641 
1167.3725 
1166.92272 
1166.27115 
1165.65729 
1165.01225 
1164.23137 
6| 1162.75581 
{  1167.1155 
1162.49039 
1161.95399 
116157054 
1160.!11262 








1160.7946 i 
1159. 46S5E5 
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-0.29066| 
-0.30681] 
















-0.30436| 
-0.31263} 
-0.718291 





-0.32166i 
-0.32309} 
-0.22599} 


Label 

































TABLE G.1 CONTINUED 


















A | 8 Cc 
3 653E-8] 45.984) 1155.9102 
8.621E-8} 46.558] 1154 47388 
8.626E-8i 47.218] 115407846 
8.687E-8} 47.999] 1158.17008 
8 648E-al 48.464) 1155.58259 
8.598E-3) 49.223{  1152.2272 
8.588E-al 49.989} 1151.57979 
8.580E-al 50.653} 1151.0558 
8.616E-3} Si.321)  1153.4352 
3.61 9E-8} 51.962] 115362245 


8.5896-3} 
8.501E-8l 


52.641) 
55.3521 


1151.61744 
1152.41556 








8.S90E-al 57.996] 1151.68542 
38.55 1E-al $4662] 1149.1969¢ 
8.575E-3i 55.329] 1150.67028} 


8.48 4£-8} 55.973] 114455866 








8.478E-6i 56.666] 1144.1764 
3.497E-8} 1145.4z%14 
8.531E-8l 57.979) 1147.73342 
3.542£-al 58.5643} 1148.50161 
8.520E-3} 59.209] 1146.97564 
8.49 1E-3) 59.963] 1145.04673 
8.425E-a} 60.636] 114058997 
3.426E-3! 61.293) 1141,272875 
8.41 7€-a) 61.945} 1129.77071 
8.418E-3 62.631] 1140.10566 
8.465E- 3) 63.28] 1143.28282 
8.299E-al 64.379] 113880927 
8.298E-al 65 078) 1138.74222 
8.254E-8] 65.593] 1135.74022 
8.41 4E-9} 66.258]  11239,8722 
3.406E-3i 66.999) 1139.21632 
3.2!2E-3} 57549] 1122.93717 
3 6.3 1127.55476 
§.342E-8] 66.974] 1134.97596 
3.758E-8) 59.612] 1176.01875 
8.337E-3} 70.289] 1136.02357 
é 70.9 1136.76147 
3a 1.567| 1132.70025 
2.272) 1132.77157 
"2.891! 1129.0i552 
3.255E-5 3.555] 1127.561396 
8.2 I5E-3/ 741907} 3127 $4296 
74.6 1125.35489 
6.8 
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| 0 | 
-0.71853) 


-0.39853] 
| -0.38375 
| -0.31983} 
-0.2189 
| -0.51946i 
-9.31738 
-0.31352) 





-0.313261 
-9.31027] 
{ -9.30518 
-9.30844) 
-9.30245] 
-9.353881 









-9.30119} 
-0.39213 
-0.38891} 
-0.38509} 
| -0.38467 
-9.38284 
| -0.38709 
-0.38414 
-0.27127 
-0.39798 
-0.28788 
-0.3702! 





{ -0.3767. 


-0.36997 
-0.36667 
~9.36842} 
-0.374491 
| -9.75204} 
-0.171 321 
-0.37352I 






-0.37227] 
-0.26628} 
| -0.27202 
-O 771611 
-9.36178} 














TABLE G.1 CONTINUED 


8.179E+8 77.48] 1123.77838 
Po 285-8] 78.133] 1130.93.476| 
> eon 8.160E-8 78.803] 1122.47441 
[C8 204E-8f 79.467) 1125.49742 


bY 8.146E+8 81.404] 1121.51693 
PCT 8.202 + 8 1125.37.46} 
8.1S7E+8 82.732] 1122.28236 
iso] 33E-8] 83.393] 112065684) : 

iif S5E-8f 82.427] 1122. 14604{ —-0.38.486) 
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TABLE G.2 POLYCARBONATE FIRST FLEXURAL MODE DATA 

































































































Es Label A i 5 i c i 
POLYC YOUNG'S (CORRTEMP (C) IRES FREQ INPHASE VOLTA 
1_JFLEXURAL 2.642E-9| 13.427) 127000391, -0.21429) 
2 [MODE DATA 2.637E-9] -12.328! 126.877084 -0.222156. 
3 [3-19-93 2.632E-9|  - 11.731) 126.773926] -0.228169 
4 2.629E+9 -11.291 126.696091| -0.231512 
5 2.627E-9 =10.936} 126639234) -0.233968| 
6 2.624E-9 -10.546! 126.584049] -0.236.452 
7 -10.132} 126.516773]  -0.237755 
8 2.618E+9 -9.648! 126.436356| _-0.227368| 
9 2.615E-9 -9.141) 126.353982) -0.229945 
10 2.611£-9 -8.563| 126.267921) _-0.233122 
1 2.608E-91 -8.046: 126.192124| _-0.236328] 
12 2.605E-9 -7.49} 126.107905| _-0.242368 
13 2.601E-9| -6.992! 126.026784| -0.246815 
i4f. sd 2 598-9} -6.45ii 125.951162) -0.249454 
18 -0.252972 
16 -0.257177 
17 2 587E-9} -4.632) 125685811} _-0.259309 
18 2.584E-+9} ~4.009} 125.596844) -0.264811 
i -3.327! 1255022991 -0.26675 
20, id S257 7-91 =2.792! 125 426931] _-0.267575 
21 2.573E-9! -2.157! 125.331196] -0.268595 
22] | 2 56 9F-9! -1.451)  125,.240551 ~0.27233 
23 2.565E+9i -0.791! 125.154687] -0.275119 
24 2.562E+9 -0.23) 125.068641| _-0.279032 
25[. sds 559-9! 0.514 124.984882| _-0.264202 
26 2.555£-9} 1.188: 124904864! -0.283489 
27 2.553E-9 1.794) 124642905} -0.286021) 
28 2.549E-9! 2.486! 124.75755} -0.288916! 
29 2.546E-9| F133! 124.679181| _-0.291323 
30 2.543E-9} 3.72) 124601762] -0.298152 
3) 2.53 98-9 4.421} 124517856 -0.30487! 
32 2.536E-9} 5.025! 124.439624 -0.30729 
| 124.368327 -0.30635} 

34 2530-9 6.348! 124.296501 -0.3112 
35 [_2.528E-91 6.993! 124.227743| 0.31594) 
36 2.525E-9) 7.624) 124.159476 -0.32253 
27, sid 2.522£-9} 6.314! 124.093972 -0.32685) 
38 2.520E-9) 8.996} 124.028582 -0.3315} 
39 2.517E-9) 9616: 123.964034 -0.33747 
40 2.514E-9} 10.273! 123.899811 -0.34123| 
ai) id ' 123.827324 -0.34563! 
42 2.509E-9} 11.56) _123.76906! -0.3511 
43 2.506E-9 12.2391 123.699074 -0.35691! 
44 2.504E-9} 12.886: 123.629647 -0.36307! 
45 2.501E-9 13,524! 123576923 -0.36493 









~0.37141 





46 2 499E-9| 14.2} 123.522434 
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|_| 
47 
48 
49 
50 
51 
52 
53 
54 
55 
56 
57 
58 
59 
60 
61 
62 
63 
64 
65 
66 
67 





81 
82 
83 
84 
85 
86 
87 
88 


89 


90 
91 
92 
93 





TABLE G.2 CONTINUED 
























122.047595 
121.932692 
121.823426 
121,921701 


32.792! 
33.462) 





i 





Label A { 5 i c | D i 
|S ststé‘ié‘LSCC*‘é«CAQGEE- QI 14.856: 123.458805 - 0.38072! 
ee ee 15.54: 127.40346 -0.38335] 
| ssid: 49 2E +9 16.181) 123.348461 -0.38767 
| sd 16.866 123.294596 - 0.39262} 

17.519) 123.235722 -0.39996 
cl 18.167! 123.178167 -0.39983 
2.483E+9 18.65) 123.115736 -0.40691 
19,492) 123.071968 -0.411 
20.176) 123.017078| _-0.41765 
[Cd 76-9 20.839) 122.9573831 0.4266 
2.475E-9] 21.509} 122.916744 -0.43011 
2.472E-9 22.139) 122666088] -0.4353) 
| id 2. 70 +9 22.823; 122.614112 -0.44088 
2468-9 23.478! 122.747125 -0.44543| 
2.466E-9 24.15} _ 122.703381 -0.45089 
2.464=+9 24.796, 122644444 -0.45669 
2.46 2E+9} 25.472) 122.596565 -0.46061 
2.459E-9 26.152! 122.539603 -0.46553 
(=. — 2 457E-9! 26.807) 122.47987| — -0.47061| 
2 455E-9) 27.463; 122.422573 0.47455 
2.452E-9 28.128) 122.365124 -0.48164] 

2.45 1£-9] 28.769) 122.325329 

2.447E-9} i 122.239238 

122.176972 

122.127921 


-0.48777 


-0.48064 
-0.48739 
-0.48912! 














2.441£-9 34.126} 122,079184 

2.440E- 9! 34.796! 122.055235 

ee ee ese 35.451) 121.777248 
2.418E-91 36.126: 121507604 

2.430E-9} 36.784) 121.601276 
Pd 2434-9] 7.428) 121.91.4386) 
2410E-9] 38.106) _121.300468) 

2.407E-9 38.772) 121.240164 

eae! 2.430E-9! 39.4171 121,808303 
40.082} 121.538649 
40.76) 121,084137 

2.414E-9! 41.442 121.401275 

2.415E-9} 42.085! 121.433903 

2.404E-9 42.781; 121.13976 

2.402E-9! 43.425) 121.111747 

ee er 44.063) 120.947438 
45.405} _120.788192 
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-0.48111 
~0.4793} 
-0.47127 
-0.477121 
~0.47637 
-0.47266 
-0.47439 
-0.48116 
-0.48265! 
- 0.48626) 
-0.48191 
~0.47932 
-0.4975 

~0.463 





TABLE G.2 CONTINUED 














| is Label A | B ; € D 

94 2.369F-9! 46.077) 120.763558 -0.501838} 
95 2.388E+9} 46744!  120.75862 -0.49943 
96 -0 49404 
97 2.395E-9 48.072} 120.924066 -0.47946 
98 2 400E+9 46.735! 121.058061 -0.4789 
99 2.394E+9} 49.392) 120.908332 ~0.48401 
100 2.388E+9 50.064! 120.756667 -0.48031 
101 2.378E+9 50.716} 120.490645 -0.48812 
to2y Cd 2 574E- 9] 51.376) 120.397957] =0.4964) 
103 2.373E+9 52.051} 120.375247 -0.48714 
1o4f 2.382E-9) 52.711} 120.609069| -0.4934 


105] t—~<“‘t«~‘*:*SC*‘ EDF! 52.365! 120.480052 -0.49199 


106{ sid 2B OE-9 54.037} 120.539759 -0.48619 








107 2.382E-9 54.697) 120.585249 -0.47684 
108 2.379F-9! 55.365) 120.516878 -0.46473 
109 2.379E-9 56.024! 120.522072 -0.46798 
110 3 i 56.688: 120451713 -0.49324 
nif 26 4E- 9} 57.3711 120.134821 -0.49446 
12 2.366E-9! 58.044) 120.182231 -0.49044! 
13 2.266€-9; 58.685: 120.196239 -0.48324 
14 2.366E-9} 59.345! 120.197276 -0.46414 








i 
1 

i 

] 

m5[ 2S 7E +9} 60.019, 119951637] -0.48158] 
116] —i«dEC ES E-O 60.711; 119.962119, _ -0.47913| 
] 
1 








7 

18 2.350E-9) 62.012? 119.785193 -0.47557 
119 2.347E-9 62.677: 119.70865 - 0.48729} 
120[ sd] sa 357-91 63.344 119.951269] _—-0. 488 
121 2.349E-9) 64.413) 119.764163 -0.46064 
122 2.355E-9| 65.087! 119.910461 ~0.45743 
123 2.336E-9! 65.742! 119.429209 -0.4912! 
124 2.326E-9i 66 408! 119.170168 -0.47085) 
125 2.349E-9} 67.052! 119.762901 -0.4527) 
126 -0.481 43 
2d at : 66.369! 119.399575 

126 -0.490.43 
129 -0.45481 
30] s—“‘~*:*C‘ AE“ 70.367; 118.853069 -0.48013} 
+z 2.328E-9| 71.022! 119.231173 -0.48014! 
132 2.235E-9} 71.691: 119.388544 -0.45177} 
133 2.31 2E-9| 72.243! 118.804162 -0.45707! 
134 2.305E-9} 73.014) 118643191 - 0.4669! 


135 : 72.628:  118.59607 
136 2.307E+9| 74.286} 118.681741 -0.47157| 


74.96) 118.932521 


75.638! 118615157 
140 2.3065-9! 76.92: 118.6485 -0.48032! 





























2.307E-9) 76.296} 118693026 
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FREQ RESP 


CURVE _ FIT 
10.0 





10.0 


dB 


dB 


























Curve Fit 
Pslas And Zeros 


POLES 20 


162. 04 -917.106 
10. 9253 1. 13209k+) 1. 23059 
-14.7513 1. O9516k -1.28672k+) 2. 56731k 
-22. 0523 2. 2045 1. 14602k+)j 3. 28775k 
-28. 3637 3. 32774k -1.25748k+) 4. 88251k 
-39. 0655 4. 45974k 964. 421 5. 4483k 
-45. 9193 5S. 50757 -944. 807 7.21364k 
—S5. 976 6. 74894k 677.097 7.67721k 
-65.0111 7. 88372k 49.5549 7. G$4962k 
07 So. s003 7. 950S4k -38. 6579 S. 82967 
117-96. 0727 9. O3Sk 


L 
2 
3 
4 
s 
6 
7 
8 
i=} 
1 


a a a a a 


Time dalay= 0.0 S Gain= 23. 86 Scala= 





Figure G.1 Polycarbonate frequency response curves and 
pole zero information for the torsional 
mode of the hollow sample at T=21.7°C. 
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FREQ RESP 
CURVE _ FIT 
0.0 


0.0 


dB 


dB 












































Curve Pit 


Poles And Zeros 
POLES S ZEROS 8 
1 -1.49901 +j 132. 666 —590.074 
2 -2.73361 +j 362.252 34.6517 
BS -5.15614 +j 699.516 353. 353 
4 7. 55234k 
5 -94.3374 +j 535.34 
6 64.2728 +} 562.902 














Time delay= 0.0 S Gain=-SOS. 3p Scale 1.0 





Figure G.2 Polycarbonate frequency response curves and 


pole zero information for the flexural 
mode of the hollow sample at T=21.8°C. 
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APPENDIX H. RLC CIRCUIT DATA 
CONTENTS: 
- FREQUENCY RESPONSE CURVES AND POLE ZERO INFORMATION FOR 


THE RLC CIRCUIT AT VARIOUS RESISTANCE LEVELS (FIGURES H.1- 
H.7) 
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Figure H.1 Frequency response curves and pole zero 
information for the RLC circuit (R=21M). 
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Figure H.2 Frequency response curves and pole zero 
information for the RLC circuit (R=71M). 
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Figure H.3 Frequency response curves and pole zero 
information for the RLC circuit (R=121M). 
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Figure H.4 Frequency response curves and pole zero 
information for the RLC circuit (R=171M). 
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Figure H.5 Frequency response curves and pole zero 
information for the RLC circuit (R=221N). 
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Pigure H.6 Frequency response curves and pole zero 
information for the RLC circuit (R=261N). 
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Figure H.7 Frequency response curves and pole zero 
information for the RLC circuit (R=371N). 
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